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PREFACE 


The 1988 Johnson Space Center (JSC) National Aeronautics and Space Administration (NASA) / 
American Society for Engineering Education (ASEE) Summer Faculty Fellowship Program was 
conducted by the University of Houston and JSC. The 10-week program was operated under the auspices 
of the ASEE. The program at JSC, as well as the programs at other NASA Centers, was funded by the 
Office of University Affairs, NASA Headquarters, Washington, D.C. The objectives of the program, 
which began in 1965 at JSC and in 1964 nationally, are 

1 . To further the professional knowledge of qualified engineering and science faculty members 

2. To stimulate an exchange of ideas between participants and NASA 

3. To enrich and refresh the research and teaching activities of participant’s institutions 

4. To contribute to the research objectives of the NASA Centers 

Each faculty fellow spent at least 10 weeks at JSC engaged in a research project commensurate with 
his/her interests and background and worked in collaboration with a NASA / JSC colleague. This 
document is a compilation of the final reports on the research projects done by the faculty fellows during 
the summer of 1988. Volume 1 contains reports 1 through 14, and volume 2 contains reports 
15 through 26. 


CONTENTS 

Volume 1 

1. Ahmed, Selina: "Comparison of Soviet and U.S. Space Food 

& Nutrition Programs” 1-1 

2. Allen, Robert: "Production System Chunking in SOAR: Case Studies 

in Automated Learning" 2-1 

3. Bartram, Peter N.: "Development of a Software Interface for Optical 

Disk Archival Storage for a New Life Sciences Flight Experiments 

Computer” 3-1 

4. Bell, Christopher and Lachman, Roy: "Computer Technologies and 

Institutional Memory” 4-1 

5. Berry, Frederick C. : "Development of Parallel Algorithms for Electrical 

Power Management in Space Applications” 5-1 

6. Botros, Nazeih M.: "Automatic Voice Recognition Using Traditional 

and Artificial Neural Network Approaches” 6-1 

7. Casserly, Dennis M.: "Development of an Atmospheric Monitoring Plan 

for Space Station” 7-1 

8. Chang, Ing: "Model Formulation of N on-Equilibrium Gas Radiation 

for Hypersonic Flight Vehicles" 8-1 

9. Davis, John E.: "A Model for Plasma Volume Changes During Short 

Duration Spaceflight” 9-1 

10. de Korvin, Andre: "An Evidential Approach to Problem Solving 

When a Large Number of Knowledge Systems is Available" 10-1 

1 1 . Elam, Reid P. : "Effects on Motor U nit Potentiation and Ground 

Reaction Force From Treadmill Exercise” 11-1 

12. Flynn, George J.: "An Assessment of the Micrometeoritic Component 

in the Martian Soil” 12-1 

13. Geer, Richard D.: "Evaluation of Available Analytical Techniques 

for Monitoring the Quality of Space Station Potable Water 13-1 

14. Kelly, Frederick A.: "Model Evaluation, Recommendation, and 

Prioritizing of Future Work for the Manipulator Emulator Testbed” 14-1 

Volume 2 

15. Knopp, Jerome: "Optical Joint Correlation Using the Deformable Mirror Device” 15-1 

16. McLauchlan, Robert A.: "Intelligent Control of Robotic Arm /Hand Systems 

for the NASA EVA Retriever Using Neural Networks” 16-1 


PRECEDING PASS W.ANK NOT FIT, VIED 


iii 



17. Nechay, Bohdan R.: "Maladjustment of Kidneys to Microgravity: Design 

of Measures to Reduce the Loss of Calcium” 17-1 

18. Nelson, Dennis O.: "Evaluation of Crustal Recycling During the Evolution 

of Archean- Age Matache wan Basaltic Magmas” 18-1 

19. Nelson, Kerri L.: "Lithospheric Buckling and Intra-Arc Stresses: 

A Mechanism for Arc Segmentation” 19-1 

20. Nerheim, Rosalee: "Toward an Image Compression Algorithm for the 

High-Resolution Electronic Still Camera” 20-1 

21. Odell, Patrick L.: "Equating an Expert System to a Classifier in Order to 

Evaluate the Expert System” 21-1 

22. Randhawa, Manjit S. : "Feasibility of U sing High Temperature Superconducting 

Magnets and Conventional Magnetic Loop Antennas to Attract or Repel 

Objects at the Space Station” 22-1 

23. Starks, Scott A.: "Visual Perception and Grasping for the Extravehicular 

Activity Robot” 23-1 

24. Tan, Arjun: "Estimation of Velocity Perturbations in Satellite 

Fragmentation Events” 24-1 

25. Tolliver, C.L.: "Time-Domain Imaging” 25-1 

26. Willey, Ronald J.: "Arc Jet Diagnostics Tests” 26-1 


Pfr 


PA$!5 WLANE NOT FE/tfED 


v 



N89-20074 


OPTICAL JOINT CORRELATION USING THE 
DEFORMABLE MIRROR DEVICE 

Final Report 

NASA/ASEE Summer Faculty Fellowship Program 
Johnson Space Center 


Prepared by: 

Academic Rank: 

University and Department: 


NASA/JSC 

Directorate: 

Division: 

Branch : 

Contract Number: 


Jerome Knopp 

Associate Professor 

University of Missouri- 
Columbi a/Kansas City 
Department of Electrical 
and Computer Engineering 
Truman Campus 
600 West Mechanic 
Independence, MO 64050 


Engineering 

Tracking and Communications 
Tracking Techniques 
NGT 44-005-803 


15-1 



ABSTRACT 


An experimental investigation of the Deformable Mirror 
Device (DMD) developed by Texas Instruments at Dallas for 
use in object identification was completed. The DMD was 
tested as a joint correlator. The DMD was used as a spatial 
light modulator on which the squared modulus of the Fourier 
transform of test object pairs was written. The squared 
modulus was phase encoded on the DMD after it had been 
thresholded and rewritten as a binary phase function. The 
thresholding was found to produce a sharp peak in the 
autocorrelation when the test objects were matched and no 
significant peak in the case of distinct objects. It was 
concluded that the use of the DMD as a joint correlator 
looks promising and further studies should be carried out. 
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INTRODUCTION 


* 


The location and recognition of objects belonging to 
small well defined sets ( e.g., handtools drifting in space) 
or objects belonging to large less well defined classes 
(e.g., boulders on the surface of Mars.) are typical parts 
of mission requirements on several NASA programs including 
the planned unmanned expedition to Mars, the EVA (Extra 
Vehicular Activity) retriever program and autonomous 
rendezvous and docking. In most cases the speed of the 
recognition process must rival that of the human eye-brain. 
This is difficult competition. Even the fastest available 
digital computers would not be able to process data fast 
enough to keep pace with human vision. An interesting 
alternative to digital computing is optical computing. 
Optical computing takes advantage of massive parallel 
processing. Within the last 10 years, optical processing 
techniques have matured to the point that they are worthy of 
consideration in tasks where speed is a premium. Although it 
lacks the flexibility of digital processing at the present, 
it can perform certain fixed tasks, such as Fourier 
transforms in two dimensions, almost instantaneously (i.e., 
essentially in the length of time it takes the light to pass 
through the optical processing system.) At present, the 
Tracking and Communications group at NASA is studying the 
use of a hybrid system that exploits the best of both 
digital and optical processing for artificial vision 
problems. The key element in this approach is the Deformable 
Mirror Device (DMD) , a new type of integrated optics device 
under development by Texas Instruments (TI) in Dallas. 

The DMD consists of a 1/4" square array containing 128 
x 128 mirror elements. Each mirror element consists of a 
set of four individual cantilevered mirrors arranged as in a 
2X2 matrix (roughly in the form of a cloverleaf.) Each of 
the mirror elements can be deflected as a group using the 
electrostatic forces determined by capacitatively stored 
charge between a mirror and the array back plane. The charge 
stored, and thus the deflection, is determined by on-chip 
electronics and addressing. The DMD is. in essence, a 
deformable mirror surface whose shape can be digitally 
controlled. The mirror is also fast. Each element can be 
deflected at an 8 kHz. rate, although the current 
experimental setup operates the DMD at 50 Hz. to maintain 
compatibility with a 50 Hz. European standard video. 

At the present, the DMD is being studied for use as an 
optical correlator. Correlation strength measures can be a 
key feature in vision identification tasks [2]. 

Correlation, however, has not been normally used in 
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practical vision systems since it is computationally 
intensive and slow. This is not the case using the DMD where 
correlation using optical processing is attractive. The two 
classical approaches to correlation using optical processing 
are Vender Lugt Filtering (VLF) and Joint Correlation (JC) . 
The use of VLF has been one of the prime thrusts of previous 
research on the DMD and is described elsewhere [3] . The JC 
approach was developed and studied this summer as a part of 
research conducted under the ASEE summer faculty program. 

The key results of this work are described here. 

The basic approach to JC involves two Fourier transform 
operations and one square law (intensity detection) 
operation. First two objects to be correlated are optically 
Fourier transformed in the same operation. This is done by 
placing the objects side by side in front of a lens and 
transil luminating them with coherent light. The lens 
produces a joint Fourier transform of the objects in terms 
of the electric field 15] . The field is then detected with a 
square law detector. The squaring operation produces 
crossterms containing the Fourier transform of the joint 
correlation of the objects. An additional Fourier transform 
of the squared field (using another lens) will yield two 
crosscorrelation terms. Optically this results in 
autocorrelation peaks in the form of intense points of light 
when the objects are matched. The degree of correlation as 
well as the details of the object structure determine the 
"sharpness" of the correlation point. The implementation of 
the JC process on the DMD required some modifications of the 
traditional JC process. In conventional optical 
implementations the square law operation is usually carried 
out by detecting the intensity of the joint Fourier 
transform on photographic film. A positive transparency of 
the transform then serves as an input for the second Fourier 
transform operation using a lens that yields the JC. In the 
approach used here the film is replaced by the DMD. The DMD 
is not primarily an amplitude modulating device; it is for 
the most part, a phase modulating device. In this case, the 
intensity was phase encoded. 

In previous research, using a computer simulation [5], 
it was shown that directly writing the intensity as a phase 
would yield good quality correlations. This was attempted 
initially by writing the squared Fourier transform patterns, 
calculated by a computer, directly on the DMD. The DMD was 
then illuminated with collimated light from a helium neon 
laser. The reflected light was then Fourier transformed with 
a 50 cm. lens. The resulting inverse transform was then 
imaged on a CCD array camera and observed on a monitor. 

These first attempts were negative and no correlation peaks 
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were observed. The process was reviewed for flaws. One of 
the major problems in the previous computer simulation is 
the assumption that the DMD could be accurately modeled 
using simple reflecting surfaces whose mechanical deflection 
was directly proportional to the charge pattern written into 
the DMD. In the real device, this does not appear to be the 
case. Based on data obtained by TI [61, it appears that 
there is a nonlinear relation between charge patterns 
written on the DMD and the resulting pattern of mechanical 
deflection. Therefore, an alternative approach was tried 
that appeared to be much more robust. 

The key identification feature in the NASA tasks is a 
sharp autocorrelation peak that shows strong correlation 
between an object to be identified and a reference object. 

A sharp peak is due primarily to a strong single frequency 
sinusoidal component occurring in the joint Fourier 
transform. The strength of this single frequency component 
is proportional to the number of equally spaced point pairs 
on the objects being correlated. If this sinusoidal 
component can be emphasized in the case of a good match, 
then the correlation spot will be bright. A simple approach 
that appears to accomplish this, and completely avoids the 
nonlinearity of the DMD, is to threshold the modulus squared 
of the Fourier transform and rewrite the squared transform 
as a binary object. If the sinusoidal component is present, 
it will be emphasized as the fundamental component of a 
square grating. Typical results using this approach are 
described next. 


FINDINGS 

The use of JC on the DMD was tried using thresholding. 
The Fourier transform of test object pairs were calculated, 
squared and then thresholded by replacing intensities above 
a constant value (the threshold) by a 1 and intensities 
below the threshold by 0. In the cases studied, the < 
threshold was set at one percent of the peak value. The 
resulting pattern was written on the DMD. Then an inverse 
Fourier transform was carried out optically and the 
resulting "correlations" were recorded. Quotes are used 
here to indicate that correlation is used loosely since a 
nonlinear mapping has replaced the joint Fourier transform 
intensity with a binary object. 

Figs. 1 and 2 show two sets of objects that were used 
to test the DMD as a recognition correlator. Fig. 1 shows 
the pair SS and Fig. 2 shows the pair ST. In this case, the 
S and T can be considered as reference objects that are 
being compared with a test object, i.e., another S. A sharp 
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Fig. 1. Tha tast pair SS , 
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Fig. 2. The t*st pair ST. 
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correlation peak is desired when the objects are matched to 
indicate recognition. However when the objects are not 
matched then a weaker peak (or a peak too weak to observe) 
is desired. 

Figs. 3 and 4 show the computer generated Fourier 
transforms after threshholding and converting to binary 
images for the SS and ST pairs, respectively. Note the 
strong single spatial frequency component for the SS 
transform in the diagonal direction. This frequency 
corresponds to the separation between the objects and will 
determine where the correlation peak will appear in the 
output plane. Figs. 5 and 6 show respectively, the 
resulting correlations made using the DMD. The images were 
taken from a CCD array camera and stored using a frame 
grabber. Note the presence of strong off-axis correlation 
peaks in the case of a match as shown in Fig. 5 and the lack 
of distinct peaks in Fig. 6. It should be pointed out that 
all the figures shown are only crudely grey scaled using a 
10 level grey scale plotter in conjunction with a laser 
printer. The grey scale code is shown below the figures. 

CONCLUSIONS 

The results presented suggest that the DMD can be used 
to correlate images for identification. with good results 
provided that proper thresholding is used. Since only a 
small set of objects was tried it is not possible to offer a 
general scheme for thresholding. It is not known how well 
the particular threshold scheme used here will work in 
general. Furthermore, more testing is needed to evaluate 
scaling, rotation and noise effects using a more interesting 
object set. 

Also, the similarity between the thresholding approach 
used here and in optical neural networks has not gone 
unnoticed. Indeed, the only element missing to convert this 
system to an optical neural net system is a feedback 
mechanism. This is being considered in future studies using 
a system the author is presently developing (7] . 

ADDITIONAL COMMENTS 

During the preparation of this final report, but after 
the author's ASEE summer epoch, some further results were 
obtained on the JC during a joint "overnight research 
adventure at NASA JSC " involving Richard Juday and a summer 
intern student from The University of New Hampshire ,Joe 
Bailey. The DMD was tested using photographic 
transparencies as the correlator input. The squared Fourier 
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Fig. 3. The computer generated Fourier 

transform of SS after thresholding. 
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Fig. 4. The computer generated Fourier 

transform of ST after thresholding 
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Fig. 5. The joint correlation of the 
SS pair. 
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Fig. 


. Th« joint eorralation of tha 
pair ST . 
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transform was obtained optically using a lens and directly 
recorded on a CCD camera. The camera image was written into 
the DMD using a frame grabber. Good quality correlations 
were obtained. These results will be reported in detail at a 
later date. It would appear that real-time JC is possible 
provided a fast spatial light modulator (SLM) is used as an 
input to the DMD correlator. Such a real-time device is in 
the process of being installed. The only additional 
requirement is a fast frame grabbing routine for 
transferring the Fourier transform intensity to the DMD. 

This should not be a problem at the 50 Hz. frame rate now 
being used to write to the DMD. 
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ABSTRACT 


Adaptive/general learning algorithms using varying neural network 
models are considered for the intelligent control of robotic arm 
plus dextrous hand/manipulator systems. Results are summarized and 
discussed for the use of the Bar to /Sut ton/ Anderson neuronlike* un- 
supervised learning controller as applied to the stabilization of an 
inverted pendulum on a cart system. Recommendations are made for the 
application of the controller and a kinematic analysis for trajec- 
tory planning to simple object retrieval (chase/approach and cap- 
ture/grasp) scenarios in two dimensions. 


i 


Over v i ew 


INTRODUCTION 


The research work reported herein is important to the future 
development of the NASA/JSC EVA Retriever. This highly autonomous* 
f ree --flying robot or robotic system is comprised of MMU * arm and 
smart hands. It is being developed to aid crewmen in the perfor- 
mance of EVA tasks including the chase* capture and return capabili- 
ty required for adrift crewmen or station equipment. The ultimate 
goal of the work in developing this system is to enhance the effec- 
tiveness of EVA crewmen II, 233 . 

The intelligent control of robotic ar m/hand systems using 
neural network learning controllers is very relevant to EVA Retrie- 
ver development. This follows because of the need for autonomous* 
adaptive behavior in both planned and unplanned contexts in the 
space environment. Neural networks and related advanced learning 
controllers offer such capabilities C233. 

The work reported herein is concerned with the i nvest i ga t i on 
and development of neural networks or other types of advanced 
learning controllers as: 

(a) Supervised controllers with training which because of 
their connective, associative memory structure can 
develop significant controller generalization capabili- 
ty. Such generalization can lead to similar performance 
of the retriever ar m/hand controller in different but 
analogous physical system situations and in stochastical- 
ly related loading/excitation environments. 

<b) Unsupervised controllers which can self train/adapt to 
new learning situations and also exhibit significant 
gener al i zat ion capability. As learning develops* and un- 
familiar situations become familiar ones, these neural 
networks should provide feedforward compensation with 
less compensation via the feedback path 17, 11, 15, 263. 
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Neural Networks for Intelligent Control 

Neural networks are massively parallel, distributed processing 
systems. They have the ability to continuously improve their perfor- 
mance via dynamic learning C7, 9, 15-18, 363. As used in this re- 
port, neural networks refers to "artif icial" , i.e., programmable 

systems of processing elements. As such they form a research area 
of intense interest in artificial intelligence. 

Initial neural network research concentrated on the computa- 
tionally intensive areas of adaptive signal processing, as, e.g., 
pattern recognition, real-time speech recognition and image inter- 
pretation. Recently there has been a resurgence of interest in 
neural networks because of (a) Advances in training algorithms for 
networks, and (b) Availability of extremely fast, relatively inex- 
pensive computers for implementing these algorithms. These develop- 
ments have lead to the consideration of neural networks for the 
real-time i den t i f i c a t i on and control of large flexible/articulated 
aerospace and robotic systems [7, 87, 883. 

Neural networks can provide mechanisms for (a) Associative 
memory, (b) Pattern recognition, and <c) Abstraction. These are 
emergent properties of networks of neuronlike units with adaptive 
synaptic connections CIO, 1 4- , 88, 89, 383. These mechanisms arise 
from the neural network being a system of interconnected "neuron- 
like” elements modeled after the human brain. This system operates 
on input data in an "all at once” mode rather than in a conventional 
computer’s "step by step” algorithmic approach [7, 9, 893. Differ- 
ent learning architectures can be used in training for intelligent 
control. This is done to provide appropriate inputs to the system so 
that the desired responses are obtained. Uncertainty and noise can 
be handled by a neural network via the Hebbian type of associative 
learning arising from adaptively modified connection strengths [21, 
893. Kawato el al C15-183 indicate that a neural network model can 
be used to control voluntary movement with applications to robotics. 
Implemented as a multilayered, h i er arch i c a 1 1 y intelligent control 
system, neural networks can be implemented to effect the following: 

(a) Pattern recognition/ condition matching 

(b) Trajectory and approach, grasping, etc. opertation 

<c) "Point of view” transformations - as, e.g., visual to 

sensor/end effector to object, etc. 

(d) System (robot, object, etc.) state observer or model 
synthesis and simulation behavior 

(e) Generation of mo t i on/ac tua tor commands. 

Adaptive control is useful for systems which perform over the 
large ranges of uncer t a i nt i es which result from large variations in 
physical and operating parameter values, environmental conditions, 
and signal inputs. However, adaptive control as such (i.e., without 
unsupervised 1 ear n i ng /unant i c i pa ted problem solving features) has 
difficulty with the fallowing generic problems in designing 
controllers: 

* Sensor data overload - arising from (a) Data redundancy 
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per se* and (b) Specialized, rarely required data 

* Nu 1 t i -spec tr al , multi-sensor data fusion and mapping/use 
in the proper feedback control law 

* Need for system robustness to handle large parameter 
excursi ons 

* Required h i gh~speed , real-time control degradation resul- 
ting from time consuming artificial intelligence calcula- 
tions 

* Unsolved sensor choice and placement problems for 
r obo t i c / 1 ar ge control systems. 

It should be noted that human intervention is used in traditional 
control systems operating with large uncertainty. Such interven- 
tion is unacceptable in many real-time applications. This is espec- 
ially true for the hostile space environment in which the NASA EVA 
Retriever is to operate Cl, 233. It means that automatic techniques 
for handling uncertainty must be developed. Neural networks show 
great promise for the intelligent, unsupervised control of the mul- 
tiple arm plus dextrous robotic hands of the Retriever. The next 
section of the report describes the author’s research work with the 
Barto et al intelligent controller which is a special kind of neural 
network with associative search and associative critic neuronlike 
e 1 events . 

ACE/ASE NEURONLIKE LEARNING CONTROLLER 

The Bar to/Sutton/Anderson adaptive learning controller is 
composed of two types of neuronlike elements with significant 
unsuper’v i sed problem-solving capacities. These elements are the 
associative search element ( ASE ) and adaptive critic element < ACE > . 
Barto et al 1983 used a single element of each type. Their ASE 
element exhibits a learning strategy which is similar to the earlier 
"BOXES" adaptive problem solving system of Michie and Chambers 
C29, 253. The ASE/ACE elements embody refinements discussed in the 
literature by Barto and colleagues £2-6, 30-313. They evolved from 
the heterostatic brain function and adaptive systems work of Klopf 
£19, 203. Adding a single ACE element improves the learning perfor- 

mance over that of a single ASE alone. This can be clearly shown by 
comparing the problem-solving capabilities of BOXES with those of a 
single ASE/single ACE learning system and solving the control prob- 
lem of balancing an inverted pendulum on a cart. It is interesting 
to note that strong analogies exist between the behavorial inter- 
pretations of the ASE, ACE adaptive elements and animal behavior in 
instrumental learning. There are also strong parallels with the 
"bootstrap adaption" systems work of Widrow et al £33-353. This 
work considered the (a) pun i sh / rewar d critical learning and (b) pat- 
tern recognizing control problems. Relevant artificial (i.e., 
programmable) neural networks the ASE, ACE neuronlike elements are 
significant. This follows because they indicate that if adaptive 
elements are to learn effectively as network components, then they 
are constrained to have adaptive capabilities at least as robust as 
these Barto et al learning controller elements £23. 

Figure 1 depicts the inverted pendulum on a cart system which is 
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Figure 1. Representative Model for Cart and Inverted Pendulum 
System . 
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Figure 2. ASE and ACE Controller for Cart Plus Inverted Pendulum 
System 
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to be controlled. Here the cart can move within the bounds 
indicated on a one-dimensional track. The pendulum can move only in 
the vertical plane of the cart and the track. The applied force 
F ^ t ) results from the output of the learning controller. It is 
applied in a bang-bang <+/-) manner and acts with a fixed magnitude 
to the left or right at discrete time intervals. The pendulum-cart 
system is described by a four state variable model in the time 
domain C83. The four state variables are as follows: (a) y(c. - the 
position of the cart on the track, (b)^p - the angle of the 
pendulum with the vertical, (c) X'c - the cart velocity, and (d) 

0p - the rate of change of the pendulum angular displacement. The 
state variable model for this system can be written as 



CMp*Ly*Co 
Opy'y CJp 




5(0 P )J 

+Mp*l^3 


EMp*ly>si n < & p ) * 
|C-Mp*g*y*sin< © 


^c*sgn< Vj») 
p> - Mp*u ) P : 


+ Fapp 




( 1 ) 


Physical parameters in the above equations specify pendulum length 
and mass, cart mass, the coefficients of friction between the cart 
and the track and at the pin connection between the pendulum and the 
cart, the applied control force, the force due to gravity, and time. 
Table 1 defines the notation used in equation 1. 

The system of first order equations has been solved using second 
order numerical intergration procedures which have been implemented 
in the FORTRAN computer program NRLNET. In implementing the 
learning controller algorithm the state space has been partitioned 
based on the following 252 quantization interval thresholds: 

< 1 ) 


( 2 ) 


(3) 

i <4> 

Figure 2 depicts the ASE plus ACE adaptive learning controller 
of Barto et al C23. The neuronlike learning system can be described 
; by the following equations: 

| Element output v(t) which is determined from the decoded stats 

1 quantization interval vector input 

i y = f L ( w ( I , t ) * x ( I , t ) ) + n < t ) 3 (2) 

j 

i 
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Cart position ■tfc : +/- 0.8, +/- 2.4 m ,(4 quantiza- 
tion intervals including failed regions above and 
below 2.4m) 

Pendulum angular displacement : 0,+/- 1, +/- 6, 

+/- 12 degrees, (7 quantization intervals including 
failed regions above and below 12 degrees) 

Cart velocity yc : +/- 0.5, + /-0^m/s, (3 quantiza- 
tion intervals) * 

Pendulum angular velocity ^p: +/- 50, + /- OO degrees 
per second, <3 quantization intervals) 
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Here the noise n(t) is a real random variable with probability 
function p < x ) and f is either a threshold, sigmoid, or identity 
transfer function. For the work reported herein, p(x) is the zero 
mean Gaussian distribution with standard deviation O* , and f is 
the bang-bang type thr esho 1 d func t i on ; 

f(x) = + 1 , x .EG. 0 (applied force action 

to the r ight ) ( 3 ) 

-1, x .LT. 0 (applied force action 
to the left) 

ASE weights w(I,t) , 1 .LT. I .LT. N which change over discrete time 

as f o 1 lows : 


w(I,t+l) = w(I,t) + ALPHA *- r ( t ) * e(I,t) 


(4) 


In equation <4: 
ALPHA 

r ( t ) 


positive constant determining the rate of 
change in w ( I , t ) 

real-valued reinforcement at time t 


e(I»t) = eligibility at time t via the input pathway I. 


Eligibility traces for the ASE weights which exponentially decay 
with increasing time, given in equation 5 as: 

' e ( I , t+ 1 ) = DELTA * e(I,t) + (1-DELTA) * y(t) * x<I,t> (5) 


in which. 


DELTA = the eligibility decay rate. 


ACE weights v(I,t) . 1 . LT . I . LT . N which change over discrete t i me 

as foil ows : 


v( I , t+1 ) = v(I,t> + BETA * rhat(t) * xbar ( I , t ) (6) 

In equation 6, 

BETA = positive constant defining the rate of change 

of v < I , t ) 

rhat(t) = r(t) + GAMMA * p(t) - p(t-l), the improved 
internal r e i nf or cement signal for the critic 
e 1 ement 


xbar ( I , t ) 


p ( t ) 


= LAMBDA * xbar(I,t) + (1-LAMBDA) * x(I,t), 
the eligibility traces for the ACE weights 

= I , t ) * y(I,t), the prediction of eventual 

reinforcement 
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GAMMA 


reinforcement learning rate 


LAMBDA = xbar<I,t) trace delay weight 

Barto and Sutton C2, 5] explain the derivation of the ACE learning 
rule as used above. Additional discussion of the ASE , ACE adaptive 
learning controller can be found in references C30, 313. 

ASE/ACE LEARNING CONTROLLER RESULTS 

This section of the report discusses representative results 
obtained by the author with his FORTRAN computer program NRLNET131 
implementing the Barto et al ASE/ACE neuronlike learning control- 
ler. This program is the result of several modifications by the 
author to incorporate general data file input and the file and prin- 
ter plot output of the applied control force and the four state 
variables as functions of time. The original FORTRAN program 
NRLNETOO was the author’s implementation of a PASCAL program written 
in 1988 by Doug Walker of GHG in support of the Special Projects 
Branch ( EC5 ) in the Crew and Thermal Systems Division at NASA/JSC. 


TABLE 1. SUMMARY OF PHYSICAL PARAMETER VALUES FOR CART PLUS INVERTED 
PENDULUM SYSTEM 


Me 

Mp 

Lp 

MUc 

MUp 

Fapp 


Cart Mass, 1.0 kg 
Pendulum Mass, 0.10 kg 
Pendulum Length ( 0.50 m 

Cart Coef f i cient of Coulomb Friction, 0. 005 
Pendulum/Cart Pin Coefficient of Friction, 0.00002 
N m sec/rad 

Magnitude of Force Applied to Cart in x Direction, 
(+/-) ION 


TABLE 2. SUMMARY OF THE ASE/ACE NEURONLIKE LEARNING CONTROLLER 
PARAMETERS 


ALPHA 

BETA 

DELTA 

GAMMA 

LAMBDA 

M 


<r 


Rate Constant for ASE Weights? 1000.0 
Rate Constant for ACE Weights? 0.50 
Decay Rate for ASE Eligibility Traces? 0.90 
Learning Rate for Improved Internal Reinforce- 
ment ? 0 . 95 

Decay Rate for ACE Eligibility Traces, 0.95 
Mean Value for Gaussian Normal Distribution Used 
to Define ASE/ACE Output Noise Function? 0.00 and 
0. 10 

Standard Deviation Value for Gaussian Normal Dis- 
tribution Used to Define ASE/ACE Output Noise 
Function? 0.01? 0.05? 0.10? 0.15? 0 . £0 ? and 0 . £5 


Table 1 gives the physical and control parameter values used in 
the simulation work with NRLNET131 for the cart plus inverted pendu- 
lum system depicted in Figure 1. Values used for the ASE/ACE neuron- 
like learning controller parameters are summarized in Table £. These 
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Figure 3. Example Simulation Re- 
sults Showing Learning 
Performance of ASE/ACE 
Learning Controller 
System 
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Figure 4. Average Number of Trials 
for Five Runs as a Func- 
tion of Standard Devia- 
tion With Mean Value as 
Par ameter 
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Figure 5. Applied Force Fapp(t)> N 
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Figure 8. Pendulum Angular Dis- Figure 9. Pendulem £ngular 

p 1 acement t ) » r ad Velocity < t ) » rad /sec 
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parameter values were used to generate the simulation performance 
results plotted in Figures 3 and 4. Figure 3 plots curves for the 
number of time steps until failure versus the trial number. These 
are typical curves for individual runs of the ASE/ACE learning con- 
troller system. Figure 4 gives plots of the average number of trials 
which are required to stabilize the cart plus pendulum system for 
5000 time steps (100 seconds with t = 0.02 seconds). The average 
number of trials are given as a function of the standard deviation 
for the Gaussian normal random noise process with the mean for the 
process as parameter. The mean equal zero curve indicates a trend 
toward an increasing number of trials as the standard deviation is 
increased from 0.01 to 0.25. The other curve for the mean equal to 
0.10 shows relative constancy over the same range in standard devia- 
tion. These runs were originally made to examine the sensitivity of 
the ASE/ACE learning controller performance to variation in the 
noise process used in generating its output function. An additional 
objective was to develop a base from which the generalization and 
robustness properties of the controller weights could be investi- 
gated. Five runs were used to generate each point plotted in Figure 
4. The results shown in Figures 3 and 4 are in general agreement 
with those published by Barto et al C23. However, the author has 
found that his NRLNET131 implementation of the ASE/ACE controller 
usually takes a lesser number of trials for successfully learning to 
stabilize the cart plus pendulum system for both the 5000 ^ t (100 
seconds) cases shown here and the 200,000 4 t (66.7 minutes) cases 
which the author ran to directly compare his results with those of 
Barto et al. Extensive runs were not made for the 200,000 A t (66.7 
minutes) stabilization period because of the excessively long 
elapsed* time rerequired for the VAX system available to the author 
to return answers for a single run. 

Figure 5 shows the controller output force which is applied to 
the cart in stabilizing the inverted pendulum. Here the applied 
force is plotted as a function as a function of time over the first 
100 4 1 intervals (2 seconds). Extensive runs have been made with the 
ASE/ACE controller system and all exhibit the characteristic + /- 10N 
on— off or bang— bang behavior with At = 0.02 sec. This value of the 
time increment should be adequate, based on physical system oscilla- 
tion behavior, for the second order numerical integration scheme 
used . 

Figures 6-9 plot the four state variables: cart displacement 
, cart velocity ijCj * pendulum angular displacement <<9) , and 
pendulum angular velocity , respec t i vel y . They are also plotted 

as functions of time over the first lOO^At intervals (2 seconds). 
Consideration of these and similar time domain results for the state 
variables and the applied force indicates that (a) significant in- 
efficiencies can occur with respect to the input force and its im- 
pact on the actual state variable behavior of the cart plus pendulum 
system, (b) with^ t = 0.02 sec there may be some interaction be- 
tween the numerical integration method used and the dynamics of the 
ASE/ACE learning controller. To investigate (b) above, additional 
runs were made in which Al t was reduced ( At = 0.01 , 0.005, 0.001 
sec). These results although not included here did show significant 
reduction with decreasing At in the bang-bang nature of the input 
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force and the higher frequency oscillations present in the state 
behavior over time (especially for the cart linear and the pendulum 
angular velocities) . 

The author has extended the single ASE/single ACE learning con- 
troller system to include two search and two critic elements. The 
elements in each pair work in parallel. Since the outputs are aver- 
aged in the 2 ASE/2 ACE learning controller system, it has -Fapp , 0, 

+Fapp as three possible outputs. This extension was implemented in 
the author '5 FORTRAN computer program NRLNET20. Initial runs indi- 
cate that the new controller as implemented has good performance up 
to a maximum learning point (maximum time for stability as a func- 
tion of number of trials). Beyond this point the learning is severe- 
ly degraded with increasing trials, or a form of limit cycle behav- 
ior occurs. These results indicate that the split-decision nature of 
the 2 ASE/2 ACE system as implemented in its averaging form may 
cause the observed behavior. In this case using a 3, 5? etc. ( i .e. ? 
odd number of elements) in the ASE/ACE system may be warranted. 

These controllers would also have a “smoother" (i.e., less hang- 
bang) control action. Another alternative to improve performance is 
to more richly connect the elements both within and across the . 
search element and the critic element layers. This would give the 
ASE/ACE neuronlike controller system a counter pr opagat i on/Gr ossberg 
layer plus Kohonen layer type of neural network structure C12, 133. 

CONCLUSIONS 

An examination has been made of the use of neural networks for 
the intelligent control of robotic arm-plus hand /man i pu 1 ator sys- 
tems for the EVA Retriever. Based largely to the present time on a 
review of the literature and computer simulation work, this examina- 
tion has indicated that a hierarchical , multi-layer neural network 
system can be used for intelligent control. Baseline feedforward 
control is used in conjunction with trajectory planning in these 
systems. Joint torque feedback provides the correction signal. These 
systems have the character ist ic that as additional response behav- 
iors are learned, much of the control action passes to the feedfor- 
ward path. 

Additional investigation into neural networks for intelligent 
control has focused on the use and extension of the Barta et al 
neuronlike ASE/ACE intelligent controller. A FORTRAN family of com- 
puter programs ( NRLNETXX ) were developed by the author as extensions 
of a previous Pascal language implementation of the controller at 
NASA/JSC. Work with these programs has concentrated on the follow- 
ing: (a) Verifying published results for convergence to stable solu- 

tion (number of trials for a specified period of stability), (b) 
Developing graphics, etc. feedback tools to monitor system behavior 
(as, e.g., given by the applied control force and the four state 
variables as functions of time), (c) Investigate learning control 
behavior as a function of the number of unsupervised trials required 
to obtain stability and the random process parameters (Gaussian pro- 
cess mean and standard deviation), and <d) Basic extensions to the 
learning controller network i ncorpora t i ng two adaptive search ele- 
ments (ASEs) and 2 adaptive critic elements (ACEs) in its structure. 
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RECOMMENDATIONS 


This section of the report presents recommendat i ons for the 
intelligent control of smart robotic arm plus hand systems using 
neural networks. These recommendations are based on the results 
presented above and on additional related work done by the author 
during Summer 1988. They are presented in the form of a research 
and development program plan. The R & D program plan gives activ- 
ities that can continue the author’s research begun during the 1988 
summer program. 

(1) Investigation of two dimensional graphics as a kinematic 
si mul a t ion tool for planning EVA object retrieval in terms 
of the approach to and grasping of objects using an artic- 
ulated two-link arm/scissor hand system. 

(2) Implementation of dynamics* sensing* and control models of 
the articulated two-link arm/scissor hand system. It is 
desired to mount this arm/hand system on a cart to repre- 
sent the EVA Retriever in two dimensions. 

( 3 ) Examination of hierarchical neural networks with fuzzy 

logic reasoning as adap t i ve/gener a 1 learning systems com- 
prised of (a) Network ar ch i tec tures * (b) Transfer func- 

tions* and (c) Dynamic learning rules. These systems can 
employ joint torque and state vector feedback to control 
the arm/hand system(s) in object retrieval as discussed 
above . 

(A) Investigation of extensions to the Bar to /Anderson neuron- 
like learning system and counter propagation/back propaga- 
tion type networks to the related problem of stabilizing/ 
controlling the motion of simple and compound (articulated 
linkage) pendulums on a cart. Successful employment here 
can lead to similar use with the arm/hand retriever sys- 
tems. 

REFERENCES 

1. ’’Extravehicular Activity (EVA) Retriever Program Plan*" NASA/ 
Johnson Space Center* Houston* Texas* Report No. JSC-22199* 

May 1987. 

2. A.G. Barto* R.S. Sutton* C.W. Anderson* “Neuronlike Adaptive 

Elements That Can Solve Difficult Learning Control Problems*" 
IEEE Transactions on Systems* Man* and Cybernetics * Vo 1 . SMC- 
13, 1983* pp 839-896. 

3. A.G. Barto, C.W. Anderson* R.S. Sutton* “Synthesis of Non- 

linear Control Surfaces by a Layered Associative Search Net- 
work*" Biological Cybernetics * Vo 1 . 93* 1982* pp 175-185. 

9. A.G. Barto* R.S. Sutton, “Landmark Learning: An Illustration 
of Associative Search*" Biological Cybernetics * Vo 1 . ^2, 1981* 

pp 1-8. 

5. A.G. Barto, R.S. Sutton, “Simulation of Anticipatory Responses 
in Classical Conditioning by a Neuronlike Adaptive Element*" 
Behavioral Brain Research * Vol. 9, 1982* pp 221-235. 

. A.G. Barto* R.S. Sutton* P.S. Brouwer, “Associative Search 


OF PQ-K QUALITY 


6 



I 

I 


Network: A Reinforcement Learning Associative Memory," 

Biological Cybernetics , Vol, 4Q , 1981 » pp £01-211. 

7. Behnam Bavarian, "Introduction to Neural Networks for Intelli- 
gent Control," IEEE Control Systems Magazine , Vol. 8, No, 2, 
April 1988, pp 3-7. 

8. Robert H. Cannon, Jr., Dynamics of Physical Systems , McGraw- 
Hill, New York, 1967, pp 703-710. 

9. M. Caudill, "Neural Networks Primer: Parts I, II, and III," 

AI Exper t , Part I, Dec. 1987, pp 46-52, Part II, Feb. 1988, pp 
55-61, Part III, June 1988, pp 53-59. 

10. M. Fujita, "Adaptive Filter Model of the Cerebellum," Biologi- 
cal Cybernetics , Vol. 45, 1982, pp 195-206. 

11. Allon Guez, James L. Eilbert, Moshe Kam, "Neural Network Archi- 
tecture for Control," IEEE Control Systems Magazine , Vol. 8, 

No. 2, April 1988, pp 22—25. 

12. Robert Hech t-Ni e 1 sen . "Counterpropagation Networks," Proceed- 
inos, IEEE International Conference on Neural Networks , 1987. 

13. Robert Hech t-N i e 1 sen , "Applications of Coun ter propag a t i on Net- 
works," Neural Networks, Volume 1 , 1988. 

14. J.J. Hopfield, "Neural Networks and' Physical Systems With 
Emergent Collective Computational Abilities," Proceedings of 
the National Academy of Sciences , Vol. 79, April 1982, pp 
2554-2558. 

15. Mitsuo Kawato, Yoji Uno , Michiaki Isobe, Ryo j i Suzuki, 

"Hierarch ical Neural Network Model for Voluntary Movement 
With Application to Robotics," IEEE Control Systems Magazine , 
Vol. 8, No. 2, April 1988, pp 8-16. 

16. Mitsuo Kawato, "Adaptation and Learning in Control of Volun- 
tary Movement in the Central Nervous System," Advanced Robot- 
ics , Vol. 2 (In Press). 

17. M. Kawato, M. Etoh, Y. Oda, N. Tsukahara, "A New Algorithm for 

Voltage Clamp Iteration: A Learning Control of a Nonlinear 
Neuronal System," Biological Cybernetics , Vol. 53, 1985, pp 

57-66. 

18. M. Kawato, M. Isobe, Y. Maeda, R. Suzuki, "Coordinate Trans- 
formation and Learning control for Visually-Guided Voluntary 
Movement With Iteration: A Newton-Like Method in Function 
Space," Biological Cybernetics (In Press). 

19. A.H. Klopf, "Brain Function and Adaptive Systems - A Hetero- 
static Theory," Air force Cambridge Research Laboratory 
Research Report AFCRL-72-0164 , Bedford, Massachusetts, 1972. 

(A Summary Appears in Proceedings of the International Con- 
ference on Systems, Man, and Cybernetics , 1974) . 

20. A.H. Klopf, The Hedonistic Neuron: A Theory of Memory , 

Learning, and Intelligence , Hemisphere, Washington, D.C., 

1982. 

£1. Daniel S. Levine, "Neural Networks: Growth Field for the 1990s," 
SIAM News , Vol. 21, No. 4, July 1988, pp 8-9. 

22. Richard P. Lippmann, "An Introduction to Computing With Neural 
Nets," IEEE A5SP Magazine , IEEE Acoustics, Speech, and Signal 
Processing Society, Volume 4, Number 2, April 1987, pp 4-22. 

23. R.A. McLauch lan , "Dynamics, Control and Sensor Issues Pertinent 
to Robotic Hands for the EVA Retriever System," Participant’s 

16-14 


OF d^ AL PAGE 'S 

0F P0 °R quality 



24 . 


25. 


26. 

27. 

28. 

29. 

30. 
31 . 

32. 

33. 

34. 


35. 

36 . 


ORIGINAL 

OF POOR 


PAGE IS 
QUALITY 


Final Report - 1987 MASA/ASEE Summer Faculty Research Program, 
Department of Civil and Mechanical Engineering, Texas A & I 
University, 3 September 1987. 

D. Michie, R.A. Chambers, "BOXES: An Experiment in Adaptive 
Control," in Machine Intelligence 2 , E. Dale and D. Michie, 
Editors. Oliver and Boyd, Edinburgh, 1968, pp 137-152. 

D. Michie, R.A. Chambers, "BOXES as a Model of Pattern-Forma- 
tion," in Toward a Theoretical Biology, Vo 1 . 1, Prolegomena , 

C. H. Waddington, Editor. Edinburgh University Press, Edin- 
burgh, 1968, pp 206—212. 

W.T. Miller, et al , "Application of a General Learning Algorithm 
to the Control of Robotic Manipulators," International Journal 
of Robotics Research , Vo 1 . 6, No. 2, Summer 1987, pp 84-98. 
Demetri Psaltis, Athanasios Sideris, Alan A. Yamamura, "Neural 
Controllers," Proceedings IEEE First Conference on Neural Net- 
works , San Diego, California, Vol. 4, June 1987, pp 551-558. 
Demetri Psaltis, Athanasios Sideris, Alan A. Yamamura, "A 
Multilayered Neural Network Controller," IEEE Control Systems 
Magaz i ne , Vol. 8, No. 2, April 1988, pp 17-21. 

D. E. Rumelhart, J.L. McClelland, Parallel Distributed Proces- 
sing: Explorations in the Microstructure of Cognition, Volume s 
I and I I , MIT Press, Cambridge, Massachusetts, 1986. 

R.S. Sutton, A.G. Barto, "Toward a Modern Theory of Adaptive 
Networks: Expectation and Prediction," Psychology Review , Vol. 
88, 1981, pp 135-171. 

R.S. Sutton, A.G. Barto, "An Adaptive Network That Constructs 
and Uses an Internal Model of Its World," Cognition and Brain 
Theory , Vol. 4, 1981, pp 213-246. 

David W. Tank, John J. Hopfield, "Collective Computation in 
Neuronlike Circuits," Scientific American , Volume 256, Number 
12, December 1987, pp 104-114. 

B. Widrow, M.E. Hoff, "Adaptive Switching Circuits," in 1 960 
WESCQN Convention Record , Part IV, 1960, pp 96-104. 

B. Widrow, N.K. Gupta, S. Maitra, "Punish/Reward: Learning 
With a Critic in Adaptive Threshold Systems," IEEE Transac- 
tions on Systems, Man, and Cybernetics , Vol. SMC-3, 1472, 
pp 455-465. 

B. Widrow, F.W. Smith, "Pattern-Recognizing Control Systems," 
in Computer and Information Sciences , J.T. Tow, R.H. Wilcox, 
Editors. Clever Hume Press, 1964, pp 288-317. 

B. Widrow, J.M. McCool, M.G. Larimore, and C.R. Johnson, 
"Stationary and Nonstationary Learning Char ac ter i st i cs of the 
LMS Adaptive Filter," Proceedings of the IEEE , Vol. 64, 1976, 

pp 1151-1162. 


16-15 


N89 - 2007 6 

^ t 

MALADJUSTMENT OF KIDNEYS TO MICROGRAVITY: 

DESIGN OF MEASURES TO REDUCE THE LOSS OF CALCIUM 

Final Report 

NASA/ASEE Summer Faculty Fellowship Program— 1988 
Johnson Space Center 

Bohdan R. Nechay, D.V.M. 

Professor 

University of Texas Medical Branch 
Department of Pharmacology & Toxicology 
Galveston, Texas 77550 

Space and Life Sciences 
Medical Sciences 
Biomedical Laboratories 
Donald M. Foster, Ph.D. 

October 25, 1988 
NGT 44-005-803 


Prepared By: 

Academic Rank: 
University & Department: 

NASA/JSC 

Directorate: 

Division: 

Branch: 

JSC Colleague: 

Date Submitted: 

Contract Number 


17 -1 



ABSTRACT 


Losses of skeletal calcium and body fluids occur during prolonged exposure to 
microgravity. The kidney plays a major role in regulating the physiological functions 
involved. Relative to this regulatory function, the kidney performs three operations: 
filtration of blood plasma through the glomeruli, reabsorption, and secretion of fluid and 
electrolytes so that needed components are retained and only waste is eliminated in the 
urine. 


Using data published in "Biomedical Results from Skylab" (1), we performed new 
calculations that reflect more directly the operations of the kidney in the handling of 
calcium, sodium, chloride, potassium and phosphate during space flight. These 
calculations revealed that the fraction of filtered calcium that was rejected by renal 
tubules and excreted in the urine increased by 71%, from 1.77% (preflight) to 3.02% (in- 
flight) of the filtered load. This represents a large absolute increase because the total 
filtered amount is huge. Because the tubular rejection fraction of other ions increased 
relatively less than that of calcium, we postulate the inflight development of a specific 
renal defect that causes an excessive loss of calcium in urine and thereby contributes to 
the weakening of bones. 

Because the reabsorption of glomerular filtrate (the conservation of body fluids and 
electrolytes by the kidney) depends upon the energy of adenosinetriphosphate (ATP) and 
enzymes that liberate its energy, we propose that an induction of the major renal enzyme 
system involved, sodium and potassium-dependent adenosinetriphosphatase, may improve 
the reabsorption of sodium directly and calcium indirectly and thereby lead to improved 
hydration and reduced calcium loss. 
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INTRODUCTION 


Because of the known detrimental shifts of fluids and electrolytes, particularly 
calcium losses, that occur under conditions of microgravity, we undertook an evaluation 
of renal function studies during the Skylab flights and of renal handling of important 

ions. Using data published in "Biomedical Results from Skylab" (1), we performed new 
calculations that reflect more directly the operations of the kidneys of astronauts under 
conditions of microgravity. 

To make this report clear to the nonspecialist, we now give a brief summary of 
pertinent renal physiology, in part derived from Pitts (2). The two kidneys weigh about 
300 grams and thus constitute 0.4% of total body weight. In humans, each kidney is 
composed of about one million units called nephrons, which function in parallel. Each 
nephron begins with a glomerulus (microscopic sieve) attached to a tubule. Fluid filtered 
by the glomeruli passes through the renal tubules, to collecting ducts and eventually to 
the urinary bladder. The kidneys receive the greatest blood flow, in proportion to 
weight, of any organ of the body. Renal blood flow is equal to 20-25% of cardiac output 

or about 2000 liters per day. Urine formation starts with the filtration of blood plasma 

through the glomeruli at the rate of about 200 liters per day. The renal tubules 

reabsorb needed components of the glomerular filtrate, which contain most constituents 

of blood plasma excluding protein, and add to the filtrate undesirable products by 

secretion. The net result is excretion of waste in urine at an average rate of two liters 
per day. Thus, about 99% of glomerular filtrate is reabsorbed. Given the fact that an 
average adult has 40 liters of fluid in the body, the internal autodialysis of all body 

fluids through the kidneys takes place five times every 24 hours (200 liters of glomerular 

filtrate/40-5). The net result of renal filtration, reabsorption and secretion is the 
regulation of composition, within narrow limits compatible with life, of body fluids, often 
referred to as the constant internal environment. 


METHODS 

Glomerular filtration rates (GFR) and renal handling of ions were estimated (2,3) 
based on plasma and urinary concentrations obtained from published data from the nine 
crewmen who participated in Skylab flights (Tables 1 and 2). Where appropriate, plasma 
concentrations of substances obtained at different times were averaged to match the 
periods of urinary excretion data available. 

The renal clearance of a substance is defined as the number of milliliters of plasma 
completely cleared of that substance in one minute. If the substance is completely 

filterable through glomeruli, that is, it is neither reabsorbed nor secreted, if it has no 
effect on renal function, and if it can be accurately measured in plasma and urine, then 
its clearance can be taken as a valid measure of GFR. Normal GFR per 1.73 square 
meters of body surface area, standard for an adult, averages 125 ml/min in men and 110 
ml/min in women. In healthy persons, GFR is remarkably stable from day to day over a 
period of years. 

The formula for calculation of renal clearance is U x V/P, where U and P represent 
concentrations of the substance in urine and plasma (mg/ml) respectively, and V 
represents urine formation rate in ml/min. If the clearance of a given substance is 
larger than the simultaneous GFR, that substance is said to be "secreted" by the tubules 
of the kidney. If clearance of a substance is smaller than GFR, the substance is said to 
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be "reabsorbed.” It should be noted that, in the absence of more detailed studies, the 
clearance of a substance is a net figure and does not take into account simultaneous 

secretion and reabsorption that might also be taking place. 

The standard method of determination of GFR is by measurement of inulin clearance 
following intravenous infusion. Inulin is a high (5,000) molecular weight polysaccharide 
that is completely filterable and is neither secreted nor absorbed by the kidney. 
Obviously, this method cannot be used during space flight. Instead, measurement of 

endogenous creatinine clearance, which agrees reasonably well with the inulin clearance, 
is usually taken as an approximate measure of GFR and was used in our calculations. 
Creatinine is a product of muscle creatine phosphate; its production is constant and is 

proportional to muscle mass. Renal elimination of creatinine is also constant. The 

agreement of GFR as measured by inulin and by creatinine clearances results from the 
balance of two errors: 1) the analysis method slightly overestimates the creatinine 

concentration in plasma (but not in urine), and this gives the appearance of a lower 
creatinine clearance; and 2) the true clearance of creatinine is slightly higher than GFR 
because of tubular secretion. These are the limitations in using creatinine clearance as 

an approximate measure of GFR in these studies. 

Filtration rate of an ion (mg/min) through glomeruli was calculated by multiplying its 
filterable concentration in plasma (mg/ml) by the GFR (ml/min). Because 40% of plasma 
calcium is bound to plasma proteins, which are too large to pass through the glomeruli, 
the filterable fraction of calcium was taken as 60% of total plasma calcium. Phosphate 
was assumed 90% filterable. Virtually all plasma sodium, potassium and chloride are 
freely filterable. 

The tubular reabsorption rate (mg/min) of an ion is its filtration rate (mg/min) minus 
its excretion rate (mg/min). Tubular rejection fraction or fractional excretion (%) was 

calculated .by the formula: Excretion rate (mg/min)/ filtration rate (mg/min) x 100. This 
yields the percentage of the filtered quantity of a given ion that escaped tubular 
reabsorption. 


RESULTS 

Renal function data calculated from Tables 1 and 2 are shown in Table 3. Creatinine 
clearances proved remarkably stable throughout the missions, indicating intact renal 
vascular autoregulatory mechanisms that control the GFR in the face of the well known 
circulatory disturbances and body fluid shifts that occur during space flight (4-6). In 
flight, GFR’s held within 2-3% of the preflight values and fell only a modest 9% during 
the 60-85 day flight period. 

Tubular reabsorption of individual ions was depressed during flight and reverted 
postflight to baseline, preflight values (Fig. 1). Peak fractional excretion rose 71% (from 
1.77% to 3.02%) for calcium, 12% for phosphate, 38% for chloride and 40% for sodium. 
This indicates a generalized increased tubular rejection of several ionic components of 
the glomerular filtrate with a relative selectivity for calcium. 

The data for renal handling of calcium are plotted in Figure 2. The curves for 
fractional excretion and urinary excretion rate of calcium were similar the former 
exhibited a smoother natural progression, suggesting that this was the mechanism for 
increased calcium loss in the urine. Increased fractional excretion coincided with a 
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slight but significant rise in plasma calcium concentration (peak 6%). Filtered and 
reabsorbed calcium increased (peak 8-9%) during flight and paralleled the creatinine 
clearance curve throughout. Thus, the increased fractional excretion of calcium during 
flight was the dominant finding. 

To see more clearly the mechanism of tubular transport of calcium, in Figure 3 we 
have plotted renal calcium handling adjusted per 100 ml of GFR. The curve for the 
urinary calcium excretion rate resembled more closely the curve for fractional calcium 
excretion than that shown in Figure 2. The adjusted filtration and reabsorption rates of 
calcium remained parallel but became more stable than those shown in Figure 2, and peak 
increases during flight were only 4-5%. This is further evidence that the calcium 
excreted in urine appears there as a result of diminished tubular reabsorption. 

The urinary excretion of 0.29 gm of calcium per day during the Skylab flights was 
nearly twice that of control values of 0.16 gm per day. Thus, 7.8 to 11.05 gm of calcium 
were lost during 59 day and 84 day Skylab flights respectively. These amounts are equal 
to more than the total of extraskeletal intracellular and extracellular calcium stores. 


DISCUSSION 

The outstanding finding of our study is that renal tubular reabsorption of calcium 
abruptly diminished during space flight and rapidly returned to control levels after 
landing. This resulted in urinary excretion of 0.29 gm of calcium per day during flight, 
up from 0.16 gm per day control rate. The net total calcium loss attributable to the 

kidney during the Skylab flights of 59 to 84 days equaled 7.8 gm to 11.05 gm 

respectively. The crewmen had an adequate average dietary intake of 0.73 gm calcium 
per day but went into negative calcium balance because fecal calcium excretion also 
increased (1,7). The calcium -wasting by the kidney is further demonstrated by the lack 

of renal conservation of calcium in the face of the increased fecal calcium losses. 

To place the magnitude of the calcium loss in perspective, it is essential to realize 
that all extracellular body fluids contain only about 1 gm of calcium, cells contain 5 gm, 
and bone contains a total of 1,200 gm (3). Thus, any sizeable calcium loss, like that 
which occurred during the Skylab flights, must come eventually from bone. Calcium 
homeostasis involves continuous mineral turnover in bone, absorption from food and 
reabsorption by the kidney. Overall calcium metabolism is primarily set to maintain a 

constant calcium concentration in extracellular fluid, even at the expense of bone 

density. Absorption from the gut is controlled by vitamin D, bone resorption and renal 
reabsorption are promoted by parathyroid hormone and inhibited by calcitonin (3,8). 

Based on the above considerations, options for prevention of calcium loss are: 1) 
increase intake and/or decrease losses of calcium from the gut; 2) decrease calcium loss 

from bone; and 3) decrease calcium loss from the kidney. Relative to option 1, stable 
and adequate calcium intake during Skylab flights did not prevent negative calcium 
balance (1,7). In fact, there was an increased amount of calcium lost in feces (7). 
Enhancement of calcium absorption from food deserves further study. 

Decreased egress of calcium from bone appears to be under most intensive 

investigation at present as a solution to bone demineralization during space flight. In 

our opinion, producing a halt to bone resorption without increased absorption from the 
gut and reduced loss via the kidneys would lead to a reduced intracellular and 
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extracellular calcium concentration. Since calcium is required for muscular contraction 
and nerve conduction, inadequate calcium levels may lead to increased weakness and 
atrophy, already a problem during space flight (9). Diminished supply of calcium to the 
heart could lead to impairment of cardiac output and to irregularities of heart rhythm 
(arrhythmias) that could result in sudden death. 

Relative to the third option, reversal of the acquired defect in tubular reabsorption 
of calcium would be a desirable solution to the excessive calcium loss during space flight. 
Because renal calcium reabsorption is in part dependent upon the transport of sodium, we 
propose that an induction of renal sodium and potassium-dependent 
adenosinetriphosphatase (Na and K ATPase), an enzyme system that allows the use of 
adenosine triphosphate (ATP) for ion transport (10), would improve tubular reabsorption 
of sodium directly and that of calcium secondarily. In addition to calcium conservation, 
this maneuver would improve reabsorption of sodium, chloride and water, all of which are 
also excessively lost during space flight. Methods for induction of renal Na and K 
ATPase deserve further investigation in an attempt to prevent osteoporosis and fluid and 
electrolyte loss associated with space travel. The magnitude of renal-related calcium 
losses, coupled with the fact that Skylab astronauts had adequate calcium intake and 
excessive fecal calcium losses, make it wise to consider ways in which renal losses of 
calcium can be curtailed. 
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TABLE 23- V . — Skylab Summary, Plasma Biochemical Results (9 Crewmen) 







(Mean ± Standard error) 




No. 


Sodium 9 

Potassium 

Chloride 

Creatinine 

Glucose 

Osmolality 

Calcium 

Phosphate 

36 

Preflight 

■moq /liter 

141 ± 0.7 

meq/liter 
4.12 ± 0.04 

meq/liter 
97.7 ± 0.5 

mg pet. 
1.26 ± 0.03 

mg pet. 
86.6 ± 0.03 

mO smote 8 
290 ± 0.8 

mg pet. 
9.7 ± 0.05 

mg pet. 
3.4 ± 0.1 

9 

Mission day 
3, 4 

139 ±2 

4.26 ± 0.08 

96.8 ± 0.7 

1.31 ± 0.03 

90.3 ± 2.4 

289 ± 1 

•10.4 ± 0.1 

3.7 ± 0.3 

8 

5, 6 

•137 ± 2 

4.30 ± 0.14 

96.9 ± 0.8 

1.27 ± 0.03 

86.7 ± 1.8 

•287 ± 1 

•10.2 ± 0.1 

•3.6 ± 0.3 

6 

13, 14 

137 ± 1 

4.41 ± 0.15 

•94.7 ± 1.1 

1.28 ± 0.03 

86.7 ± 1.8 

286 ± 2 

•10.2 ± 0.1 

■3.9 ± 0.3 

6 

20, 21 

140 ± 1 

4.25 ± 0.11 

95.7 ± 0.8 

1.35 ± 0.03 

87.0 ± 1.8 

289 ± 2 

■10.1 ±0.2 

•3.4 ± 0.1 

6 

27, 30 

* 138 ± 0.8 

4.25 ± 0.10 

•95.2 ± 0.8 

1.27 ± 0.03 

84.3 ± 2.3 

•287 ± 2 

•10.4 ± 0.1 

•3.9 ± 0.3 

6 

38 

•136 ±2 

•4.05 ± 0.15 

93.5 ± 1.2 

1.31 ± 0.07 

•80.1 ± 2.5 

» 280 ± 4 

10.1 ± 0.2 

•3.1 ± 0.5 

6 

45, 48 

• 137 ± 2 

4.30 ± 0.13 

94.5 ± 0.7 

1.34 ± 0.03 

•84.4 ± 1.4 

287 ± 3 

•10.1 ± 0.1 

•3.8 ± 0.1 

6 

58, 59 

*137 ± 2 

4.19 ± 0.13 

94.0 ± 1.5 

1.38 ± 0.12 

•81.8 ± 2.2 

286 ±4 

•10.1 ± 0.2 

•3.8 ± 0.2 

3 

73 

139 ±2 

3.75 ± 0.20 

94.6 ± 1.2 

1.51 ± 0.05 

80.9 ± 2.2 

284 ± 2 

10.1 ± 0.3 

• 3.9 ± 0.2 

3 

82 

137 ± 0.6 

4.19 ± 0.06 

95.8 ± 0.2 

1.54 ± 0.03 

•81.0 ± 1.2 

•285 ±2 

10.1 ± 0.1 

3.6 ± 0.1 

9 

Recovery (R) 
R + 0 

139 ± 1 

4.18 ± 0.05 

•96.2 ± 1.0 

1.28 ± 0.05 

• 100.5 ± 2.6 

289 ± 1 

• 10.0 ± 0.1 

*3.9 ± 0.2 

9 

R + 1 

139 ± 1 

4.10 ± 0.08 

•96.4 ± 1.0 

1.31 ± 0.06 

92.3 ± 2.8 

289 ± 1 

* 10.1 ± 0.1 

•3.6 ± 0.03 

9 

R + 3.4 

139 ± 1 

4.02 ± 0.13 

96.9 ± 1.0 

1.26 ± 0.06 

•90.5 ± 1.4 

*294 ±2 

9.8 ± 0.1 

3.4 ± 0.2 

6 

R+ 14 

141 ± 0.8 

4.05 ± 0.05 

97.7 ± 1.6 

1.33 ± 0.09 

85.4 ± 0.7 

289 ±2 

•9.4 ± 0.1 

2.8 ± 0.2 

i 

Corrected for Na-EDTA. * P<0.05. 








Table 23- VI. — Skylab Summary, Plasma Biochemical Results (9 Crewmen) 




(Mean 

± Standard error) 




No. 


Cortisol 

Angiotensin I 

Aldosterone 

ACTff 

Insulin 

HGH 

PTH 

39 

Preflight 

ug/100 ml 
122 ± 0.7 

ng/ml per hour 

0.77 ± 0.14 

pg/100 ml 

180 ±25 

Mr* 

nU/mX 

17 ± 0.6 

ng/ml 
1JS ± 0.2 

r f'si 


Mission Day 








9 

3, 4 

12.7 ± 1.6 

1.09 ± 0.24 

176 ± 58 

X 1S2± 4.9 

15 ±2 

•2.1 ± 0& 

17 ±2 

8 

5, 6 

•14.8 ±1.0 

1.75 ± 0.42 

163 ±75 

26.5 ± 9.2 

18 ±6 

1.2 ± 0.3 

16 ±3 

6 

13, 14 

13.4 ± 1.7 

.91 ± 0.28 

252 ± 65 

33.0 ± 8 

18 ±3 

1.5 ± 0.2 

14 ±1 

6 

20, 21 

12.3 ± 1.5 

.52 ± 0.12 

163 ±90 

*11.9 ± 4 

*8 ± 1 

1.2 ± 0.3 

20 ±4 

6 

27, 30 

13.6 ± 2.1 

.45 ± 0.16 

204 ±88 

32.0 ± 7 

20 ±3 

3.2 ± 2.0 

14 ±2 

6 

38 

13.7 ± 1.0 

.72 ± 0.36 

94 ±17 

17.7 ± 11.6 

‘10 ± 1 

1.1 ± 0.3 

15 ±2 

6 

45, 48 

14.3 ± 1.3 

.37 ± 0.10 

118 ±7 

*12.1 ± 5.3 

* 9 ± 2 

1.5 ± 0.5 

18 ±4 

6 

58, 59 

*13.5 ±0.7 

*1.11 ±0.51 

148 ± 31 

32.3 ± 18.7 

*9 ± 2 

1.6 ± 0.4 

18 ±3 

3 

73 

14.5 ± 3.4 

.27 ± 0.08 

117 ± 39 


9 ±4 

0.6 ± 0.1 

24 ±2 

3 

82 

*16.1 ±0.6 

.32 ± 0.04 

142 ± 17 


11 ±4 

0.7 ± 0.1 

25 ± 2 


Rtcovery (R) 







17 ±2 

9 

R + 0 

13.2 ± 2.1 

.71 ± 0.23 

215 ± 74 

23.8 ± 6.3 

20 ±3 

* 2.9 ± 0.6 

9 

R + l 

10.8 ± 1.0 

*2.15 ± 0.55 

*478 ± 77 

* 24.0 ± 7.5 

20 ±2 

*2.8 ± 0.8 

19 ±3 

9 

R + 8, 4 

13.7 ± 3.0 

.86 ± 0.45 

*357 ±65 

* 23.3 ± 2.4 

18 ±2 

*2.6 ± 0.8 

19 ±3 

9 

R + 13, 14 

10.6 ± 0.7 

*.14 ± 0.05 

153 ±35 

38.2 ± 13.9 

17 ±3 

1.2 ± 0.2 

18 ±4 


* P< 0.05. 
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Figure 1. Renal function of nine Sky lab crewmen: Fractional 

excretion of calcium, phosphate, sodium, and chloride. The data 
indicate filtered fractions of these ions that escaped reabsorp- 
tion and were excreted in urine. Peak percent changes shown in 
the figure are relative to preflight values. 
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TIMES PREFLIGHT IN-FLIGHT POSTFLIGHT 


Figure 2. Plasma calcium concentration, GFR, and renal handling 
of calcium in nine Skylah crewmen. Peak percent changes indicated 
in the figure are relative to preflight values. 
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TIME: PREFLIGHT INFLIGHT POSTFLIGHT 


Figure 3. Renal handling of calcium per 100 ml/min GFR in nine 
Sky lab crewmen. Peak percent changes indicated in the figure are 
relative to preflight values. 
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ABSTRACT 


The Matachewan-Hearst Dike swarm consists of plagioclase 
megacryst-bear ing tholeiite dikes that were emplaced over 
an area of approximately 250,000 km2 of the Superior 
Province near the end of the Archean (2.45 Ga). Invaded 
supracrustal rocks vary from metavolcanic dominated 
granitoid-greenstone belts to subprovinces dominated by 
metasedimentary lithologies. The dike compositions vary 
from depleted (iLa/Smln < 0.8) to moderately enriched 
([La/Smin > 1.80) and define linear arrays in compositional 
space (data can be obtained from the author). Neither the 
composition of the immediate country rock nor the 
individual subprovince in which a given dike is located 
exerts significant control on dike composition. Although 
variations in individual element concentrations were pro- 
duced by shallow combined rep 1 eni shment-f ract i ona 1 crystal- 
lization (RFC), the range in incompatible element ratios, 
e.g. ELa/Smin, insensitive to RFC, requires additional 
processes. Relative abundances of the rare earth (REE) and 
high field strength (HFSE) elements suggest either crustal 
contamination or involvement of a slab-derived component. 
Distinguishing between these two alternatives is important 
in determining whether or not crustal formation/evolution 
at the end of the Archean was by subduction-re lated 
processes or some alternative such as under p 1 at ing. 
Combined assimilation-fractional crystallization (AFC) 
models, using crustal rock data from xenoliths and from the 
adjacent Kapuskasing Structural Zone (KSZ) in which lower 
crustal rocks are exposed, suggest that assimilation can 
accommodate the incompatible element variation. The 
assimilant was dominated by silicic granitoids (tonalite to 
granod i or i te ) although mixing of parental dike magmas with 
partial melts of previously underplated basaltic rocks best 
explains some dikes. Assimilation rate, expressed as mass 
of assimilant to mass of cumulate phases, was approximately 
0.5$ total relative mass assimilated varied from 5 to 20 
percent. Although the Matachewan-Hearst magmas possessed 
compositional signatures often attributed to a slab-derived 
component, e.g. high REE/HFSE, Ta/Ta# * 0.17-0.41 and 
falling within the ’destructive plate margin' field on a 
Th-Hf-Ta diagram, modeling indicates that these character- 
istics can be produced by AFC involving observed depleted 
dike compositions and granitoids of the KSZ. Partial melt 
models of a hypothetical amphibo 1 e-bear ing peridotite 
(slab-component enriched C ? 3 ) produce variations in Ti/Ti* 
and (Tb/Yb)n that are not compatible with the MD data. AFC 
processes during underplating of the existing Archean crust 
produced the variations in Incompatible element ratios, 
including the relative depletions of the HFSE. 
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Basaltic magma, derived through the partial melting of 
planetary mantles, is a fundamental product of all the 
terrestrial planets (1). The ascent of basalt magma is one 
of the principal mechanisms of transfering 'juvenile’ 
mantle material to the surface. Subsequent chemical 
differentiation of basaltic magma in that environment is 
important in the generation of new crust and the 
modification of pre-existing crust. Data from missions to 
the other terrestrial planets (i.e. Mercury, Venus, Moon 
and Mars) suggests that all of these bodies have undergone 
differentiation into crust and mantle domains to varying 
extents; Earth and Moon are endmembers representing the 
greatest and least extents of this differentiation. 
Although there is broad agreement that the Earth’s crust 
formed through chemical processing of the mantle and its 
derivatives, the mechanism(s) by which this was 
accomplished is(are) a matter of debate (2, 3). Although 
estimates vary, it is clear that the bulk of the Earth’s 
crust formed prior to the Late Archean (2.5 billion years 
ago (Ga) ) (3). Of interest is whether crustal growth and 
evolution during the Late Archean occurred in a manner 
similar to more modern environments, (e.g. involving 
lateral plate motions and in particular, subduction 
processes), or was early crustal genesis dominated by 
vertical tectonics, involving the process of underplating 
(e.g. 4). In the latter process, large volumes of basalt 
magma pool at the base of the low-density crust. The 
subsequent crystallisation and differentiation of this 
magma adds to mass of the crust. Distinguishing between a 
modern plate tectonic model and the underplating process 
requires knowledge of the chemical nature of the mantle 
source (i.e. is it comparable to modern sources of 
convergent zone magmas) and the extent and nature of 
crustal recycling (i.e. crustal melting and assimilation) 
during the underplating process. In this particular study, 
the major and trace element compositions of dikes from the 
extensive Matachewan-Hearst swarm are evaluated in order to 
characterize their mantle source(s) and to identify the 
process(es) of magmatic differentiation that operated 
during their evolution. 

BACKGROUND 

The Matachewan-Hearst dikes (MHD) represent very large- 
scale basaltic magmatism that occurred in a stable 
continental mass (5) at 2.45 Ga (6). The dikes exhibit a 
north to northwestward trend for a strike distance of over 
700 km and an areal extent of approximately 250,000 km2. 
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The large size of the swarm results in the dikes crossing 
1 i tho-tectonic boundaries along which Archean blocks Ci.e. 
subprovinces) of markedly different lithologic character 
have been juxtaposed. The subprovinces can be described as 
metavo 1 can i c- r i ch belts, e.g. the Abitibi, Wawa and Wabi- 
goon subprovinces, separated by intervening metasedimentary 
belts such as the Quetico and English River Gneiss belts 
(7, 8). A nor theast-trending structural discontinuity, the 
Kapuskasing Structural Zone (KSZ), forms a western boundary 
along the central portion of the swarm and has been inter- 
preted as an oblique cross section through the Archean 
crust (9). Studies of KSZ rocks indicate that the Archean- 
age crust is a complex mixture of mafic to silicic rocks 
that has been metamorphosed to conditions of amphibolite 
to granulite facies (5, 10, 11). Geochemical studies indi- 
cate that the granulitic portion of the Archean crust in 
the KSZ has not been depleted to the extent of more 
'typical' lower crust elsewhere (10, 11). 

Most of the MHD host plagioclase megacrysts whose compo- 
sition (AN85+5) is identical to the plagioclase in Archean 
anorthosites of the Superior Province (12, 13). Trace 
element modeling (14, 15) of a subset of the MHDs suggests 
that a hierarchy of controls on MHD magma composition 
occurred. On a local scale, flow differentiation and 
multiple magma injection produced chemical variability 
observable within single traverses across dikes. Varia- 
tions in the concentrations, observed in both intra- and 
interdike suites, were produced in part by combined 
replenishment-fractional crystallization (RFC). This pro- 
cess resolves the apparent decoupling of major- and trace 
element variations (14, 15) and is consistent with the 
uniform (i.e. non-zoned) plagioclase megacrysts (12). The 
RFC process, however, is incapable of producing the 
observed range in values of incompatible element ratios, 
e.g. (La/Sm)n varies from 0.62 to 2.02. As illustrated in 
figure 1, the processes of simple fractional crystal- 
lization (FC) and RFC, although capable of affecting a 
significant change in the concentration of a given element, 
La in this case, are Incapable of producing detectable 
changes in the (La/Sm)n ratio. A pre-RFC stage of evolu- 
tion, involving other processes, must have operated. Such 
processes, which are not mutually exclusive, include 
combined assimilation-fractional crystallization (AFC) and 
source effects. 

In the AFC process, basaltic magma either melts, 
dissolves or otherwise assimilates crustal material. 
Energy required for the assimilation is derived from the 
latent heat of crystallization released as the magma simul- 
taneously crystallizes. The exact trajectory of composi- 
tional evolution (e.g. fig. 1) will depend on (a) the 
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composition of the parent magma, (b) the composition of 
crystallizing phases, (c) the partitioning of elements 
between the crystallizing phase and the melt, (d) the 
composition of the assimilant and (e) the rate of assimi- 
lation to crystallization. Clearly, the AFC process cannot 
be described in terms of a simple mix between the parent 
magma and the contaminant. 

Source effects refer to the impact that chemical 
heterogeneity in the mantle source and/or variations in the 
melting processes can have on the compositions of the 
derivative melts. Numerous studies have shown that the 
mantle is non-uniform, both on a local and regional scale. 
If different sources are tapped or if a single markedly 
heterogeneous mantle is variably melted (16), variations in 
incompatible element ratios can be produced (fig. 1). 

GEOCHEMICAL MODELS 

The MHD magmas inherited a certain amount of 
compositional variability, either from the mantle source(s) 
from which they were derived, from interaction with lower 
crust through AFC, or from both. In this section, these 
two alternatives are explored in more detail. Algorithms 
of the AFC and melting processes generate hypothetical 
magmatic evolution paths. These calculated trajectories are 
compared to the actual data to develop and test various 
geochemical models. In addition, the calculations help 
constrain models regarding the chemical nature of the lower 
crust and mantle source(s). This Information is pertinent 
to arguments of whether or not plate tectonics operated 
prior to the generation of the MHD magmas. 



Figure 1.- A plot of chondr i te-norma 1 i zed La against 
(La/Sm) for the MHD data (squares). Also shown are 
vectors indicating the direction of compositional 
evolution via processes discussed in text. 
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A compositional characteristic of the MHD rocks 
pertinent to this argument is the relative depletion of the 
high-field-strength elements (HFSEs), e.g. Ta, Ti, Nb, 
relative to the rare earth elements (REEs). Values for 
(La)n/Ta in the MHD range from <40 to >160. Both La and Ta 
are incompatible elements in basaltic systems; Such a range 
would not be expected during FC/RFC or melting processes. 
Fractionation of La from Ta is believed to be charac- 
teristic of magmatic rocks, including more silicic differ- 
entiates, produced by subduction zone processes (17, 18). 

For this study, data for volcanic rocks from a number of 
tectonic environments was compiled and compared with the 
MHD. The fields of oceanic island basalt (OIB), ridge 
basalts (RB) from modern spreading centers, volcanic rocks 
erupted on continents either as subduction zone related 
continental arcs (CA) or as continental flood- and rift 
basalts (CV), and a single example of an oceanic arc, the 
Finger Bay center in the Aleutians (FB) are plotted in 
figure 2. The bulk of the MHD plot within the field of 
CA/CV, distinct from OIB, FB and RB. It would appear that 
the source of the MHD magmas was not similar to that of 
modern oceanic islands. Although the MHD lie within a 
field that includes subduction zone magmas (CA), crustal 
contamination is possible in all plotted cases. Further, 
the MHD data define a linear trend from the RB field 
towards that of crustal rocks (CR). 

Asslml lation-Fractional Crystal 1 ization 

Condie et al. (19) argue that crustal contamination was 
not important in the evolution of the MHD magmas. This 
conclusion is based on two assumptions that, on close 
inspection, do not appear valid. They assume that 1) 
parental magmas of the MHD were identical to modern MORB, 
and 2) the lower crust in the Superior Province was 
depleted in LILE, similar to "uplifted Archean granulite 
terrains". In fact, the more primitive of the MHD magmas 
have some significant differences with modern MORB, e.g. a 
low Zr/Nb ratio (<8.6 in contrast to >20 as assumed by 
these authors). Further, studies of lower crust exposed in 
the KSZ clearly indicate that these granulite facies rocks 
are not depleted in a style typical (?) of other such ter- 
rains (10, 11). Also implicit in the qualitative arguments 
of Condie et al. (19) is that contamination is only a two 
endmember process. As discussed above, it is not. In this 
study, AFC models were developed from crustal rock data, 
e.g. a granodiorite xenolith and tonalitic compositions 
from the KSZ (10). In addition to considering bulk assimi- 
lation of silicic granitoids, models reflecting mixing of 
basaltic magma with partial melts of basaltic precursors 


original 
of pooh 


18-6 



150 


CO 



WIGINAL PAGE | S 
0F P00R QUALITY 


(Ce/Sm)n 

Figure 2.- A discriminant diagram, plotting (Ce/Sm)n 
against (La)n/Ta for rocks from the environments of 
oceanic islands (OIB), spreading centers (RB), an 
oceanic arc (FB) and rocks emplaced on the continents, 
either above subduction zones, as flood basalts or in 
rift zones (CA/CV). The MHD data shown as squares. 
Arrows point in the direction of crustal rocks (CR>. 


garnet in the residuum) are considered. These latter 
assimilants are considered because of the possibility that 
underplating is an episodic process (4), and that the 
partial melting of previously underplated basalts might 
provide a potential contaminant. Experimental results 
suggest that the major element compositions of such melts 
will be tonalitic (20). Remelting of the tonalites 
(partial melts of earlier underplated basaltsC?]) would 
yield liquids in the granodiorite to granite range (21). 

In figure 3, the MHD data are shown as a field on the 
Th-Hf-Ta triangular diagram of Wood (22). The data fall 
within the field of "destruct i ve plate-margin basalts and 
differentiates” (22, p. 12). Menzies et al. (23) point out 
that magmas assimilating Th-rich crustal rocks will be 
displaced towards the Th apex of the triangle. AFC 
calculations in this study substantiate this claim. The 
symbols in figure 3 represent the trajectories of AFC, 
assuming a depleted parental magma (symbols in field N) and 
either silicic or mafic granulites as assimilants. 
Although the various assimilants can not be discriminated, 
it is clear that the AFC process can explain the MHD data. 

The relative abilities of the four assimilants are 
evaluated In figure 4. Incompatible element ratios are 
utilized here to essentially filter out the effect of 
subsequent RFC processes. RFC will not produce significant 
variations in these ratios (fig. 1). Figure 4 indicates 
that the MHD data can best be accommodated by AFC processes 
involving silicic granitoids. The assimilation rate used 
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Hf/3 



Figure 3.- Triangular Th-Hf-Ta discriminant diagram of 
Wood (24). Fields are N-type (N) and E-type (E) MORE, 
alkaline within plate basalts and differentiates (AWP), 
primitive arc tholeiites (PA) and calc-alkaline basalts 
and differentiates (CA). MHD data are shown as shaded 
field and the two points in the N field. Symbols track 
the trajectories of AFC models. 


AFC MODELS (Ma/Mc = 0.5) 



(La)n/Ta 


Figure 4.- Plot of (La)n/Ta versus Ti/Sm for MHD data 
(squares). Vectors represent AFC models (Ma/Mc = 
0.5) for various assimilants as labeled. 
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was 0.5 ( = mass of ass iml 1 ant/mass of cumulates), a rate 
consistent with conditions of the lower crust. Several 
factors contribute to the scatter of the MHD data. The MHD 
data considered in this study occur over an area of 60,000 
km2. Given the lithologic diversity of the lower crust 
(10, ID, and the possibility of variable assimilation 
rates, one can envision a splay of vectors originating from 
the parental composition. Heterogeneities within the 
mantle source are also possible, although isotopic data 
suggest that variations were small and the mantle was 
depleted in character (24, 25). The isotopic support for a 
depleted source is consistent with the choice of the light- 
REE depleted MHD compositions as parental liquids. In 
figure 5 the impact of AFC, incorporating small variations 
in source compositions, has been evaluated. Compositions 
18B and 15B represent a MHD and an older pillowed flow 
rock, respectively; both, therefore, reflect variations in 
mantle chemistry below the Superior Province. The para- 
meters Ta/Ta* and Ti/Ti# (Ta and Ti anomalies) analogous to 
Eu anomalies, are calculated from the extended REE diagrams 
of Briqueu et al. (18). The AFC calculations portrayed in 
this figure indicate that taking into account minor hetero- 
geneities in both the mantle source and the lower crust 
will account for the bulk of the MHD data. 

Further examination of the models and data of figure 5 
indicate that some of the MHD compositions appear to 
require an assimilant derived from a mafic (i.e. basaltic) 
precursor. AFC calculations indicate that this mafic 
source rock must have been light-REE enriched; the use of 
partial melts of depleted compositions such as 15B and 18B 
as contaminants were unsuccessful. Figures 4 and 5 also 
indicate that if the mafic rock had a garnet-free mineral- 
ogy, partial melts of this assemblage would be unsuccessful 
contaminants. The implication of a garnet-bearing residuum 
is that the crust must have been thick enough to stabilize 
this phase. The requirement for a mixed silicic-mafic 
assimilant package emphasizes the heterogeneous nature of 
the lower crust (5, 9 and 11). 

The relative success of the AFC models is illustrated in 
the extended REE diagram given in figure 6. Shown in this 
figure is a potential parental composition (15B), a grano- 
diorite contaminant (3G), and an evolved MHD composition 
(25A). The stippled field is the range of compositions 
produced through AFC (Ma/Mc - 0.5) after the original magma 
mass has been reduced by 10 to 20 percent. It is evident 
that the trace element composition of 25A is consistent 
with such a process. Condie et al. (19) argued that signi- 
ficant assimilation of silicic crust would change the 
silica content of the evolved magma to that uncharacter- 
istic of basalt. Recalling that the AFC process involves 
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Figure 5.- Anomalies of Ta versus Ti (see text and 20) 
for MHD data (squares), various AFC models involving 
possible parental melts 15B and 18B, and partial melt 
trajectory from a hypothetical metasomat ized amphibole- 
bearlng peridotite (AP). Tonalitlc compositions (10) 
from the KSZ shown as outlined area (SG). 


three endmembers (parental magma, assimilant and cumulate 
phases), it is noted that the liquidus assemblage of the 
MHD magmas must have been dominated by c 1 inopyroxene and 
plagioclase (14), both of which have silica contents higher 
than the parental magma. Fractionation of such phases 
offset the higher silica content (in the case of 3G, 
approximately 70 percent) of the assimilant. Qualitative 
calculations indicate that the AFC model portrayed in 
figure 6 produces an evolved liquid having a Si02 content 
consistent with that of 25A. 

Enriched Mantle Source(s) 

In the previous section it was demonstrated that AFC 
processes can produce the compositional characteristics of 
the MHDs. In this section, the possibility that these 
characteristics could alternatively reflect an enriched 
source is evaluated. Of particular interest is whether or 
not a source, previously enriched by subduction zone 
processes, is indicated. The depletion of HFSE relative to 
the LILE and REE, observed in the MHD data, is interpreted 
by some to be unique to subduction related volcanism (e.g. 
18, 27), resulting from the transfer of the LILE and REE 
from the slab as a result of dewatering processes (17). The 
fractionation of LILE and REE from the HFSE, however, may 
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Figure 6.- Extended REE diagram for samples and models 
of the MHD. All compositions normalized to values 
given by Briqueu et al. (20). Sample 15B represents a 
possible parental MHD composition, 3G is a granodiorite 
xenolith considered as potential crustal contaminant, 
and 25 A is intermediate MHD composition. Stippled area 
represents range of compositions produced by AFC (Ma/Mc 
= 0.5} 15B as parent and 3G as contaminant), with F = 
0. 9 to O. 8. 


also be a function of restite mineralogy, in particular Ti- 
rich phases such as sphene or magnetite (18). Further, 
Arculus (28) has argued that the LILE enrichment relative 
to HFSE may not be restricted to subduction zones, but may 
also occur in continental terrains (see figure 2), being 
more a function of process than tectonic environment. In 
agreement, the partial melt calculations discussed above 
indicate that HFSE anomalies can be produced during the 
underplating process by the melting of basaltic precursors. 

In order to evaluate whether or not the MHD magmas may 
have been derived through variable melting of mantle, 
previously enriched above a subduction zone, partial melt 
calculations were performed on a hypothetical amphibole- 
bearlng peridotite. Metasomatic fluids would likely 
produce modal amphibole at the expense of c 1 inopyr oxene 
(29). Progressive melting of such a source (Fig. 5) would 
result in decreases in both Ti and Ta anomalies, with Ti 
being affected to the greatest extent. Few of the MHD data 
conform to such a trend. In figure 7a, the ability of 
amphibole to fractionate middle-REEs (Tb) from heavy-REEs 
(Yb) is noted. The MHD data clearly do not correspond to 
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the trend of variable melting of an amphibo 1 e-bear ing peri- 
dotite. Rather, the MHD compositions are subparallel to 
AFC trajectories involving silicic granitoids as contami- 
nants. Finally, we note that a consistent relation exists 
between (La)n and Ta/Ta* and Ti/Ti* (Fig. 7b). Recall that 
sha 1 1 ow- 1 eve 1 RFC processes affect the (La)n value but 
would have no impact on anomaly values. If variable melt- 
ing of an enriched source was responsible for the anomaly 
values, then, in order to maintain the relationship shown 
in figure 8b, the shallow RFC would have to be related to 
the extent of mantle melting, an improbable relationship. 
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Figure 7.- (a) Plot of (Tb/Yb)n against (La/Yb)n, 
showing partial melt and AFC trajectories. Symbols as 
in figure S. (b) Anomalies of Ta and Ti plotted 
against (La)n for the MHD data. 

CONCLUSIONS 


The simplest model for the MHD magmas is AFC, presumably 
occurring at the base of the crust during underplating. 
Subduction zone enriched mantle sources are not required. 
Trace elements suggest that the mantle sources for the MHD 
were depleted, but possessed a degree of heterogeneity. 
Rates of assimilation were approximately 0.5 ( = Ma/Mc); the 
contaminant mass was <20%. The contaminant was dominated 
by tonal ltes-granodiorltes, similar to xenoliths and rocks 
in the KSZ. Assimilation of partial melts of light-REE and 
garnet-bearing basaltic precursors may have produced some 
the MHD magmas. Apparently, previous underplating-AFC 
processes had already produced a thick crust. The silicic 
granitoid assimilant for the MHD magmas was probably pro- 
duced by earlier processing of underplated mafic crust (4, 
5, 10, 21 and 30). Calculations suggest that the derived 
silicic rocks possess negative Ta and Tl anomalies even 
though they were not the product of subduction. 
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ABSTRACT 


Modern arcs have been subdivided into distinct segments based on features such 
as changes in the alignment of arc volcanoes, variations in seismic parameters (e.g. 
frequency, magnitude and distribution of earthquakes) and/or variations in arc 
lava chemistry. Models for segmentation are varied and include subduction of 
discontinuites on the ocean floor (transforms; ridges) and weaknesses in the 
overriding plate. However, these models are not generally applicable between or 
even within a given arc. Our investigation of 9 ocean-continent convergent 
margins has shown, however, that there are consistent patterns to segment 
development. Segments tend to form in 100-300 km lengths regardless of 
tectonic parameters such as convergence rate or subduction of ocean-floor 
features. Further, the are consistent relationships between the distribution of 
calderas in an arc and segment configuration: calderas are concentrated at segment 
margins and decrease in both size and number towards segment interiors. We 
model arc segmentation in terms of intra-arc stresses caused by buckling of the 
downgoing slab as it is subducted. Deformation and/or tearing of the subducted 
lithosphere induced by buckling provides sites for upwelling of asthenospheric 
material. As the ashenospheric material underplates the overriding slab, crustal 
thinning and extension occur, favoring the emplacement of caldera-forming magma 
chambers. Additionally, the large volumes and short residence times of magmas 
produced at these points precludes extensive modification of their chemistry by 
petrogenetic processes (e.g. fractionation, replenishment, assimilation and/or 
mixing) that are believed to be important in the evolution of more normal "arc- 
type" calcalkaline compositions such as those found at segment interiors. 
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INTRODUCTION 


Volcanic arcs are often broken into discrete segments which are recognized by 
1) changes in the volcanic front, 2) seismological data and/or 3) variations in 
chemistry along the arc (Stoiber and Carr, 1973; Issacks and Barazangi, 1977; 
Hughes et al., 1 980). Although the presence of segments is well documented in 
individual arcs, little work has been done in the systematic comparison of 
segments among a number of arcs to determine what process(es) is (are) 
controlling the formation of segments within an arc. Factors such as discontinuites 
in the subducted slab, subduction of aseismic ridges, and/or crustal weaknesses 
have been invoked by a number of individuals to explain segmentation of given arcs 
(e.g. Issacks and Barazangi, 1977; Carr, 1984; Kay et al., 1982; Hall and Wood, 
1985). To better delineate what factors are involved in segment development, we 
have evaluated segments from a number of arcs, using both a tectonic and volcanic 
perspective. 

INVESTIGATIONS 

Using data from the literature, we compiled a list of pertinent tectonic and 
volcanic paramters for each arc. The majority of the tectonic data was obtained 
from Jarrad's (1986) comparison of stresses at convergent margins. To 
minimize the number of variables influencing segment development, we chose arcs 
which had been active for at least 1 0 ma and were characterized by an oceanic plate 
being subducted beneath a continental plate. Arcs that satisfied these requirements 
were the Aegean, Ryuku, Central America, the Alaska Peninsula, Cascades, Japan, 
Kamchatka, Mainland Alaska, and the Central Andes. We were particularly 
interested in the relationship between stratovolcanoes and calderas within a 
segment as these volcanic features represent sites of focused magmatism and form 
in response to differing tectonic stresses (Hildreth, 1981; Nelson and Wood, 
1987). 

Using the volcanic data base, we noted the size and location of each caldera and 
stratovolcano for each arc. After plotting the volcanic data on maps, we determined 
the volcano/caldera distribution, and located the individual segments within an arc 
based solely on changes in strike of the volcanic front. As mentioned earlier, 
several criteria have been used to identify arc segments. While use of seismology 
and chemical data may be useful in refining segment boundaries in arcs, Issacks 
and Barazangi (1977) noted that the alignment of volcanoes within an arc best 
defines the individual segments. Further, the alignment of volcanoes is preserved 
long after subduction ceases, thus providing information on segmentation on 
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inactive and ancient arcs. This method may also be used where chemical data is 
insufficient to document geochemically anomalous volcanoes and calderas. 

After establishing segments within each arc, we concentrated on isolating 
factors which may control segment location/formation. It became clear that there 
was no consistent relationship between location of segment boundaries and 
subduction of seafloor features. In fact, at the location of segment boundaries, 
there were no ocean floor features currently being subducted which could be called 
upon to cause segmentation. This does not preclude previous subduction of 
ridges, seamounts, transforms etc. which may have initiated earlier 

segmentation of an arc. However, subduction of such features has been called 

upon to alternatively focus (Hall and Wood, 1985), terminate (e.g. Issacks and 
Barazangi, 1977) or have no affect (e.g. Fisher et al., 1981) on volcanism and/or 
segmentation in an arc. Therefore, while subduction of seafloor features may be 
important in facilitating segmentation of some arcs, we do not believe it is a 
prerequisite to segmentation of all arcs. 

Tectonic parameters which relate the interrelationship between both the 
subducting and overriding plate such as convergence rate and strain class were 
initially considered to be important in controlling the segmentation history of an 
arc. Strain class (SC) is a quantitative measure of the amount of extension or 
compression an arc is undergoing, with larger values for SC indicating a more 
strongly compressive arc. Intuitively, it would seem that the larger the magnitude 
of interaction between the two plates (e.g. faster convergence, higher SC values) 
the greater the number of segments that should form. This relationship is clearly 
not supported by our data set. In fact, the poor correlations observed for the 

number of segments vs. strain class (r 2 ~ 0.344; Figure 1) and convergence rate 



Figure 1 .- The variation between number of segments and strain class.factors 
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(r 2 = 0.048; Figure 2) would suggest that plate interactions are not primary 
factors in segmentation. Not only is segment number not controlled by SC, 
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Figure 2.- Variation between segment number and convergence rate. 

segment length is also independent of arc stress. Figure 3 illustrates that there 
is no systematic change in segment length with SC. Although SC 2 is unusual as it 
is represented by only one arc (Ryuku) which has only one defined segment, there 
is no clear relationship between SC and segment length: more extensional arcs do 
not form longer or shorter segments than compressional arcs. Most segments tend 
to form in lengths between 100 and 300 km (figure 4). This is the same range 



Figure 3.- Distribution of segment lengths as a function of strain class. 
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Figure 4.- Frequency of segment lengths. 

that has been noted for individual arcs such as Middle America(Stoiber and Carr, 
1973) and the Eastern Aleutians (Kienle and Swanson, 1983), indicating that 
segment length is relatively constant among arcs. 

Other factors suggested as important in controlling segments are arc age 
and crustal thickness. Again, however, the relationship between these factors and 
segment development is not clear. The poor correlation between number of 

segments and arc age (r 2 »0.197; Figure 5) indicates the segmentation history of 
an arc is not affected by maturity of an arc. In fact, it would seem that 
segmentation is established early in the history of an arc and is maintained 
throughout its development. There is a moderate correlation between crustal 

thickness and segment number (r 2 -0.621; Figure 6) which may indicate that 
thicker crust may be weakened (fractured ?) by subduction and serve to more 
effectively segment an arc. However, careful examination of this relationship 
reveals that arcs with an average crustal thickness of 35-40 km can support 
between 1 and 7 segments. Clearly, crustal thickness is not a sensitive control or 
monitor of segmentation. Our statistical analysis (both correlation and stepwise 
regression) of the available data indicate that the dominant control on the number 
of segments was neither tectonic style nor character of the overriding plate. 
Rather, arc length (hence lateral extent of the subducted slab) is the dominant 

factor in controlling the number of segments (Figure 7; r 2 - 0.769; Correlation 
coefficient - 0.877;). This, coupled with the earlier observation that segment 
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Figure 6.* Variation between segment number and average crustal thickness. 
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lengths cluster at 100-300 km implies segmentation is characteristic of arcs and 
segment number is controlled by the horizontal extent of the downgoing slab. 

As large-scale tectonic features are apparently not important in determining 
the length and number of segments in an arc. we next evaluated what small-scale 
(intrasegment) variations in stresses may be recorded in arcs. Hildreth (1981) 
noted that stratovolcanoes formed when extension was shallow and/or subordinate 
while calderas formed during crustal extension. Previously, we determined that 
caldera diameter reflected the amount of extension experienced in an arc: larger 
calderas formed as extensional stresses increased (Nelson and Wood, 1987). 
Therefore, examining the distribution of stratovolcanoes and calderas within 
segments should yield information on stress variations. To evaluate these 
variations, we divided each segment into percent lengths by effectively "folding" 
the segment in half: the geographic center of each segment was assigned a value of 
50% length, and the percentage length decreased toward segment margins (0%). 
This allowed for easy comparison of segments of varying lengths. Some consistent 
relationships between calderas and percent segment length were observed. At 
segment boundaries, calderas can be located great distances behind the volcanic 
front while at segment interiors, they are located primarily on the front (Figure 
8a). Also, the largest diameter calderas are located at segment boundaries and the 
diameter of calderas decreases toward segment interiors (Figure 8b). Finally, the 
highest percentage of calderas (38%) are located at segment boundaries and the 
percentage progressively decreases toward segment interiors (Figure 8c). These 
consistent relationships, which are seen for all 35 segments, argue for segment 
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boundaries being more extensional than segment interiors. Further, the relatively 
smooth decrease in both size and percentage of calderas from segment boundaries to 
interiors indicates that the change in stress along a segment are gradual and not 
abrupt as often implied. 
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Figure 8.- Distribution of calderas within segments with respect to (a) distance 
behind the volcanic front; (b) caldera diameter and (c) percentage calderas. 
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SEGMENTATION MODEL 


The similarity of segment development in length, number and stresses among 
the various arcs investigated is compelling evidence for a common process or 
processes controlling arc segmentation. We model segmentation as occurring in 
response to subduction of a spherical oceanic lithosphere beneath an overriding 
plate . Frank (1968) was first to note that the subduction processes should be 
influenced by the spherical nature of the subducted lithosphere and Strobach 
(1973) noted that there was a relationship between the curvature of the trench 
and the shape of the subducted lithosphere. Bayly (1982) suggested that the 
arcuate nature of volcanic chains in convergent margins could be produced by 
buckling of the lithosphere as it was subducted. Experimental modeling of 
lithospheric buckling reported by Yamaoka et al. (1986; 1987) and Yamaoka 
(1988) on an arc-arc scale has determined that the subducted lithosphere will 
deform through processes of buckling, with the wavelength of buckling 
corresponding to one arc length. Overall, the experimental models predicted the 
behavior of modern arcs fairly well. Interestingly, however, Yamaoka et al. 
(1987) noted that, for two arcs, the Aleutians and Middle America, the data 
achieved a better fit to experimental predictions after being subdivided 
(segmented?) into separate arcs. 

We believe this large-scale (arc-arc) buckling is applicable to, and provides a 
mechanism for, small-scale (segment-segment) features observed within an arc. 
As the lithosphere is subducted, it will buckle, creating sites of inflection and/or 
tearing in the downgoing slab (Figure 9). Once deformed, portions of the slab can 
act independently of one another, descending at different rates and dips into the 
mantle with the overlying volcanic chain reflecting the changing geometry of the 
subducted model. A similar model was suggested by Carr et al. (1982) to account 
for the transverse boundaries observed in Central America. Further, the points of 
inflection/tearing in the lithosphere can provide sites for asthenospheric 
upwelling. Underplating by asthenospheric material at segment margins results in 
crustal thinning and extension, favoring the emplacement of large caldera-forming 
magma chambers at relatively shallow levels. The heat source for the crustal 
thinning/extension would decrease in intensity away from segment margins. 
Therefore, the distribution of calderas and stratovolcanoes within a segment is a 
natural consequence of buckling of the subducted slab. 

It would seem, then, that lithospheric buckling can explain both the segmented 
nature of arcs and the distribution of calderas and stratovolcanoes within segments. 
Additionally, buckling provides a mechanism for explaining the observed chemical 
variations along segments. Several individuals (e.g. Hughes et al., 1980; Kay et 
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Figure 9.- Schematic model illustrating the relationship between lithospheric 
buckling and segmentation. 

al., 1982; Carr et al., 1982; Wood and Moberger, 1983) have noted the calderas 
and volcanoes located at segment boundaries erupt material that tends to be 
chemically atypical for a given arc. Kay et al. (1982) modeled the chemical 
distribution (calc-alkaline vs. tholeiitic) observed for Aleutian Arc magmas as a 
function of segment boundary magmas (large tholeiitic centers) rising more 
quickly through the lithosphere than those magmas (calc-alkaline trend) which 
were located at segment interiors. Their general model of ascent rate affecting 
the chemistry of Aleutian magmas has implications for magma evolution at other 
arcs and can be accommodated within the framework of lithospheric buckling. 

As mentioned earlier, points of tearing or inflection produced by lithospheric 
buckling would create more extensional environments. This is recorded by the 
presence of large calderas (e.g. Fisher et al., 1981; Carr, 1984; this study) and 
in the transverse alignment of volcanic cones (e.g. Carr et al., 1982) at segment 
boundaries. Extensional environments such as segment boundaries, would favor 
the rapid ascent of magma to shallow crustal levels. In contrast, magmas 
produced at the more compressional segment interiors would ascend more slowly. 
The differences in crustal residence times for magmas produced in these two 
environments is an important control on their chemical evolution. 
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Magmas produced in compressional environments are envisioned to pool, 
perhaps several times at different crustal levels, prior to eruption. In these 
chambers, the magmas will have ample opportunity to undergo processes of 
fractionation, replenishment, mixing and assimilation. These combined processes 
will produce the more typical range of arc-type magma chemistry associated with 
convergent margins. At segment margins, where significant amounts of subcrustal 
heating occurs in response to underplating of asthenospheric material, more 
voluminous magmas can be produced. Important in the evolution of these segment 
boundary magmas is their rapid rise through the crust. Although similar 
petrologic processes are operative throughout the segment (e.g. fractionation 
certainly can occur in segment boundary environments), the shorter crustal 
residence time and larger volumes of segment boundary magmas will place 
restrictions on the degree of chemical modification and evolution they can 
experience. As a result, they will tend not to have chemistries typical of arc 
magmas produced at segment interiors.. Rather than being chemically 
"anomalous", however, we believe that these segment boundary magmas are 
reflecting a normal tectonomagmatic processes at convergent margins which is 
produced by lithospheric buckling. 

As a final point, Yamaoka et al. (1987) noted that the age and length of the 
subducted slab controlled the wavelength of lithospheric buckling. Further, 
Yamaoka (1988) noted that lithospheric bending is "strongly constrained by the 
lateral (along arc) continuation of the lithosphere". While we observed no 
correlation between age or length of the subducted slab, and segmentation history, 
our work on segmentation of modern arcs clearly indicates that the number of 
segments is controlled by the lateral extent of the subducted slab. As segment 
length is relatively constant (100-300 km) it would seem that segmentation 
should occur in all arcs, regardless of discontinuities on the subducting slab 
and/or crustal weaknesses. Features such as these may affect the distribution of 
segment lengths about the 100-300 km range and represent a second-order 
control on segmentation. 

CONCLUSIONS 

Comparison of segment development of a number of arcs has shown that 
consistent relationships between segmentation, volcanism and variable stresses 
exists. We have successfully modeled these relationships using the conceptual 
model of lithospheric buckling of Yamaoka et al. (1986; 1987). Lithospheric 
buckling (deformation) provides the needed mechanism to explain segmentation 
phenomenon: offsets in volcanic fronts, distribution of calderas within segments, 
variable segment stresses and the chemical diversity seen between segment 
boundary and segment interior magmas. 
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Abstract 


Taking pictures with a camera that uses a digital recording 
medium instead of film has the advantage of recording and 
transmitting images without the use of a darkroom or a courier. 
However, high-resolution images contain an enormous amount of 
information and strain data-storage systems. Image compression 
will allow multiple images to be stored in the High-Resolution 
Electronic Still Camera. The camera is under development at Johnson 
Space Center. Fidelity of the reproduced image and compression 
speed are of tantamount importance. Lossless compression 
algorithms are fast and faithfully reproduce the image, but their 
compression ratios will be unacceptably low due to noise in the front 
end of the camera. Future efforts will include exploring methods 
that will reduce the noise in the image and increase the compression 
ratio. 
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Background 


A conventional camera uses film to store images. Light enters 
the lens and strikes the film. A digital camera uses a digital 
recording medium such as memory chips to store images. It converts 
the light entering its lens into a voltage corresponding to the 
intensity of the light. The analog voltage is then converted into a 
number and the numbers are stored in the camera’s memory 
cartridge. A digital camera has the advantage of transmitting images 
to a remote location without the necessity of physically sending film. 
Compared to conventional video, a digital camera has superior image 
quality and provides flexibility and confidence in image 
transmission. Because the image is represented by a file of numbers, 
it can be sent at any transmission rate. Further, error detection and 
correction methods used to transmit digital files virtually guarantee 
that the images arrive distortion-free. 

Developing such a camera is the goal of the Electronic Still 
Camera project, which is under the direction of Don Yeates at the 
Johnson Space Center 1 . It will supplement conventional 35mm 
cameras on the Space Station. Current plans of the Space Station call 
for long intervals between shuttle visits. Using a digital camera to 
transmit images to earth will relieve ground support from waiting 60 
to 90 days for high quality pictures. 

The digital camera will be contained in a housing similar to a 
conventional 35mm camera body, and crew members on the Space 
Station will use it in the same way they would use a conventional 
35mm camera. After taking pictures, the crew member will remove 
the camera's memory cartridge on which the pictures are recorded. 

To preview the pictures, the crew member will use the 
Feedback/Downlink Unit. (See figure figure 1.) The same unit will 
transmit images to earth. The ground station will receive the images. 
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Figure 1.- Electronic Still Camera 
Block Diagram 
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and will have the options of displaying them on a monitor, printing 
them on film and archiving them for future use. 

The amount of memory the camera will need to store an image 
is tremendous. In order to approach 35mm film quality, the 
camera’s resolution will be 2048x2048 pixels. Proper color response 
dictates that a pixel have eight bits of information for each of its red, 
green and blue components for a total of 24 bits. Recording a single 
image requires 12 Mbytes of memory. 

The dimension of the camera body restricts its memory 
capacity to one image. Even with the one Mbit chips now 
commercially available, 12 Mbytes of memory requires 108 chips 
and occupies a volume approximately 4 inches by 4 inches by 6 
inches. There is no room inside the camera for additional memory. 

A memory cartridge with a one-picture capacity is 
unsatisfactory. A crew member inside the Space Station would find 
it irritating to change memory cartridges after every camera shot. 

For an extravehicular activity, crew members would find this 
problem extremely troublesome. A solution to the memory capacity 
problem is to use a compression algorithm in order to store multiple 
images on the memory cartridge. 

Compression is the process of representing an image file in 
fewer bits 2 . The original image file is compressed into a smaller file. 
(See figure 2.) Reconstructing the image from the compressed file 
yields a viewable picture. Important properties of compression 
algorithms are fidelity, compression ratio and speed. Fidelity is the 
degree to which the reconstructed image matches the original image. 
The compression ratio is a ratio of original file size to compressed file 
size. Speed is either expressed as a proportion relative to a function 
of the image resolution or in the number of seconds required to 
compress a file of specified dimension 3 * 4 . 
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Figure 2.- The Compression Process 
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Goals 


The compression algorithm for the Electronic Still Camera must 
have high fidelity, a constant compression ratio and high speed. The 
images from the Space Station must be good enough to be used for 
public affairs. A constant compression ratio is necessary because the 
crew member must know in advance how many shots can be stored 
on a single memory cartridge. The speed is important because the 
logic in the camera must compress over four million pixels between 
shots. Waiting more than five seconds between shots would be 
annoying. For this application, fidelity and speed are more important 
than a high compression ratio. The compression algorithm will most 
likely be implemented as a specialized processor in the camera. 

Method 

I am using all three properties to evaluate the merit of 
compression algorithms. Objective and subjective measure are 
necessary to evaluate fidelity. Such objective criteria as root-mean- 
squared error and signal-to-noise ratio are helpful, but they do not 
adequately measure fidelity because two pictures with the same rms 
error may appear to have drastically different visual qualities 5 . A 
useful subjective method is the pair-comparison method, where 
observers are shown two images at a time and are asked to give a 
preference. For measuring speed, I chose to record the number of 
arithmetic operations per pixel, and the number of separate passes 
through the image file. The first gives a machine-independent 
measure that is more precise than O-notation and the second gives a 
good indication of the amount of parallelism in the algorithm. 

The algorithm will be used to compress pictures similar in 
composition to those taken on the Space Shuttle. I chose nine NASA 
photographs as being “typical" and used a Howtek scanner to convert 
them into computer files. The pictures include shots of payloads, the 
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orbiter, people, equipment and experiments. In addition, I scanned a 
piece of photoblack paper and a piece of white paper to serve as 
baseline images. The eleven files serve as test data for the 
algorithms. 


Results 

Compression algorithms are called either "lossless" 
(information-preserving) or "lossy" (entropy-reducing) 6 . Lossless 
algorithms allow exact duplication of the original image from the 
compressed form. Lossy algorithms lose some information and an 
exact duplicate cannot be reconstructed. 

I chose to concentrate on the lossless algorithms for the initial 
study because fidelity is important and because lossless algorithms 
are generally faster than lossy algorithms. However, the results 
were disappointing. Figure three is a graph giving the maximum 
ratios for compression pixel values in the images. For lossless 
algorithms, there is an upper limit on the compression ratio based on 
entropy 7 . For a second attack, I compressed the relative differences 
between adjacent pixels instead of the pixel values. The results were 
better, but not significantly better. (See figure 4.) 

It was troublesome that the image of uniform photoblack did 
not compress very well. When the image's pixels were magnified on 
a monitor, I found a lot of noise. For lossless algorithms, high 
compression ratios occur when long strings of pixels have the same 
value. Noise corrupts any long string into many shorter strings of 
slightly differing values. The Howtek scanner used to convert the 
photographs into computer files uses a charge-coupled device to 
convert the light into a voltage. The camera will use the same 
technology, so noise is a problem that must be addressed. 
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Figure 3.- Compression Ratios for Hues, Best Case 

X-axis: picture 
Y-axis: compression ratio 
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Figure 4.- Compression Ratios for Differences, Best Case 

X-axis: pictures 
Y-axis: compression ratio 
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Future Work 


One approach to increasing the compression ratio is to "clean up" the 
image before compressing it. One way to clean up the image is to try 
calibrating the camera's charge-coupled devices. Another method is 
to use statistical techniques to compare a pixel with its neighbors and 
adjust its value if it meets certain requirements. 

The act of cleaning up an image loses information. Therefore, I 
want to explore some lossy algorithms, especially the transform 
methods, because they reduce noise and compress the image as one 
process. 
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ABSTRACT 


A strategy to evaluate an expert system is formulated. The 
strategy proposed is based on finding an equivalent classifier to an 
expert system and evaluate that classifier with respect to an opt im al 
classifier, a Bayes classifier. 

This paper shows that for the rules considered that an equivalent 
classifier exists. Also, a brief consideration of meta and meta-meta 
rules is included. Also, a taxonomy of exjDert systems is presented and 
an assertion irwrte that an equivalent classi f ier exists for each type of 
expert system in the taxonomy with associated sets of underlying 
assumptions. 
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1. INTRODUCTION 


It is the purpose of this paper to formulate a formal mathematical 
relationship between an expert system and a classifier. If indeed the 
relationship is unique and constructable, then one will be able to 
evaluate an expert system by comparing that system to a Bayesian clas- 
sifier which is by definition the best estimate. 

In order to formulate a general result one must develop a defini- 
tion of an expert system and a taxonomy of special types of expert 
systems. I have selected my glossary from a document prepared for 
Electric Bower Research Institute entitled, Approaches to the 
Verification and Validation of Expert Systems [1], 

The major purpose of this paper is to develop a theoretical 
foundation for evaluating expert systems. 

2. DESORPTION OF EXPERT SYSTEMS 

A number of definitions of expert systems are found in Artificial 
Intelligence literature, tut the most informative, yet concise one 
seems to be the following. First, they perform large, tediously 
complicated and sometimes difficult tasks at expert levels of per- 
formance. Second, they emphasize domain specific problem-solving 
strategies over the more general 'weak methods' of AI. Third, they 
employ rules of self knowledge to dynamically consider their own 
inference process and provide explanations or justifications for 
conclusions reached. 

An expert system usually consists of four main parts, as shewn in 
Figure 1. 

1. Knowledge Base(s) of structured domain facts and their 
relationships and heuristics (problem-solving rules) , 

2. Inference Engine (s) for controlling the application of facts, 
relationships, and heuristics in solving the problem (s) at hand, 
usually under the control of an overall Meta-Controller module, 

3. Problem Data Base of information about the problem being 
solved and the history of the solution process, and 

4. User Interface providing results/status displays, explana- 
tions, and interaction facilities for user inputs. 

Either AI or Conventional software techniques can be used for the 
storage management of the problem data base and the control of the use 
interface, but AI prog ra mming techniques are always used for repre- 
senting the knowledge base(s) and developing the inference engine (s) . 

In the sections that follow, the differences between expert 
systems and conventional software and the types of expert systems, with 
respect to the ease of verifying and validating them, will be dis- 
cussed. 

2.1 Differences Fran Conventional Software 

Ccnputer scientist claim that expert systems differ from conven- 
tional software programs both in the types of problems they solve and 
their internal structure. Like a human expert, the expert system has 
to accomodate information that is incomplete, erroneous, or mislead- 
ing. Yet, a choice of an action or decision among several alternatives 
must be made. They also claim that a conventional software program 
works correctly only when inputs are complete, of the proper syntax. 
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and the problem is anbiguous. These claims are not at all obvious to 
others in other discipline who view AI and AI tools as simply conputer 
language and an approach to helping solve important engineering 
problems. However, expert systems and their rule base does pose sane 
new specific problems in evaluation of software, The problems expert 
systems are designed to solve generally fall into the following 
categories: interpretation, prediction, diagnosis, design, planning, 

monitoring, debugging, repair, instruction and control. Notice, each 
of these categories contain a function of observing, recognizing and 
may include a call for action. 

Expert systems have many structural similarities with conventional 
software in the modules which perform conventional tasks such as input 
problem data processing, data management, and user interface display 
processing. However, the core of expert systems, including the 
knowledge bases and inference engines, can be different. Symbolic 
representation techniques are primarily used to represent knowledge 
versus conventional software's numeric or table— based techniques for 
representing information. These symbolic representations include 
information about relationships between data items and data item 
aggregations and information about heuristics (procedures and rules) 
for problem solving, which are not represented in conventional data 
bases. Expert systems are developed to solve problems which may not 
have easily formalized or algorithmic solutions, so the problem-solving 
apparatus (Inference engines) must use urxxawentional methods, in- 
cluding heuristic-guided search, symbolic inferencing, generation and 
test of solutions, and constraint-based reasoning. 

It is these differences from conventional software which make 
expert systems interesting and applicable, and yet also less easily 
tested, verified and validated. Because their solution methods are 
often potentially unbounded and not prescribed in a strai^itforward 
"recipe" format like conventional solution methods, they are therefore 
harder to "prove" correct. 

2.2 Expert System Types 

There are many ways of dividing expert systems into types, by 
application, by method of reasoning, etc. Examples of these groupings 
are found in Table 2-1 and Figure 2-1. In this discussion, expert 
systems will be divided into types relative to the complexity and 
difficulty of performing V&V on each type. The types are listed in 

Table 2-0. 

TABLE 2-0.- EXPERT SYSTEM TYPES 

NUMBER NAME 

1 Simple, based on Codified Knowledge 

2 Simple with Uncertainty Handling 

3 Simple, based on Elicited Knowledge 

4 Elicited with Uncertainty Handling 

5 Ocnplex 

6 Ocnplex with Uncertainty Handling 

The first type of expert system. Simple, is developed through the 
straightforward encoding of validated and verified decision tables 
and/or procedure trees. Its search space is small and examined with 
exhaustive search techniques or large and factorable and examined with 
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Figure 1-1. Typical Expert System Components 






Table 2-1 


GENERIC CATEGORIES OF KNOWLEDGE 
ENGINEERING APPLICATIONS 


Category 

Problem Addressed 

Interpretation 

Inferring situation descriptions from sensor data 

Prediction 

Inferring likely consequences of given situations 

Diagnosis 

Inferring system malfunctions from observables 

Design 

Configuring objects under constraints 

Planning 

Designing actions 

Monitoring 

Comparing observations to plan vulnerabilities 

Debugging 

Prescribing remedies for malfunctions 

Repair 

Executing a plan to administer a prescribed remedy 

Instruction 

Dianosing, debugging, and repairing student behavior 

Control 

Interpreting, predicting, repairing, and monitoring 
system behaviors 
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Figure 2-1. Expert System Types Based on Reasoning Methods 
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enumerative search techniques, for which proofs exist of the best 
answer. A single line of reasoning is used and reasoning is monotonic, 
approaching a single, best answer. This type rou^ily corresponds to 
types 1,3, and 4 in Figure 2-1. 

The most important feature besides simplicity for this type is 
that the knowledge to be encoded into Idle expert system has already 
been written down in some form and tested for c o rrec t ness. If the 
codified decision tables or procedure trees must be augmented with 
nan-proved knowledge, then the expert system does not fall into this 
type category. So, V&V of this type of expert system involves merely 
proving that it correctly represents and uses the knowledge, not that 
the knowledge itself is correct. 

The second type of expert system. Simple with Uncertainty 
Handling, has the features of simplicity and proven knowledge in com- 
mon with the first, but also includes certainty factors, fuzzy logic, 
probabilities, or some other method of dealing with uncertain infor- 
mation. This type of expert system roughly corresponds to type 2 in 
Figure 2-1. There are two general types of uncertainty, uncertainty 
on the existence or value of knowledge items or their relationships 
and uncertainty on the rules, reflecting the expert's imxertainty of 
the applicability of the rule on the antecedent conditions or the 
appropriateness of the conclusions. Expert systems may incorporate 
either or both forms of uncertainty. There are various methods of 
combining certainty factors which may be used alone or in combination 
in an expert system. The decision of which combination function(s) 
to use is ccnplex and requires analysis and judgment on the part of 
the knowledge engineer and expert and effects the reliability of the 
system. 

Even if there exist pre-proven uncertainty values on knowledge 
items, or rules, the method of combining those factors most-likely is 
not pre-proven and tested. And, usually, the values themselves do 
not pre-exist and must be elicited from experts. An example is the 
Probability Risk Assessment (IRA) prior event probabilities. Even 
though these probabilities pre-exist and have been determined to be 
correct, to the best knowledge of experts, hew the probabilities can 
be combined and what probabilities to place on rules have not been 
predetermined and tested. Therefore, uncertainty handling compli- 
cates the V&V of simple expert systems by adding the need to V&V the 
certainty factors and their combination method(s) . 

The third type of expert system, Simple based on Elicited 
Knowledge, is probably the most common type in existence. The sim- 
plicity of the previous two types is still evident, probably because 
an expert system building tool or environment was used to develop the 
system. However, the knowledge (items and rules) , or sane portion of 
it, does not pre-exist in a tested form and must be elicited from 
usually one expert. Successive prototypes are developed and tested 
to verify that the knowledge required to solve the problem is actually 
being elicited and that it is encoded properly. 

The fourth type of expert system. Elicited with Uncertainty 
Handling, is a combination of the second and third. Its knowledge has 
been elicited from an expert and it includes certainty factors. An 
example of this type of expert system is the Mycin system in the field 
of m wriirai diagnosis developed by Stanford University. V&V of this 
type of expert system would include testing the correctness of the 
knowledge, the correctness of the knowledge implementation, the 
correctness of the uncertainty factors, and the correctness of the 
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uncertainty combination functions (s) . 

The fifth type of expert system, Complex, eliminates all or some 
of the simplicity assumptions of the previous four, as listed in the 
description of the first type of expert system. The search space (s) 
can be infinite or large and unfactorable, subproblems within the 
problem may interact in various ways, constraint-based reasoning may be 
used to limit search and multiple lines of reasoning may be pursued to 
independently produce candidate answers. This type of expert system 
roughly corresponds to types 5-11 in Figure 2-1. Usually the knew-, 
ledge Is elicited f rom more than one expert, so conflicting heuristics 
may arise. Also, visually expert system building tools are not used to 
implement these systems because the tools do not allow the complexities 
to be built in. 

Examples of the fifth type of expert system include autonomous 
vehicle control and battle management applications, such as those in 
the Strategic Computing Program funded by the Defense Advanced Research 
Projects Agency (DARPA, 1986) . Expert systems of this type are just 
now beginning to work in a primitive fora, and V&V of these systems is 
still a research issue. It is very unlikely that this type of expert 
system will be developed for non-military application ary time in the 
near future, so V&V of this type of expert system will not be discus- 
sed. 

The sixth type of expert system. Complex with Uncertainty 
Handling, arHa the complicating factor of uncertainty handling to the 
fifth, making V&V even more impossible and even more of a research 
issue. This type also will not be discussed. 

3. THE MAIN RESULT 

Consider a simple expert system represented by the following rules 
involving two conditions C, and C_ which can occur or fail to occur and 
four actions Ag, A^, Aj, and A^. 


RO 

If CjAC^, 

then Aq, 

Rl 

If ^ , 

then A^ 

R2 

If (^AOj, 

then A^ and 

R3 

If ^ACj, 

then A-j, 


which can also be represented in vector notation 


RO 

If 

(0, 

0), 

then 

(0) 

Rl 

If 

(0, 

1), 

then 

(1) 

R2 

If 

(1, 

0), 

then 

(2) 

R3 

If 

(If 

1), 

then 

(3) 


Note that the set of rules { RO, Rl, R(2) , R(3) } can be thought of and 
modeled by a function R 

R: (r^ r 2 ) { 0,1,2, 3 ) 

where r- and r 2 € { 0,1 ) 
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Suppose further that 


Pq-pt [C^Cj] , 

P x = Pr [^A^] »■ 

p 3 - Pr [^AC^] , and 

P4 = Pr [C^A^] 1 

are the a priori probabilities that the system are in the respective 
states, 

S 0 ” ' 

S 1 = c i aG 2 ' 

5 3 = Cj^a Cj , and 

5 4 = c 1 aC 2 ' 

Also let data be taken to make a decision as to which state the 
system belongs. The data can be vector valued, dependent or indepen- 
dent, known exactly (deterministic) or observable (random) . In order 
to determine if C x occurs, we follow the usual logic of statistical 

classification theory [2]. Here if 

X = (X^ X^ ..., Xp) T and Y * (Y r Y 2 , ..., Y q ) T 

denotes the observations, then the recognition rules are 

KR 1 If X # 1^, then - 0. 

If X ^ 1^, then = 1. 

RR 2 If Y ^ Rj, then Cj = 0. 

If Y $ Rj, then = 1. 

We \ose the notation Cj to denote that we estimate the value of Cj 

and do not know that value exactly. Ihe classification regions Rl and 
R 2 are selected by the expert or opti m a l ly as Bayes regions. 

Note that the set of rules when described as above is equivalent 
to a classical statistical classification problem. We knew the solu- 
tion to problem; that is, the best classifier is a Bayes classi- 
fier. A Bayes classifier is one in which the regions and Rj are 

selected judiciously to assure that the expected costs of misdassi- 
ficaticn are minimized. 

Probabilities that an expert system ocomits errors can now be 
rmpri-Pri when the probability density functions f(x) and f(y) of X and 
Y are known. That is, 

Pr [C x = 0 = 1] = ^ f x (x)dx 

*1 

Pr [C 2 = 0 l C 2 = 1] = J f 1 (y)dy 

*2 
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and the probability of not making an error is given 

Pr [C x = 0 ] = 0] = I f Q (x)dx 

*1 

Pr [ C 2 = 0 ( C 2 = 0] = J f Q (y)dx 

*2 

Note that if X and Y are independent then various combined prob- 
abilities can be computed easily. If X and Y are not independent then 
the joint probability density function is required. T 

One can stack the observations into a new random vector Z = (x,y) 
and establish Z as the data vector. This form allows one to consider 
the effect of the covariance structure of X and Y. 

Suppose one introduce some "beliefs” concerning the rules. Let us 
model them in the following form. 

RB 0 Pr[R0 is true] = qQ 

RB 1 Pr[Rl is true] = q-j^ 

RB 2 Pr[R2 is true] = q 2 

EB 3 Pr[R3 is true] = q 3 

This in many applications can affect the action, that is, the 
actions are modified to 

RO If C^A C 2 , then PtCAq] = 

and Pr[A^] = l-q^ 

R1 If C^A C^, then Pr^] = q 3 

and Pr[A^] = l-q r 

R2 If C^A C 2 , then PrC^] = q 2 

and Pr[A^] = l-q 2 , and 

R3 If then PIAj] = q 3 

and Pr[A^] = l-q 3 

These rule can be modelled as stochastic actions which are applied 
with a specified probability. This formulation introd uce s the concept 
of stoc hastic actions. This in turn generates a requirement for con- 
flict resolution rules for handling uncertainties of this type. To 
exanple consider the following two types of conflict resolution rules. 

CRRl: If q > 1-q, then select action A, if q $ 1-q, select action iJ . 

CRR2: If the ejected cost of Aq is greater than expected cost of A^, 

select A£; if not Aq. 
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CKO: Generate a uniform random nunfcer N on the interval [ 0 , 1 ], if 0 < 

N£q then select action Agj if q<N51, then select A^. 

The first two rules are deterministic rules while CKO is a sto- 
chastic rule. Clearly, these rules can be extended to more than two 
actions per rule. 

What we have done is given a precise formulation of an expert 
system which in turn gives a mathematical model for evaluating the 
expert system as ccn^ared with on optimal Bayesian classifier. 

Let us new consider the problem of evaluation. Let 

= {X; X<X<X) 

1*2 = {Y; YiY 5Y) 

then we see that directly 

Y 



Ficmre 3.1. The (X,Y) plane partitioned into action regions based 
on an. Expert System if 

X £ 91x1 Y ^ ** 2 * then Aq 

X £ and Y ^ Rj, then A^ 

X ^ R^ and Y % Rj, then A 2 

X^andY*^, then *3 

«here ^ and % are the regions defined by the expert. 

However, a Bayesian classifier would give different regions R^and 
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Figure 3.2. The (X,Y) plane partition into action regions based 
on Bayesian classifier. 

Now consider the following meta rule and a meta-meta rule: 

MR1: If after K trials 


Ci^Cj occurs Kq ti m es 
CfAC^ occurs ti m es 
occurs Kj times 
C 1^2 oocurs K3 times 

then replace ( a's are preselected known parameters) , 

P 0 with a Q P 0 + (l-a Q ) VK 
P 1 with a 1 P 1 + (l-a 1 ) Kj/K 
P 2 with ajPj + (l-aj) Kj/K 
P 3 with a^ + (l-eLj) Kj/K 
and carpute new regions and Rj. 


MMR1: Let K(j) * K + j 50. 

The argument appears reasonable that an expert system can be 
modelled mathematically by finding an equivalent classifier and the 
expert system can then be evaluated by evaluating that classifier. 


4. CONCLUDING REMARKS 

The types of expert systems listed in Table 2.2 and discussed in 
Section 2 of this paper may new be defined precisely by defining an 
associated classifier when different amount and/or kind of knowledge 
differs. 

This allows evaluation of expert systems to be done in a direct 
way compatible with the theory of classification. 
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ABSTRACT 


A study was undertaken to see if magnetic forces can be 
used at the Space Station to attract or repel spacecrafts such 
as the Orbital Manuevering Vehicle (OMV) or the Orbiter. A 
large magnet, in the form of a current loop, is assumed to be 
placed at the Space Station and another one on the spacecraft. 
The expression for the force between the two dipoles (loops) is 
obtained. Using a force of 15 Newtons (3.4 pounds) in order to 
move the spacecraft, the number of ampere-turn needed in the 
current loops was calculated at various distances between them. 

The expression for the force of attraction between a 
current loop and a soft magnetic material was also examined and 
the number of amp-turn needed to provide a force of one-tenth 
of a pound at various distances is also calculated. This 
one tenth of a pound force would be used in a life line system 
for the retrieval of an adrift crewman or tool at the Space 
Station. 

The feasibility of using conventional antenna on the 
Station and the incoming vehicle for attraction or repulsion 
was also examined. 
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INTRODUCTION 


In the earth environment, strong magnetic forces are 
needed at short distances for attraction or repulsion, while in 
space, weak forces are needed at large distances due to the 
microgravity environment. The force needed to move an Orbital 
Manuevering Vehicle (OMV), weighing 20,000 pounds, with an 
acceleration of 150 ALg at low earth orbit will be about 3 
pounds, while a force of 3.4 pounds (15 Newtons) will move the 
Orbiter (weighting 225,000 pounds) with 15 >U. g of acceleration. 
These are very small forces. A study was undertaken to examine 
if a force of about 3 pounds at a distance of approximately 100 
feet, the distance of the object to be attracted to the Space 
Station, could be exerted using magnetic loop antennas, 
conventional electromagnets, or superconducting magnets. The 
temperature in the Space Station, due to background radiation, 
is about 100° K which is below the critical temperature of the 
recently discovered high temperature superconductors (HTSC). 
Provided there is sufficient shielding from the sun, the HTSC 
could be used in the Space Station environment and could use 
cryogen (liquid hydrogen and oxygen) already aboard the Space 
Station for cooling. 

Another problem examined was the retrieval of a drifting 
crewmember using a life line technique. With this technique, 
the crewmember shoots a small iron ball or plate with a 
monofilament line attached to it, at an energized electromagnet 
in the Space Station using a spring action gun. The magnet 
attracts the iron material as it gets closer and holds it. The 
crewmember then uses the line to move to the Station. The same 
technique could be used to retrieve a drifting tool. 


MAGNETIC FIELDS AND FORCES 


Magnetic Field Due to A Dipole 

The magnetic flux density due to a magnetic pole of 
strength 'm' at a distance r is given by 

d — Vn 

B k-TT h* 

" hete = i+TT X lo 7 

and has the units T.m/A or N/A 2 or W/A.m, 


( 1 ) 
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where T stands for Tesla or Weber/m 2 and N for Newton. The 
pole strength 'm' has units of ampere-meter (A.m.). The force 
on a magnetic pole of strength 'm' placed in the magnetic field 
is given by 


F - m B ( 2 ) 

The flux density at a distance much larger than the length of a 
dipole along its axis is given by 2 

vilo AM 

B= 4TT -p 3 -<*> 

where M is the magnetic moment of the dipole given by the 
product of its pole strength and length. 


Force Between Two Dipoles 

The potential energy of a dipole placed in the B field of 
another dipole is given by 

— > — > 

U = — M * 8 (4 


The translational force acting on the dipole is. 


F - - V U 


M • VB 


(5) 


In a uniform B field, there is no net translational force 
on a dipole as the gradient of B is zero. 

Using (3) in (5) we get 


p _ Ao. £ MiM a. 
M-7T f 1 * 


The negative sign in (6) indicates force of attraction. 
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Field Due to A Current Loop on Its Axis 



The flux density at a point on the loop axis duj to 
element of length dl carrying current I is given by 


,t _ M.o I <UX f 

dB “ 4 tt 


an 


( 7 ) 


where r is unit vector^along r. The field is normal to the 
vector "r and element dl. 

The dir$£tion of field due to element dl diametrically 
opposite to dl is shown. The components of the field at right 
angles to the axis cancel, while along the axis they add. 


B COS © — 


Mx> 


The net field along the axis due to a loop with N turns is 

X N C os ©• C 

r 9 ' J 

0 

XNft 

> +X »)% (8) 




for X >> a 



9JV? 

X 3 


( 9 ) 
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where M 


I N A 


( 10 ) 


and has a direction given by right hand cork screw rule. 

Comparison of (3) and (9) indicates that a current loop is 
equivalent to a magnetic dipole with its magnetic moment given 
by (10). 


Force Between Two Current Loops 


Let the magnetic moments of two loops be equal, that is, 
Mj - M 2 - I N A. 


If the two loops, each having a cross section area of lm 2 , 
are separated by a distance r along their common axis, then the 
force between the two loops using (6) is given by 


F 


AX .o 
47T 


jp 4 


which can be written in the form 


IN 


N 


Mir f r 4 

6 A o 


( 11 ) 


A force of 15 Newtons (3.4 pounds) will move the OMV 
(weighing 20,000 pounds) with an acceleration of 150 >Ug and 
the Orbiter (weighing 225,000 pounds) with an acceleration of 
about 15 M.q. The values of ampere-turns (IN) needed to pro- 
vide a force of 15 Newtons between the two loops at various 
distances are calculated from (11) and the results are shown in 
Table 1. The values of ampere-turns needed at distances of 
0.1m, 10m, and 30m are 50, 5 x 10 5 , and 4.5 x 10 6 , respective- 
ly, and they will generate a field of 80/UT, 0.8 T, and 
7.2 T, respectively, at the center of the loops. 


TABLE 1.- DISTANCE VERSES AMP-TURN VALUES 

r(m) 0.1 1.0 5.0 10 15 20 25 

IN 50 5xl0 3 1.25x10 s 5x10 s 1.12x10 s 2x10 s 3.1x10 s 

(Amp-Turn) 


30 

4.5x10 s 
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Figure 2.- Current Loop In A Uniform B Field 


The torque acting on a current loop or a dipole placed in 
a uniform B field is given by 


T = M X B (12) 

There will be a torque acting on the current loop due to 
the earth's magnetic field causing the loop to precess about 
the earth's field. The torque acting on the loop (lm ) with 
10* ampere-turn will have a maximum value of 35 N.m. If the 
loop is fixed in the spacecraft, then the spacecraft will 
precess about earth's field. This problem may be overcome by 
firing jets in the spacecraft that impart equal and opposite 
torque. Another solution to the problem of a single dipole in 
precessing in the earth's field is to use two dipoles with 
oppositely directed magnetic moment vectors ( i . e . , quadru- 
poles), such that the torques due to earth's field cancel out. 



The force that a quadrupole, consisting of a pair of 
current loops (with oppositely directed M vectors), separated 
by a distance Y, exerts on another quadrupole at a distance Z 
away is given by 
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F 


(13) 


AU 18 M) Ma.y 3 ' 

M TT 2 6 


for Z >> y. 

where is the magnetic moment of each of the dipoles in the 
quadrupole #1 and M 2 is corresponding value for quadrupole #2. 

This configuration requires 4 x 10 7 ampere-turn in each of 
the four loops to generate a force of 15 Newtons at a 
separation of 30 meters between the two quadrupole and at a 
separation of 1 meter between the two loops of each of the 
quadrupoles. Thus, there is a high price to pay in going from 
dipole to a quadrupole system in terms of size and ampere-turns 
needed. 


Force Between Current Loop and A Piece of Magnetic Material 
at A Distance 

The force on a ferromagnetic material placed in the B 
field is given by 


F = M-VB ail 


where M is the induced magnetic moment of the magnetic material 
given by 


M - cr V 


(15) 


cr is the induced magnetization in the magnetic material 
of Volume V. 

In ferromagnetic material 1 

^ — yU 0 C H + G”) (16) 


also 


A 0 H 


(17) 
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where ylLis the relative permeability of the ferromagnetic 
alloy. * 


Since 


I 


then 

<r B/^ o 


( 18 ) 


The force of attraction between the current loop and the 
ferromagnetic alloy piece can be written using equations (14), 
(15), (18), and (9). 


F 


3/U 

Mtt a 



(19) 


The force of attraction between a magnet and magnetic 
material decreases inversely as the seventh power of the 
distance . 

M - I N A - magnetic moment of the loop and V - A 2 £X, 
where $X is thickness of the iron plate. 

For a current loop (lm 2 ) and a ferromagnetic alloy plate 
(4" x 4» x 1/8") with - 10,000, the number of ampere-turns 
needed. in the loop to produce an attractive force of one-tenth 
of a pound was calculated at various distances. The results 
obtained indicated values of 5.6 x 10 3 , 1.56 x 10 s , and 1.77 x 
10 7 ampere-turns at distances of lm, 5m, and 10m, respectively. 


When the magnet is in contact with the magnetic material, 
the force of attraction is given by 1 


F 


B*A/AyU 0 


Newtons 


where A is the contact area and B is the flux density in the 
slight air-gap between the magnet and magnetic material. If 
B - 1.0 Weber/m 2 , then F - 4 x 10 5 Newtons/m 2 or 58.0 
pounds/in 2 . 


Weight of the Coil in the Current Loop 

For ampere-turn value of 2 x 10 6 which will be needed to 
exert a force of 15 Newtons at 20m between the two current 
loops as well as between a current loop and sheet of ferro- 
magnetic alloy at a distance of about 7m, let I be 10 4 amps and 
N - 200 turns. Assuming J c - 10,000 Amp/cm 2 (projected value) 
in the high temperature superconductor (HTSC): Crossectional 
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area of the wire needed, A - I/J c - 1 cm 2 


10" 4 m 2 


length of the wire in the loop - * 77 " dN 

■TT( lm) ( 200 ) ■ 200TT meters 
Volume of the coil - area x length - 62.83 x 10 ~ 3 m 3 


Weight of the coil - density x volumne ■ 440 Kgms - 970 pounds 
The density value chosen is the average of cladding material 
and the superconductor density (7000 Kgm/m 3 ) . Length of the 
loop « wire diameter x N r~ — — 

- HA x N " 11 * 3 x 10 x 200 ■ 

2.26 meters. n) T r 

The length of the loop could be reduced by increasing the 
number of layers of the coil. To make a solenoid, the 
superconductor needs to be copper stabilized in order to 
conduct the heat in case of a quench (the superconductor 
turning normal). The weight of the magnet will be 3-4 times 
the weight of superconducting coil. 


Radiation Field of A Magnetic Loop Antenna 



Let the normal to the plane of the loop of radius a be the 
z-axis. The radiated H field due to small loop (a < A/6) at a 
point is given by 8 ' 9 




~P~ M ftl SIN 9 

a a r 


( 20 ) 


Using a phase locked loop in the incoming vehicle, the 
current in the second magnetic loop in the vehicle is kept in 
phase with the B field radiated by the first loop on the 
Station. The force on the second loop due to B field of the 
first loop antenna is — > — ^ 

F = 


if e 


90 


then 


•V VB 

F = /Up 7T N| X, ft, IMala ft 
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( 21 ) 



assuming A 
(21) yields 


A 2 - lm 2 , and I 2 N x - I 2 N 2 - I N, the equation 


IN = 




ioF 


( 22 ) 


A force of 15 Newtons at a 30m separation between the two 
antennas will require 44 x 10 7 and 44 x 10 4 amp-turns at 
frequencies of 40 KHz and 40 MHz, respectively. The high value 
of amp- turns needed will require the use of superconducting 
loop. However, the eddy current losses at high frequency in 
superconductor will produce enough heat to render it normal. 

The inductance of the loop coil is given by 10 

L = A,N b log b/ci/) <23) 

where b and a are loop and wire radii, respectively. The 
reactance of the loop coil (CO L) will be too high to achieve 
the required amp-turns with a moderate sized signal generator. 


DISCUSSION 

The recently discovered HTSC with critical temperature 
above 100° K have a demonstrated ability to remain supercon- 
ducting in magnetic field of at least 30 Teslas. A review of 
high-field and high-temperature superconductors appeared in 
references 12-15. At present, the current density of bulk HTSC 
material is about 1000 amp/cm 2 but work is being done to 
increase it substantially by improving the material fabrication 
and processing techniques. The weight of superconducting coil 
in the current loop was calculated with a projected current 
density of 10,000 amp/cm 2 ; however, if the current density of 
the bulk HTSC material could be increased to 100,000 amps/cm 2 , 
the weight of the superconducting coil would reduce to one- 
tenth of the calculated value. The current density of the thin 
film HTSC is about one million amp/cm 2 . This high value makes 
one optimistic about reaching the goal for bulk material 
current density mentioned above. 

In addition to the coil, a HTSC cylinder weighing about 
20 pounds that can trap magnetic flux due to the Variant 
Incomplete Meisner Effect (VIME) can be used as a permanent 
magnet for fields up to 2.5 Teslas 5 . To make a permanent 
magnet cylinder using VIME, the HTSC cylinder is placed in the 
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magnetic field at T > T c . Then the temperature of the 
cylinder is lowered below its and the flux gets trapped in 
it. making it a permanent magnet as long as T < T c . In case 
the HTSC cyclinder becomes normal and looses its magnetism, a 
magnetizing coil can be wound around to remagnetize it. 

The use of conventional magnetic loop antennas with phase 
locked loop to provide a force of 15 Newtons at 30m distance 
was examined. It was concluded that due to inductance of the 
loop coil (0.2 mH for air core), the required value for 
ampere-turns would not be obtained with moderate sized signal 
generators. 

To retrieve an adrift crewman, use of life line technique 
is proposed. The crewman would use a spring action gun that 
releases an iron ball or plate with a monofilament line 
attached to it towards a current loop that has been energized 
or towards the HTSC permanent magnetic cylinder in the Space 
Station. The iron ball or plate will be attracted to the loop 
or the magnet as it gets closer and will be held by it. Using 
the life line, the adrift crewmember will reach the Station. 

The 2.5 Tesla HTSC permanent magnet or loop with 2 x 10* 
amp-turns will be able to attract magnetic alloy plate (4" x 4" 
x 1/8") from a distance of 15-20 feet. The life line technique 
could be used to retrieve a tool that is drifting away. The 
crewmember will shoot a small permanent magnet at the tool with 
a line attached to it. The magnet will attach itself to the 
tool which can be pulled back using the line. 

Jet propulsion will leave a lot of undesirable plumes 
hitting the Space Station when a spacecraft undocks from the 
Station. The ampere-turns needed to exert a repulsive force of 
about 3 pounds between the two current loops, one placed in the 
Space Station and the other in the spacecraft, will be 50, 5 x 
10 , and 1.25 x 10 s at a distance of 0.1m, lm, and 5m, 
respectively. The size of the current loops needed to provide 
the necessary force is not very large. 


CONCLUSIONS 

Current loops using the HTSC can be used to attract OMV or 
Orbiter with a force of 3 pounds at a distance of about 100 
feet, although the size of the loop and the magnet is fairly 
large. With advances in material processing of HTSC, the 
current density could increase in the near future to a value 
that would reduce the size and weight of loop magnet 
substantially. A series of HTSC permanent magnets employing 
VIME could be used in place of loop magnets. Conventional 
magnetic loop antennas can not be used due to high reactance of 
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the loop at radiating frequencies. A life line technique could 
be used to rescue an adrift crewmember or tool in the early 
stages . 


The size of the current loops needed for undocking a 
spacecraft from the Space Station is not very large. For 
undocking, magnetic repulsion would be used initially to move 
the spacecraft a few meters from the Station and then the jet 
propulsion resumed. This overcomes the problem of undesirable 
jet plumes hitting the Station. 

Conventional magnets yielding field of several Teslas, - 
such as a Bitter magnet cannot be used for space environment 
as the amount of power consumed by them is very large. The 
field of a Bitter magnet is given by 


B ^ II 

N 



Teslas 


where W is the power consumed in megawatts and r is the radius 
of solenoid in cms. To produce a field of 11 Teslas with a 
conventional Bitter magnet with a bore size of 4 cms, will 
require consumption of two megawatts of power! 
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ABSTRACT 


Automation and robotics have played important roles in space activities for many 
years, most notably in planetary exploration and Space Shuttle operations. There are 
several major robotic systems being designed and developed for the Space Station. One 
such system is the Extravehicular Activity Retriever (EVAR). The EVAR will be an 
intelligent free-flying robot which will have the ability to locate, track, maneuver to, and 
retrieve objects which have been detached from the Space Station. Most notably, the 
EVAR is being designed so that it will be able to retrieve crew members who have become 
stranded during extravehicular activity. The EVAR is undergoing a series of ground-based 
demonstrations in the precision air bearing floor facility located at the NASA Johnson 
Space Center. These tests serve as a means of providing proof of concept as well as 
demonstrating system level performance of the EVAR hardware and software. 

Space-based robotics systems are presented with many technical problems which can 
be avoided in robots designed for terrestrial settings. Due to the near absense of gravity in 
space, the manipulation of objects is quite different from what is normally experienced on 
earth. In contrast to the structured environments in which factory robots operate, 
space-based robots must deal with environments containing a high degree of uncertainty. 
'Die harsh nature of lighting in space creates visual perception problems as well. The end 
result is that although many robots have been developed which perform adequately on 
earth, much technology remains to be developed to support intelligent space-based robots 
such as the EVAR. 

This study considers one such area of technology, the development of an approach to 
the visual perception of object surface information using laser range data in support of 
robotic grasping. This is a very important problem area in that a robot such as the EVAR 
must be able to formulate a grasping strategy on the basis of its knowledge of the surface 
structure of the object This paper presents a description of the problem domain and 
formulates an algorithm which derives an object surface description adequate to support 
robotic grasping. The algorithm is based upon concepts of differential geometry namely, 
gaussian and mean curvature. 
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INTRODUCTION 


It is the purpose of this paper to formulate an algorithm for the calculation of surface 
shape and structure from laser range data. This algorithm was developed with the intent of 
using it to assist in the closed loop control of the grasping operation performed by the 
Extravehicular Activity Retriever (EVAR). In order to perform the grasping operation a 
robot such as the EVAR must have knowledge of the structure of the surface of the object 
to be grasped. Using such information, the robot can plan the reach of its manipulator to 
the object and can preshape its hand in such a way so as to accommodate the object and 
ensure a stable grasp. The algorithm calculates information about the surface structure and 
is based upon concepts from differential geometry, namely the gaussian and mean 
curvatures of a surface. The approach to be presented is based in part upon previous 
works by other investigators [1], [2]. 

This report is based upon preliminary results from a study initiated this summer at the 
Johnson Space Center. In addition to the presentation of the algorithm for curvature-based 
surface calculation, this report also describes some of the key criteria relating to the 
determination of grasping strategies. It also presents a framework which allows a robot to 
store essential surface characterization and grasp configuration information in a 
computerized library. Under such a framework, a robot can make use of previously 
generated knowledge about objects and hand configurations to simplify the complexity of 
grasp planning. 


EXTRAVEHICULAR ACTIVITY RETRIEVER 

The use of robots in hazardous environments such as space is rapidly expanding. 
Numerous studies and symposia have been conducted on this very important topic [3], 

The EVAR is a voice-supervised, intelligent, free-flying robot which is currently being 
designed, developed, and demonstrated in a ground-based laboratory facility at the Johnson 
Space Center [4]. When operational, the EVAR will be used to retrieve objects such as 
construction materials, equipment, and tools which have been accidentally separated from 
the Space Station. The EVAR will be in a stand-by mode whenever astronauts perform 
extravehicular activity (EVA). If an astronaut becomes detached from the Station, the 
EVAR will be summoned to locate, fly to, grapple and return the astronaut to the Station. 
The EVAR is being designed so that it can accommodate an unconscious or uncooperative 
astronaut. 

The EVAR consists of six integrated subsystems: the Manned Maneuvering Unit 
(MMU), two robotic arms with grippers, a video and tracking system, a 3D laser imaging 
system, an onboard computer system, and a data and control network. The MMU provides 
the EVAR with the ability to fly about in a manner similar to that of an astronaut during 
EVA. The two robotic arms are anthropomorphic in nature. One of the grippers is a 
dexterous hand while the other is a parallel jaw gripper. The EVAR video and tracking 
system employs two black and white cameras, a tracking system, and a monitor to allow 
for backup supervisory control of EVAR operations by personnel at a remote location. One 
camera enables target acquisition and tracking, while the second can be used to assist 
during target grappling. The 3D laser imaging system is an Ode tics 3D laser radar ranger. 
This subsystem, often called a 3D mapper, consists of a scan unit, an electronics unit, and 
a power supply. The 3D mapper provides direct digital range measurement and is capable 
of scanning a 60 degree horizontal by 60 degree vertical area with a 128X128 raster scan. 
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The frame rate of the system is every .835 seconds. The onboard computer system 
consists of several transputers and other processors. The data and control network enables 
the routing of communication and control information among the various subsystems of 
the EVAR. A three phase ground demonstration program for the EVAR has been planned 
and scheduled. The first phase of the program was successfully completed this past year. 
The second milestone demonstration will be conducted this fall at the precision air-bearing 
floor facility at NASA Johnson Space Center. 

ROBOTIC GRASPING 

Robot manipulation is a complex process that consists of several steps: sensing, task 
planning, trajectory planning, grasping, and path following. With respect to grasping, 
there are three principal grasp selection criteria which are of major concern. They are 
described briefly below. 

Safety. The criteria of safety involves the issue of damage control as applied to three 
areas: (i) the target object must be safe in the sense that there should not be any damage to 
the object during the grasp operation; (ii) the robot must avoid damaging itself as it 
performs the operation; and (iii) the task environment should remain damage-free during 
tiie operation. Whereas the EVAR is being developed to retrieve stranded crew members 
and possibly very expensive Space Station construction tools and materials, requirement (i) 
is an overriding concern in the development of grasping strategies for the EVAR. In 
addition, the EVAR itself will be an expensive piece of equipment, therefore elevating the 
importance of requirement (ii). 

Reachability. This criteria centers on the ability of the hand to reach a particular 
target In other words, there must be a sufficiently large workspace to accommodate the 
movement of the hand during the operation. This requires that the hand be able to reach the 
object at some initial hand configuration and, with the object in hand, to follow a 
collision-free path to a final location. In addition, the workspace of the hand's fingers must 
be large enough relative to the size of the object for the fingers to secure the proper grasp. 

Stability. The grasp must be stable in the presence of external forces which could be 
exerted on the grasped object of the hand during the manipulation process. By stable, we 
simply mean that a grasp does not permit the target to move or slip with respect to the hand. 

It is clear that the above grasp criteria are closely related to the characteristics of the 
targets and the environments which they occupy. In fact, they form the bases upon which 
grasp configurations and operation strategies are determined. In the following section, the 
characteristics of targets which are most critical to the grasping operation are presented. 

CHARACTERISTICS OF OBJECTS 

Objects can be characterized by three different types of attributes (geometric, 
physical, and mechanical) [5]. Geometric attributes include the size of an object, the shape 
of the object, and the shape of the contact surface. With regard to grasping, the size of the 
object is most important in determining the accessibility of the object The shape of the 
object as well as that of its contact surface influence grasp stability and hand preshaping. 
Prior to selecting a grasp, a list of candidate grasp configurations can be extracted from a 
library. Candidate grasp configurations can then be evaluated on the basis of the geometric 
attributes of the object The set of candidate grasp configurations could then be narrowed 
by removing those grasp configurations that would lead to unstable grasps. Finally, the 
physical and mechanical characteristics of the object could be taken into consideration. 

Physical attributes include the mass, mass distribution, and inertia of the objects. The 
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center of mass affects the stability of the grasp configurations. The physical attributes of 
the target are difficult, if not impossible, to extract using solely non-contact methods of 
sensing, such as imaging or laser ranging. This being the case, it would be advantageous 
to develop a library containing the physical properties of objects likely to be retrieved by the 
EVAR [6]. By classifying the target, the EVAR could then retrieve its physical attributes 
and make use of these in the formulation of a grasp strategy. 

Target objects can be grouped into rive categories based upon mechanical attributes. 
These categories are rigid, brittle, elastic, flexibile, and slippery objects. Rigid objects 
exhibit a high degree of stiffness which makes them unlikely to deform during grasping. 
Unlike rigid objects, brittle objects have low stiffness. Potentially, objects of this type can 
be very fragile and thus easily damaged. Grasping brittle objects requires accurate 
information concerning the shape and size of contact areas in that they can be easily 
damaged. Elastic objects are characterized on the basis of behavior of their surfaces about 
the point of contact by a manipulator. When a force is exerted on an elastic object, 
deformation of the surface at the point of contact of the manipulator may appear. When the 
force is removed, the deformation disappears, restoring the object to its original form. 
Elastic objects are typically the easiest objects to grasp. Unlike elastic objects, flexible 
objects have low stiffness. As a result, deformation which result from the application of a 
force remain subsequent to the removal of the force. Normally, detailed shape information 
about the target is not as important for grasping elastic objects. Large contact areas are 
essential for grasping flexible objects. Slippery objects are perhaps the most difficult to 
grasp and handle. Whereas the coefficient of friction at the contact area is very low, the 
stability of the grasp can be diminished. 

It is important to note that the mechanical attributes of objects should be considered in 
conjunction with geometric and physical attributes in the determination of grasp strategies. 
Such design ensures that any drawback in one aspect can be compensated by advantages 
offered by the others. 

VISUAL PERCEPTION OF RANGE IMAGING 

In a broad sense, visual perception of range images can be viewed as the process of 
interpreting measurements made using any of a variety of range sensors. Because the 
success of automated devices depends critically on systems with the ability to sense and 
understand the environment, range image perception has received the attention of many 
investigators. Automatic systems have been designed to work in tightly structured 
environments such as assembly lines in factories. Unfortunately, much of the technology 
associated with systems of this type is not directly transferrable to the EVAR in that the 
EVAR must have the ability to cope with uncertainities in its environment 

The EVAR is equipped with a laser ranging device which is capable of determining a 
range image of the environment surrounding the robot A range image is merely a large 
collection of distance measurements from a known reference coordinate system to surface 
points on objects within a scene. If the distance measurements in a range image are listed 
relative to three orthogonal coordinate axes, the range image is said to be in xyz form. If 
the distance measurements indicate range along 3D direction vectors indexed by two 
integers ( i, j ), the range image is said to be in r jj form. Any image in ry form can be 

converted directly to xyz form, but the converse is not true [7]. The term image is used 
because any ry range image can be displayed on a video monitor, and it is identical in form 

to a digitized video image from a television camera. In addition the term visual perception 
is often associated with the extraction of information from a range image for the same 
reason. 
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Range images can be a source of valuable information for an automated machine. 

The EVAR can use range image information to determine the location and orientation of 
objects in its environment This information can be used to plan the EVAR's movements 
and actions. Surface fitting techniques make use of range data to produce a geometric 
description of objects. In this report we are interested in determining the underlying 
surface structure of targets which the EVAR will be responsible for grasping. An approach 
to extracting the surface structure of a target is critical to the grasp operation because 
manipulators grasp surfaces. The following section presents an approach to surface fitting 
which is based upon curvature concepts taken from the field of differential geometry. 

SURFACE FITTING 

The surface fitting problem can be stated as follows: given three m x n arrays of x, y, 
and z coordinates, we would like to determine a surface that approximates the data in the 
least-squares sense and that is smooth. In a mathematical sense, smoothness is defined as 
the twice-differentiability of the surface at all points. Fitting a surface to an L x L window 
of data is equivalent to computing a surface fit for a roughly rectangular grid of data values. 
Intuitively, a roughly rectangular grid is one that has been obtained from a rectangular 
mesh that has been deformed to fit the surface. We may represent the surface in parametric 
form by the three equations 

x = f(s,t), 

y = g(s,t), 

z = h(s,t), 

where s and t are the parameters, and the functions f, g, and h are tensor products of 
splines in tension. The fitting can be viewed as a mapping from 2D space to 3D space. 

The 2D space is characterized by the parameters s and t and the 3D space is the standard 
cartesian space. 

Principal curvatures at a point on a surface indicate how fast the surface is pulling 
away from its tangent plane at that point. Principal curvatures can be computed by 
estimating partial derivative infomation for the surface. The principal curvatures will 
achieve a local maximum in the area surrounding a discontinuity in the surface, for example 
an edge in the object By thresholding the values of the principal curvatures, we can 
declare the edge points in a range image. 

When a plane passing through the normal to the surface at a point P is rotated about 
this normal, the radius of curvature changes and will be a maximum distance rj fen- a 

definite normal section sj and a minimum 13 for another normal section S 2 - The 

reciprocals 1 <j = \fv\ and ~ U*2 316 called the principal curvatures; the directions of the 

tangents to si and S 2 at P are called the principal directions of the surface at P. Gaussian 

curvature at point P is defined as the product of the two principal curvatures. It can be 
proven that the gaussian curvature depends only upon the coefficients of the first 
fundamental form of a surface and their derivatives. The coefficients of the second 
fundamental form share the previously stated property. The first and second fundamental 
coefficients also determine the surface uniquely up to a rigid body transformation. As a 
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result, the gaussian curvature is said to be an intrinsic property of the surface. The 
principal curvatures at a point on a surface can be computed in terms of the parameters s 
and t using the following approach. Let X(s,t) represent the surface 

X(s,t) = [ x(s,t), y(s,t), z(s,t) ] T . 

Here, bold face notation is used to denote a vector quantity. If we let partial differentiation 
be represented by subscripts, the differential element dX is a vector given by the 
relationship 

dX = X s ds + X t dt, 
where 

X s = dX/ds = [ x s , y s , z s ] T , 


X t = dX/dt = [ x t , y t , z t ] T . 


If we take the scalar product of dX with itself we obtain a relationship which is known as 
the First Fundamental Form of the surface, I, which is given below: 

I = dX • dX = ( X s ds + X t dt ) • ( X s ds + X t dt ) 

= Eds 2 + 2Fdsdt + Gdt 2 > 

where the First Fundamental Coefficients (E, F, and G) are given by: 

E = I X s I 2 = x s 2 + y s 2 + z s 2 , 

F = X s • X t = x s x t + y s y t + z s z u 
G = I X t I 2 = x t 2 + y t 2 + z^. 

The unit normal is then given by 
N = (X s xX t )/IX s xX t l. 

The differential of the unit normal is 
dN = N s ds + N t dt. 

If we take the negative of the scalar product of the differentials of the surface and the unit 
normal, we obtain the Second Fundamental Form for the Surface, 

n = -dX • dN = - (X s ds + X t dt) • (N s ds + N t dt) = 

= L ds 2 + 2 M ds dt + N dt 2 * 
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where the Second Fundamental Coefficients (L, M, and N) are given by: 

L = -X s • N s , 

M = -1/2 ( X s . N t + X t . N s ), 

N = - X t • N t . 

The gaussian curvature K at any point on the surface is defined as the product of the two 
principal curvatures T<j and ^ and can be expressed in terms of the First and Second 
Fundamental Coefficients as 


K = 'k 1 l < 2 = (LN-M2)/(EG-F 2 ). 


The mean curvature H is given by the average of the two principal curvatures. It may be 
expressed in terms of the First and Second Fundamental Coefficients as well. 


H = ( LG - 2MF + NE ) / ( 2( EG-F 2 ) ). 
The principal curvatures are given by, 


kj = H - SQRT( H 2 - K ), 


k 2 = H + SQRT(H 2 - K ). 

The above formulas can be used to compute the gaussian and mean curvatures. Let us 
assume that the gaussian and mean curvatures for a particular point referenced by the 

indices ( i, j ) in a range image have been calculated and are stored in arrays K( i, j ) and 

H( i, j ) respectively. The points on a surface can then be classified according to the signs 
of these quantities: positive, negative, or zero. It turns out that there are wily eight possible 
outcomes for surface characterization based upon the signs of the gaussian and mean 
curvatures. We introduce notation for two functions (modified signum function and the 
surface classification function). The modified signum function is a mapping of a scalar 
argument into one of three values ( 1 , 0 , or -1). For positive e, it is defined as: 

sgn e (y) = 1 if y > e 
0 if I y I < e 
-1 if y <-e . 
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We take advantage of the clarity afforded by use of the modified signum function to 
introduce the surface classification function 


T( i, j ) = l + 3(l + sgn ei (H(i,j))) + (l-sgn e 2 (K(i,j)). 


The positive scalars el ande2 are chosen to be close to zero. They enable the mapping to 
zero of a number within the limits of the numerical accuracy of the machine and the 
algorithm. Figure la gives the mapping of the values of the modified sigpum function 
applied to the gaussian and mean curvatures. Figure lb gives a pictorial interpretation of 
die eight possible surface classifications afforded by the function T( i, j ). 

CURVATURE-BASED SURFACE CALCULATION ALGORITHM 

This section presents an algorithm to compute the object description in terms of jump 
boundaries, internal edges, and regions homogeneous in the sign of their gaussian and 
mean curvatures. The algorithm can be summarize as a sequence of the following steps: 

• Divide the range image into overlapping windows. 

• Detect jump boundaries and fit patches to windows of data not containing jump 
boundaries. 

• Compute the principal curvatures and extract edge points. 

• Classify each nonedge point in a patch as one of the eight possible surface classifications 
presented previously. 

• Group all points of the same type in a patch and its neighboring patches into a region. 

The details of each step in the algorithm are as follows. In the first step, the division of the 
input arrays of ( x, y, z ) coordinates into L x L windows should be overlapped to ensure 
that an internal edge is always contained in a patch. 

In the second step, the relative degree of scatter exhibited by the data is reflected in 
the standard deviation of the euclidean distance between adjacent data points within the L x 
L window. In the vicinity of a jump boundary of an object, the standard deviation of these 
euclidean distances will be relatively high. Jump boundaries can be detected by means of a 
threshold test. After the jump boundaries have been detected, smooth patches based upon 
2D B-splines can be fitted to the data in the window not containing a jump discontinuity 
[8]. In the presence of a jump boundary, the window size can be modified to ensure that 
the window excludes the jump boundary. In order to detect fine detail in the scene it is 
important that the size of the objects be much larger than the size of the L x L window. 

In the third step of the algorithm, we attempt to detect all those points that belong to 
or fall on internal edges of the object. As was mentioned previously, in the vicinity of an 
edge, the values for die principal curvatures will be high and in fact achieve a local 
maximum . As a result, we first determine all points exhibiting principal curvatures above 
threshold. Due to the presence of noise in the data and an inappropriate choice of 
thresholds, clusters of edge points may appear in the vicinity of a true edge position. To 
eliminate these clusters of computed edge points, a suppression of nonmaxima is applied at 
every edge point. That is, we only declare an edge to be present at points where the 
curvature is a local maximum. Nonmaxima suppression is applied in a direction 
perpendicular to the edge directions which is the direction associatated with the maximum 
absolute principal curvature. 

The next step is to group object points into homogeneous regions. All points in a 
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Figure la. Classification of surfaces 
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Figure lb. Surface types 
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patch are classified into one of the types indicated in Figure lb. All points of the same type 
are grouped together into a larger region. Ther merging of points does not require explicit 
fitting of a new surface to this homogeneous region since the surfaces are chosen to be 
smooth. Each homogeneous region is assigned a label depicting its type. Internal edges 
may occur within regions. The extent of regions can be delineated by jump boundaries, 
internal edges, and curvature edges, which we define as places where there is a change in 
curvature-based classification. 

The object representation in terms of regions and curvature-based properties , as 
described above, have the following advantages: 

• it is invariant under transformation of independent parameters of the surface. 

• it is invariant to rigid body transformation and is therefore independent of viewpoint 

EVAR GRASPING 

Operational constraints of the EVAR Demonstration Program demand that grasping 
be performed in a real-time, closed-loop manner. Under Phase II of the program, both 
video and range data will be available for the demonstration. In addition, the target will be 
a stationary object It is assumed that the EVAR hand will be in close proximity (within a 
foot) of the target prior to grasping. 

Probably the single most important property of the human hand is its ability to grasp a 
wide range of objects. The incorporation of a hand on the EVAR will enable the EVAR to 
grasp a wide variety of objects as well. Grasping involves the coordination of many 
degrees of freedom. In an automated system, we wish to reduce the number of degrees of 
freedom associated with robotic grasping. One approach to achieve this goal is based upon 
the notion of classification of grasp configurations. By dealing with grasp configurations 
in a packaged format through classification, one can simplify the amount of planning 
associated with the grasping operation [9]. This approach is based upon the abilty of the 
automated system to recognize surface structures based upon features derived from image 
and range measurements. Features could be constructed using the curvature-based surface 
representation calculation algorithm to describe the underlying suface structure of the object 
to be grasped. These features could be matched against those stored in a surface/grasp 
library [10]. This library would contain grasp configuration which would relate to the 
particular surface structure. In addition, the library would contain the physical and 
mechanical attributes of the associated object. Suppose the EVAR is presented with a 
surface structure similar to one previously encountered and successfully grasped. Under 
such circumstances, the EVAR can make use of information previously generated 
concerning grasping to sieze the target This information would consist of hand preshape 
geometry, grasp approach, and other relevant information. 

If die surface cannot be classified by matching its surface features with those 
contained in the surface/grasp library, it will be necessary to formulate a new grasp 
strategy. In that the hand must grasp surfaces, the information generated using the 
algorithm would be useful in the development of a new grasp configuration. If a new 
configuration is developed and successfully used to produce a stable grasp, then it is stored 
along with the underlying target surface structure in the surface/grasp library. Needless to 
say, the surface/grasp library will be dynamic. As new objects are successfully grasped, 
the library will expand. The overall approach is illustrated in Figure 2. 

SUMMARY 

This paper has presented an algorithm for the calculation of curvature-based surface 
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Figure 2. Framework for grasp operation. 
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characteristics for objects using range data. At the present time, this algorithm is being 
implemented in the programming language C. Although the approach to surface 
description based upon curvature concepts taken from differential geometry appears 
promising, it would be presumptuous to leave the reader with the impression that it will 
provide all the information necessary for secure and stable grasping for the EVAR. It is 
essential to integrate this prototype software, when available, in a real-time systems level 
demonstration in order to fully evaluate its merits. Such a demonstration would bring to 
light any deficiencies which could then be addressed by future refinements. 
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ABSTRACT 


The magnitude, variance and directionality of the velocity 
perturbations of the fragments of a satellite can shed valuable 
information regarding the nature and intensity of the fragmentation. Up 
until now, the only method used to calculate the three orthogonal 
components of the velocity change consisted of inverting the process of 
evaluating the changes in the orbital elements of the fragments due to 
velocity perturbing forces. But the traditional method failed in five 
different cases: one, when the parent satellite's orbit was circular; 
two and three, when the true anomaly of the parent was either 0° or 
180°; and four and five, when the argument of latitude of the parent 
was 90° or 270°. This report describes a new method of calculating the 
velocity perturbations which is free from the shortcomings of the 
traditional method and could be used in all occasions, provided the 
fragmentation data and the orbital elements data are consistent with 
one another. 

The method uses the parent satellite's local frame of reference 
at the time of breakup. The three orthogonal components of the velocity 
change are derived from three simultaneous equations provided by the 
changes in specific energy, specific angular momentum and plane angle 
of the fragment. The fragmentation parameters and the orbital elements 
of the fragmented satellites are taken from the satellite fragmentation 
catalog while the orbital elements of the fragments are read from the 
NORAD and PARCS data files. In order to reverse the decay of the 
fragment's orbit, the orbital elements of the fragment are propagated 
backwards to the time of fragmentation. 

The velocity perturbations of the fragments of over 20 major 
satellite fragmentation events were calculated using the new method in 
this study. This method, together with a recent technique of 
determining the masses of the fragments, have now furnished a complete 
set of information about the fragments, including the mass, velocity, 
momentum, effective cross-sectional area and diameter. The magnitudes 
and variances of the velcotity perturbations of the fragments of 
Himawari rocket and Cosmos 839 satellite have reaffirmed their 
classifications as low intensity explosion and collision induced 
breakups respectively. The preferred directionality of the fragments of 
Landsat 3 rocket indicates that the rocket most likely broke up in the 
"clam" model of low intensity explosion. The residual velocity changes 
of the fragments of Sol wind P-78 satellite provide indications of the 
general direction of the incoming interceptor. The differential 
velocity distribution of the fragments of Delta class of breakups 
could be fitted with a beta function, whereas the cumulative 
distribution of the same followed a general power law. There is 
generally a negative correlation between the mass and the velocity 
change of the fragments. A similar correlation exists between the 
effective diameter and the velocity change. 
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INTRODUCTION 


To date, nearly one hundred artificial Earth satellites have 
fragmented in orbit, thus contributing a large proportion of unwanted 
hazardous material in space (Johnson and Nauer, 1987). To determine 
retrospectively the cause of fragmentation from the debris 
characteristics is one of the principal goals of the debris researchers. 
Over the past few years, we have witnessed significant progress in this 
area. From the debris data, Culp and McKnight, 1985, developed a 
satellite fragmentation test to distinguish between explosion and 
collision induced fragmentation events. Based on the distribution of 
plane change angle and radar cross-sections of known fragmentation 
events, Badhwar, et al . , 1988, devised a scheme to distinguish between 
the three classes of satellite breakups: low intensity explosion 
induced, high intensity explosion induced and hypervelocity collision 
induced breakups. Furthermore, a technique now exists to determine the 
masses of the fragments of a satellite from the history of their orbital 
elements and radar cross-sections (Badhwar and Anz-Meador, 1988). 

One of the most important quantities in satellite fragmentation 
events is the velocity change suffered by a fragment during the breakup. 
For instance, the magnitudes of these velocity changes may be used as 
indications as to whether a fragmentation is due to low intensity 
explosion, high intensity explosion or hypervelocity impact, based on 
the study of Bess, 1975. Also, the directionality of the velocity 
perturbations may give clues as to the nature and intensity of the 
explosion as per the analysis of Benz, et al . , 1987. However, up until 
now, no satisfactory method existed to calculate all three components of 
the velocity change. The traditional method consisted of inverting the 
process of evaluating the changes of the orbital elements due to 
velocity perturbing forces from Meirovitch, 1970, and had been used by 
McKnight, 1986; Kling, 1986; and Tan, 1987. But Tan, 1987, showed that 
the calculated radial component of the velocity change became 
exceedingly large when the parent satellite's orbit was nearly circular 
or when the fragmentation took place near the apsidal points. Similarly, 
the calculated cross-range component of the velocity change became 
singularly high when the argument of latitude was 90° or 270°. 

In order to obtain a method which is free from the singularities 
of the traditional method and which could be used in all occasions, 
Reynolds, 1987, initiated a new approach of estimating the velocity 
changes. This approach parallels an earlier study by Wiesel, 1978, and 
was refined by Badhwar, NASA scientist; Tan, summer fellow; and 
Reynolds, Lockheed scientist at NASA- JSC this summer. The method is now 
fully operational and coupled with the technique of Badhwar and 
Anz-Meador, 1988, can furnish a complete set of information about the 
fragments, including the mass, velocity, momentum, effective cross- 
sectional area and diameter. A brief description of the method is 
contained in this report. The report also summarizes the results 
obtained by analyzing the velocity perturbations in the major breakup 
events calculated by this method. 
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METHOD 


In order to calculate the velocity perturbations imparted to the 
fragments, it is convenient to use the parent satellite's local frame 
of reference at the point of breakup (Wiesel , 1978; Reynolds, 1987). 

The three orthogonal directions are then defined by the radial and the 
down-range directions in the plane of the orbit and the cross-range 
direction perpendicular to the plane of the orbit and along the angular 
momentum vector (Fig.l). In this coordinate system, the velocity^ of 
the parent has the components v r , v d and 0. From elementary 
astrodynamics, we have 

V - Eu(f - j)] J . 


and 


v d = p- [ya(l-e 2 )]^ 
v r = ± (v 2 -v/)* , 


(1) 

( 2 ) 

(3) 


where a is the semi major axis and e the eccentricity of the parent's 
prbit, r is the radial distance to the breakup point and y is the 
gravitational parameter. In Eq.(3), the plus sign corresponds to the 
ascending mode of the satellite (true anomaly v < 180°) whereas the 
minus sign corresponds to descending mode (v > 180°). 

Upon fragmentation, the velocity of a fragment v 1 will have the 
components v r +dv r , v d +dv d and dv x , where dv r , dvj and dv x are the 
velocity perturbations received by the fragment during tne breakup. 

In consequence, the fragment would suffer changes in energy, angular 
momentum and orbital inclination. The change in specific energy is 
given by 


E'-E - j(v' 2 -v 2 ) • - »r (jt- y> . 


dv. 


2 ] - - £ (ir - i) , (4) 


or, K( v r + dv r r + (v d +dv d r + 

where a' is the semi major axis of the fragment's orbit. Similarly, if 
e' is the eccentricity of the fragment's orbit, the change in specific 
angular momentum is given by 

h ,2 -h 2 = | r x v ' | 2 - |r x v^ 2 = ya'(l-e' 2 ) - ya(l-e 2 ) , 

or, r 2 [(v d +dv d ) 2 + dv x 2 - v d 2 ] = ya'(l-e' 2 ) - ya(l-e 2 ) . 

Finally, the plane change angle 0 of the perturbed orbit from the 
unperturbed is given by (vide Fig. 2) 


(5) 


tan 0 = 


dv x 

V 3 ^ 


( 6 ) 
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2 




Equations (4), (5) and (6) comprise three simultaneous equations 
from which the three unknowns dv r , dvj and dv x can be solved in 
principle. After some manipulations and elimination, we arrive at the 
expressions for the velocity perturbations 

dv r - ± [y(f- lr) - \ a' (l-e' 2 )] ! - v r , (7) 

dv d = tMa'O-e' 2 )] 5 - v d . (8) 

and dv x * [wa'(1-e' 2 )] 5 . (9) 

In Eq.(7), the plus sign corresponds to the ascending mode of the 
fragment (true anomaly v' < 180°), whereas the minus sign corresponds 
to descending mode (v' > 180°). Elsewhere, the dominance of v,j over dv^ 
ensures the retention of only the positive sign of the square root in 
Eqs.(8) and (9). Equations (7), (8) and (9) can be shown to be 
equivalent to Reynolds' (1987) equations, which are expressed in 
normalized quantities. 

The plane change angle 0 can be expressed as a function of the 
inclinations i and i' of the parent's and the fragment's orbits 
respectively and X, the latitude of the breakup point (cf. Badhwar, et 
al., 1988). Applying both the cosine laws to the spherical triangle 
AA'B of Fig. 2 and the sine law to the triangles ABC and A'BC and 
substituting, we get 

2 2 i 2 2 i 

0 = + cos" 1 cos 1 cos 1,1 + (cos * ~ cos i) (cos-X - cos i') 2 

cos^X 

Here, the plus sign corresponds to i ' > i and the minus sign corresponds 
to i ' < i . 


The true anomaly v' of the fragment at the time of breakup, which 
dictates the sign of v r +dv r in Eq.(7) is determined from the argument 
of latitude u' and the argument of perigee w' of the fragment at the 
time of breakup as follows. 

V' = U' - U)' . (11) 


Applying the sine law to the spherical triangle A'BC of Fig. 2, we get 


u' = sin' 


/Sin X \ 

[ sT?TT r) 


or 180° - S in _1 (|^4r) 


( 12 ) 


where the first solution corresponds to northbound motion of the 
fragment and the second corresponds to southbound motion at the time of 
fragmentation. Since the argument of perigee is perturbed by the 
oblateness of the earth, the argument of perigee of the fragment at the 
time of observation oj' 0 is different from that at the time of 
fragmentation w' . From the rate of precession derived by King-Hele, 
1964, we get 
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(13) 


U) 1 =0)' 


4.98(5 cos 2 i'-l)it 


where At is the time of observation from the time of fragmentation. 
Here, u>' , u' 0 are expressed in degrees and At is expressed in days. It 
was verified from the successive NORAD data that the true anomaly 
predicted this way was generally accurate to within 2 degrees per year 
of the observed value, thus demonstrating the validity of this 
technique. 


The fragmentation parameters and the orbital elements of the _ 
fragmented satellites are taken from the fragmentation catalog by 
Johnson and Nauer, 1987, while the orbital elements of the fragments 
are read from the NORAD data files. The apogee and perigee heights are 
normally converted into the semi major axis and eccentricity. Thus 


a 

e 

r 

. a' 


2r +h +h 
o a p 

2 ’ 

h -h 
- a P 

2r +h +h * 
o a p 

■ r o * h • 

2r +h’ +h' 
o a p 

2 


(14) 

(15) 

(16) 

(17) 


and e' 


h 1 — h 1 
a p 

2r +h' +K* 
o a p 


(18) 


where r ? is the reference radius of the earth, h a and h p are the apogee 
and perigee heights of the parent, and h' a and h' p are those of the 
fragment. 


Atmospheric drag causes decay in satellite orbits, particularly 
at low altitudes, making both the semi major axis and eccentricity 
smaller. This introduces substantial errors in the velocity 
perturbations calculations when the debris data were measured long 
after the fragmentation event, which is quite often the case. In order 
to compensate for this error, the debris parameters are propagated 
backwards to the time of fragmentation by running the orbital decay 
code of Mueller, 1981 in reverse. For this purpose, the area-to-mass 
ratio of the fragments are first determined by the technique of 
Badhwar and Anz-Meador, 1988. The solar flux in that period, which 
determines the atmospheric densities, is read from the NOAA data base. 
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RESULTS 


In this study, over 20 major satellite fragmentation events were 
studied, each of which produced at least 40 cataloged objects. The 
fragmentation parameters and the orbital elements of the parent 
satellite were taken from the satellite fragmentation catalog by 
Johnson and Nauer, 1987 and the orbital elements of the fragments were 
read from the NORAD data files. In addition, the PARCS data for the 
Sol wind P-78 fragments were used. 

Figure 3 shows the frequency distributions of the velocity 
perturbations of the fragments of Himawari rocket body and Cosmos 839 
satellite. The frgmentation of the Himawari rocket belongs to the Delta 
class of fragmentation due to low intensity explosion, whereas Cosmos 
839 was believed to have been an ASAT target which was destroyed 
through collision (cf. Johnson and Nauer, 1987). The magnitudes of the 
velocity changes are much smaller in the latter case as compared with 
the former, which lend support to the classifications based on the 
study of Bess, 1975. Further, the velocity changes in the Cosmos 839 
fragmentation event have a much smaller variance compared with those of 
the Himawari explosion, which are again consistent with the 
classifications (vide Badhwar and Tan, 1988). 

Figure 4 depicts the two- and three-dimensional plots of the 
velocity perturbations of the fragments of Sol wind P-78 satellite which 
was intentionally destroyed through hypervelocity impact (Kling, 1986). 
The velocity perturbations show considerable scatter and their 
magnitudes are significantly larger than those in the all edged Soviet 
ASAT collisions. Conspicous in the figures are three fragments with 
large velocity changes in all three directions. 

In order to study the directionality of the fragments, we define 
two angles as follows: the colatitude 0 = cos~l (dv r /dv) , 0° < 6 < 180°; 
and the longitude 4> = tan~l (dv x /dv<j) + n 180°, -180° < <j> < 180°, where 
n = 0 for dvd >0, n = 1 for dvj < 0, dv x > 0, and n = -1 for dvj < 0, 
dv x < 0. Note that the colatitude is not the plane change angle defined 
earlier. The angles for each fragment of a breakup event are calculated 
and then plotted. Fig. 5a is an angular plot of the fragments of Landsat 
3 rocket as in a Mercator projection of the world map. The preferred 
direction in which most fragments emerged clearly shows that the 
distribution of fragments was highly anisotropic in this case. Based on 
the analysis of Benz, et al . , 1987, we infer that Landsat 3 was most 
likely broken in the "clam" model of low intensity explosion. Fig. 5b is 
a similar plot with the velocity changes in addition. Since the 
Mercator plots are area-distorted with the polar fragments appearing 
less dense compared with the equatorial, we have calculated the 
fragment flux per unit solid angle, where the elementary solid angle is 
dft = sin 0 d6 d<J>. Fig. 5c exhibits the fragment flux per unit solid 
angle, which is now free from area distortion. 
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Fig. 4. Two- and three-dimensional plots of the velocity perturbations of the fragments of Solwind P-78 
satellite. 








Analysis of the Sol wind P-78 fragmentation event with the PARCS 
data taken immediately after the fragmentation shows that the mean 
velocity changes of all fragments, i.e., the mean dv r , dv<j and dv x do 
not add up to zero but were -17.6, -4.9 and 22.5 meter per second 
respectively. For the lack of orbital element history, the masses of 
the fragments could not be determined. However, assuming equal masses 
for all fragments, we find a net transfer of momentum in the direction 
of 6 = 103°. 1 and <j> = 102°. 3, which could imply that the interceptor 
most probably arrived from the general direction of 6 = 76°. 9 and <J> = 
- 77 °. 7 in the local frame of reference of Sol wind P-78. 

Figure 6a shows the differential distribution of the fragments of 
all Delta rocket breakups (NOAA, Landsat and Himawari) versus velocity 
change. The curve, which shows considerable extension in velocity, can 
be fitted with a beta function of the form 

d m = A (dv/d W“'' (' - ^W 6 ' 1 ■ <' 9 > 

Fig. 6b is a plot of the cumulative number of fragments of the above 
category of breakups with velocity change greater than dv versus the 
velocity change. The cumulative distribution gives a smoother curve 
which is given by a power law of the form 

N(>dv) = A (dv + B)" C . (20) 

The constants a, 8, A, B, C etc. which determine the shapes of the 
distributions, could be used as characteristic parameters of an 
individual or a class of breakup events. 

The correlation, if any, between the masses and the velocity 
changes of the fragments could be found once the masses of the 
fragments are determined by the technique of Badhwar and Anz-Meador, 
1988. Fig. 7 is a plot of mass against the velocity change of the 
fragments of the Spot-1 Viking rocket. Amidst considerable scatter, the 
trend for a negative correlation is clearly noticeable. 

The change in momentum suffered by a fragment during the breakup 
can be computed from the masses and the velocity changes thus 
determined. Fig. 8 is an example of the cumulative distribution of the 
fragments of NOAA 5 rocket versus momentum change. The data points were 
fitted with an expression similar to Eq . ( 20) . 

The effective diameter of the Delta class of fragmentations have 
been computed from the effective cross-sectional area and plotted 
against the velocity change with error bars in Fig. 9. Finally, the 
distribution of the fragments of the same class of events are plotted 
against the velocity change and effective diameter in Fig. 10. A 
pyramid-shaped distribution with unequal slopes on the opposite sides 
is clearly visible. This distribution can be used in the study of the 
evolution of debris cloud after a model low intensity explosion. 
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Fig. 6a. 
Differential 
distribution 
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Fig. 6b. 
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Fig. 8. 
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Abstract 


The quest for the highest resolution microwave imaging and principle of time- 
domain imaging has been the primary motivation for recent developments in time- 
domain techniques. With the present technology fast time varying signals can now 
be measured and recorded both in magnitude and in-phase. It has also enhanced 
our ability to extract relevant details concerning the scattering object. In the past, 
the interface of object geometry or shape for scattered signals has received sub- 
stantial attention in radar technology. Various scattering theories were proposed 
to develop analytical solutions to this problem. Furthermore, the random inversion, 
frequency swept holography, and the synthetic radar imaging, have two things in 
common: (a) the physical optic far- field approximation and (b) the utilization of 
channels as an extra physical dimension, were also advanced. Despite the inherent 
vectorial nature of electromagnetic waves these scalar treatments have brought 
forth some promising results in practice with notable examples in subsurface and 
structure sounding. The development of time-domain techniques are studied 
through the theoretical aspects as well as experimental verification. The use of time- 
domain imaging for space robotic vision applications has been suggested. 
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I. Introduction 


The time-domain imaging is the synthesis of the scattered electromagnetic field 
distribution over an object plane. Over the past few years, an interest has developed 
in the use of time- domain methods or data in order to solve scattering problems. 
This approach has many attractive features, among them the large amount of in- 
formation contained in a single pulse response and simple physical interpretations; 
and improvement in time-domain technology increases the interest in such an ap- 
proach. The purpose of this paper is to present a short review of some techniques 
that have been used in the analysis of time-domain scattering. The problem of in- 
terest is the reconstruction of the shape of a convex scatter from knowledge of high 
frequency far field scattered from the object in response to a known incident field. 

The quest for the highest resolution microwave imaging and the principle of 
time-domain imaging have been the primary motivation for recent developments in 
time-domain techniques. In the last decade, modern techniques in sampling 
devices and the advent of fast pulse generators have brought new technology to 
the practicing radar engineers in the measurements of picoseconds at a greatly 
reduced cost. With the present new technology, fast time varying signals can now 
be measured and recorded both in magnitude and in-phase. 

The goals of NASA are to study the reflection coefficient at and below the sur- 
face of the tiles on the Space Shuttle with special emphasis on time-domain robotic 
vision applications and to continue research into problems related to the explora- 
tion of space and its applications of new technology and techniques. To obtain 
these goals, the Space Station shall be designed as a multi-purpose facility in which 
missions of long duration can be conducted and supported. These missions will 
include science and applications, observation, technology development and 
demonstration, commercial laboratories and production facilities, operational ac- 
tivities such as servicing/maintenance, repair of satellites, support of unmanned 
platforms, assembly of large space systems, and as a transportation node for trans- 
fer to other orbits and planetary missions. Automation and Robotic (A&R) is rapid- 
ly growing and it is an important technical area foreseen to increase productivity 
and enhance astronaut safety. The use of A&R for the Space Station can be viewed 
in two major areas: (1) teleoperated/robotic systems for servicing maintenance, 
repairs, and assembly, and (2) computerized systems to reduce the manpower re- 
quirements of planning, monitoring, diagnosis, control, and fault recovery of sys- 
tems/subsystems. 

In addition to increased productivity through autonomy, the A&R will result in 
increased operational capability and flexibility. Robotic operations for the Space 
Station will involve maintenance/repair of the entire structure including various sub- 
systems, orbiter/satellite servicing, astronaut assistance, equipment transfer, dock- 
ing and berthing, inspection, remote monitoring, rocket staging, telescience, and 
assembly of the station and large structures. To aid the astronaut in various tasks 
and replace him/her for some activities, robots must perform beyond the current 
state-of-the-art by responding to a high degree of environmental uncertainty and 
operational flexibility. In order to accommodate various performance goals in 
robots, design concepts have been proposed by Dr. Kumar Krishen of NASA and 
other researchers. 

Today’s technology enables us to extract pertinent details concerning the scat- 
tering object which conventional radar ranging lacks. The inference of object 
geometry or shape from scattered signals has received substantial attention in radar 
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technology largely due to its academic significant and the understandable commer- 
cial and military values. In the past decade, various scattering theories were 
proposed to develop analytical solutions to this problem. It is well known that ran- 
dom inversion, frequency swept holography, and the synthetic radar imaging have 
two things in common: (a) physical optic far-field approximation, and (b) utiliza- 
tion of the channels as an extra physical dimension. 

Despite the inherent vectorial nature of electromagnetic waves, these scalar 
treatments have brought forth some promising results in practice with notable ex- 
amples in subsurface and structure sounding. Recently, Dr. Kumar Krishen, of 
NASA/JSC, has proposed the use of time-domain imaging for space robotic vision 
applications. A multi-sensor approach to vision has been shown to have several 
advantages over the video approach. 

II. Theory 

The theoretical approach outlines the basic principles of time-domain imaging 
relative to the synthesis of the scattered electromagnetic field distribution over an 
object plane. It was stressed that it is the scattered field, not the total field that is 
of interest, it is known by virtue of the boundary conditions, the tangential electrical 
field must be zero at a perfect conducting surface, where as that of the magnetic 
field is discontinuous by the surface current. The literature review reveals that the 
theory will not change and it is applicable to any arbitrary source waveform. 

The scattering problems from rough surfaces is inherently different in nature 
from that of scattering by other bodies. Usually in the rough surface problem, an 
exact knowledge of the shape of the surface is neither available, nor is it of interest 
to the radar operator. Instead, only average properties of the surface shape enter 
into the problem. It is clearly understood that the last requirement rules out a bound- 
ary-value problem, since the exact boundary is not known. We are interested in 
the relationship between the average scattered field or radar cross section and the 
average surface properties. The radar cross section normalized by the area A, 
defined as, 

7(*<) = 7(*.M) (1) 

6i, 6s, and $s are shown for the bistatic case in figure 1 below. The rough surface 
is assumed to consist of height variations about a mean plane, which is taken as 
the xy-plane. The incident wave line in the sz-plane at a polar angle. 

Z 



Figure 1 - THE "BISTATIC SCATTERING GEOMETRY" 
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Figure 2- THE "SCATTERING GEOMETRY" 


Figure 1 shows the bistatic scattering for a planar rough surface whose mean 
height coincides with the xy-plane and figure 2 shows the scattering geometry. I is 
the plane of incidence and S is the scattering plane. 

The symbol < > refers to an average and the normalized average backscat- 
tering cross section are defined as 

7(*i) = 7(M,*) (2) 

Usually one has a specific purpose or reason for investigating scattering from 
a rough surface, these goals may be divided into three categories: (1) The problem 
of direct scattering; one wishes to know the average properties of the scattered sig- 
nal or cross section when the surface properties of rough surface are known, and 
the scattering information is expressed in terms of the surface properties. (2) The 
problem of inverse scattering; one wishes to obtain statistical information about the 
rough surfaces from a knowledge of the average properties of the scattered field. 
This problem is more difficult in that there appear to be many classes of rough sur- 
faces producing the same average scattering cross section as a function of the bis- 
tatic scattering angle and wavelength. 

The radar cross section (scattered field) can be expressed as 4* times the 
power delivered per unit solid angle in the direction of the receiver divided by the 
power per unit area incident at the target. The factor 4* enters from the definition 
of a solid angle. Assuming for the moment that the propagation path between the 
target and the receiving system is lossless, then this power ratio may be expressed 
as 


4*r a 


(£' • £••) 
(E* • £<•) 




(3) 


Where E s . and H s are the scattered electric and magnetic fields, respectively 
and E* and H‘are the incident fields. The scattered field is defined to be the dif- 
ference between the total field (with the target present) and the incident field. This 
is summarized as 


E s s E t - E*. 


(4) 
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Finding the cross section (a) now becomes a problem in electromagnetic field 
theory. In order that the cross section be independent of r, it is desirable to let r in 
Eq. 3 become arbitrarily large (scalar definition, not a tensor function). 


<r as 4» lim r J 

r-*oo 


(£*•£**) 

(£*£*•) 


. .. W-H”) 

4 ^ r (57T1F)* 


(5) 


To compute the scattered field, one sometimes computes the current induced 
on the target and then treats the target current distribution in terms of an equivalent 
aperture distribution. Antennae are often compared to an isotropic antenna, that 
is, anantenna which radiates uniformly in all directions. 


In the problem of clutter, returns from terrain are not wanted. However, in many 
cases the presence of clutter is unavoidable along with the desired signal. One can 
detect and analyze the desired signal significantly better if more is known about the 
properties of the clutter or noise produced by terrain scattering. It is noted that sur- 
face information is generally known, and the properties of the scattered signal are 
related to Eq. 3. 


The Fourier transform was used in this analysis and it was shown that space- 
time equations can be readily derived from their space-frequency counterparts. 
The equations derived will be applicable to field of arbitrary time variation. 


The source signal is related to the fourier transform, F(w), by: 


(<*) F{y>) = £j{t)J m *dl ( 6 ) m = 


( 6 ) 


The upper case lettering are used for the transformed functions unless other- 
wise noted. The wave number k = w/c, where c is the speed of light in the medium, 
a and p are the direction cosines of the vector r projected onto the xy-plane with 
a =sinecose, p=sinesine. 

It can be assumed that the field distribution over the antenna apertures and ob- 
ject plane is space-time separable; the field strength of each radiating element ob- 
served at some point r is e(r,t) = e(r) f(t - r/c), the field strength, and f(t - r/c) is the 
sine wave traveling at the speed of light. 

Assume the scatter in figure 2 to be placed in a far field, such that the assump- 
tion of uniform impinging planewave over the object can be obtained by the applica- 
tion of the far field Kirchhoff-Huygens principle. 

The basic principle of time-domain imaging over an object plane has been for- 
mulated and discussed by many researchers. It was not until 1983 when Wolfgang- 
M Boerner developed a finite difference-time domain method (FD-TD) with a near 
field to far field transformation using field equivalence, that a method to treat realis- 
tic scattering problems effectively was developed. In this method, the scattering 
problem is analyzed in two steps by treating the relatively complex near field region 
and the relatively simple far field region separately. The method involves first the 
determination of equivalence electric and/or magnetic current tangential to a virtual 
surface surrounding the scatter of interest by using the FD-TD method for a given 
external illumination. The computed near-field equivalent currents are then trans- 
formed to derive the far-field scattering problem and the radar cross section. Since 
the FD-TD method can deal with dielectric, permeable, and in homogeneous 
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materials in a natural manner, it is possible to incorporate most of the physics of 
wave interaction with any complex scatter of interest 

The Kirchhoff-Huygens principle deals with scattering, but the process is not 
as sensitive in analyzing scattering as the method put forth by Wolfgang-M Boerner, 
et al. TheKirchhoff- Huygens method deals with far field scattering, whereas 
Wolfgang- M Boerner treats both far field and near field scattering. 

The theory involving finite difference-time domain (FD-TD) is not without limita- 
tion. This method is based upon the physical optics approximation which is valid 
for scatter whose dimension is significantly larger than the incident wave length and 
is independent of polarization. 

: 1 



i i 


Figure 3 - SCHEMATIC DIAGRAM OF THE TARGET IMAGING SYSTEM 


a fCO *00 
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/(to) = ^-jtvF(w)T(w)E t {ka t ,kp t )exp[-jkr.]. 

Alter 0 


( 7 ) 

( 8 ) 


Where Et (/fat Apt) is the transmitting-antenna anisotropic directional filter and r 
is the distance between the transmitter aperture and an object plane taken as 
reference. 

The second application of the Kirchhoff-Huygens principle and fixed phase 
evaluation can be shown with little calculus that the received waveform, g(t), is 

g(t) = const r F(w)T(w)K(w)[E„(ka 0 , *fl.)40*w)* £,(*»„ kp t )E r (ka r , k0r)E.(ka., kp.) e - ikr ] e~> mt dw. ( 9 ) 

where the direction cosine, the constant and the distance are dependent on the im- 
aging geometry. 

For the monostatic case commonly encountered (see figure 4 below), the 
transmitter and receiver coincide. In particular, if the phase center of the object is 
on boresight which suggests that the target is closely tracked, Eqn. 9 can be con- 
siderably simplified to yield: 

g(t) = const J~ F{w)T{w)R{w)[E a (ka a , kp.) + Ow) J £,(0, 0)£ r (0, 0)£.(*a„ kp.)e- 7 > ir ^ m, dw. ^ ^ 
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The delineation of the objective response is exercised as follows: First, the time 
response in the absence of any scattering objects is recorded. 
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The antenna’s coupling filter is removed when the signal is subtracted from 
other records with the scatter present. Multiple path scattering between the object 
and surrounding obstacles can be gaited out in time when sufficiently large space 
is given. By definition, 

E,(ka„k0,) = r° r° e.(r, y)ezp[2jk(a Q z + (t,y))dxdy. 0 1 ) 

J-ooJ-oo 

The system kernel, K(w), can be obtained once the coupling term has been 
eliminated. 


K(w) = co*<(ii») a / , (w)r(ti;)ie(t«»)£,(0,0)E r (0,0)f: # (0,0)e-^‘ p V 


( 12 ) 


The use of mathematical techniques, provides an opportunity for the backscat- 
tering waveform to be simplified. To give the system kernel, K(w), and to obtain the 
electrical field strength in frequency domain and by the projection theorem, the 
response obtained is equivalent to a slice in the two dimension fourier plane (figure 
4), which is a representation 



Figure 4 - ONE-DIMENSIONAL TARGET FREQUENCY RESPONSE 


of the one-dimensional target frequency response in the two- dimensional fourier 
plane at a particular viewing direction. Equation 9 becomes 


E(ka, k0) s 


C(w)e*>* r » 

*(») 


(13) 


The application of the inverse fourier transform of equation 13 yields the object 
impulse response 

(14) 


r = «(< - 2x' tin -) 

J— CO C 


which can be considered as a series of scans across the object surface where each 
is an integral of the object field which is orthogonal to the viewing direction. The 
spatial resolution is thus achieved by the utilization of the time channel and the min- 
imal resolvable dimension ax is related to the sampling integral T, namely 

cT (15) 

“ 2«m(0) 
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a different impulse response is obtained by the rotation of the object plane about 
the z-axis, which corresponds to another slice in the Fourier plane. Further rota- 
tion over a 180 degrees span will fill up the entire Fourier plane. 

The two-dimension object field distribution is recovered, namely the recnstruc- 
tion strength, e(x, y) 

e(.,y) = - A J R ^E(ka,k0)MyK Vk{os * /,9) ^ (16 ^ 

a smoothing filter A(ka,kp) is introduced to compensate for the amplification of high 
frequency component due to the factor k in field strength equation [Eq. 16]. It is 
not necessary to have A, for noise free data. For practical numerical evaluation, 
the finite sampling rate will impose an upper limit, R, to the wave number k. To 
eliminate abrupt cutoff, which would cause ringing in the resultant image, it is neces- 
sary to introduce equation 16, to yield a smoothed image, equation 17 


A(ka,k0) = Min&£ 

(tt) 

«(*,») = E{ka,k0),in{ji'W ik W ) dk 

and is the basis of our reconstruction field strength. 


(17) 

(18) 


III. Experimental 

The government, private, and public sector have always had a need for time- 
domain imaging. This area of research did not surface until early 1962. The new 
and old researchers in electromagnetic field theory and other scientist provided 
many theoretical ideals on direct and indirect scattering. It was observed that these 
researchers did not have an insight into the circuit design, arrangement of ap- 
paratus and development of the necessary experiments to verify and determine the 
applications of the existing theory. 

Back scattering from metallic objects was experimentally observed and the sig- 
nal processing technique for extracting the object response was performed in 1985 
by Yeung and Evans. It was not until this time when Yeung and Evans put the ex- 
isting theory coupled with experimental verification together and made a projection 
of the long over due applications. After carefully reviewing the literature and the ex- 
perimental verification provided by these two researchers, it was observed that ap- 
plications of time-domain imaging is applicable to many aspects of the private and 
public sector, but the immediate and long range beneficiaries are NASA and the 
military. 

Recently, Dr. Kumar Krishen of NASA/JSC has proposed the use of time- 
domain imaging for space robotic vision application. It has been shown that space 
objects are covered with heat shields, (dielectric materials), and microwave 
penetrate through these shields. A multi-sensor approach to this robotic vision ap- 
plication has been shown to have several advantages over a video approach. 

The authors of this paper will use the necessary apparatus and experimental 
arrangement of Yeung and Evans to demonstrate the viability of impulse imaging. 
Verifying the experimental results should give us confidence in the validity of the 
theoretical method over a broad range of angles and the enormous experimental 
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potential of time-domain techniques in radar- scattering measurements. After 
scrutinizing the theory and experimental results of Yeung and Evans, the authors 
of this paper will look into time-domain robotic vision application and the effects of 
pulse width, polarization, look angles, and phase images of robotic vision coupled 
with, conducting experiments to verify the application of the existing theory and ex- 
periment. 

The target size employed in these experiments are comparable to the source 
pulse width, and it is well known that electrical shape of objects at low frequencies 
could be considerably different from the physical shape. 

Finite source rise-time (hence limited bandwidth) and insufficient source power 
are the major difficulties for achieving the highest resolution and sharp images. The 
system kernel involves the differentiation of the source waveform, and thisadds un- 
certainties in the determination of low-frequency components which is usually 
manifested as a slowly varying baseline in the deconvolved result. 

IV. Applications 

In 1986, a summary of goals and objectives were recommended by the NASA’s 
Advanced Technology Advisory Committee (ATAC) [1]. The robotic applications 
set up by the Committee are: (a) Teleoperation of mobile remote manipulator with 
collision avoidance, (b) Mobile multiple-arm robot with dexterous manipulators to 
inspect and exchange orbital replaceable units, (c) Systems designed to be ser- 
viced, maintained, and repaired by robots. 

The current remote manipulator system (RMS) will be used to manipulate ini- 
tial station assembly in the Space Shuttle. Further versions of the RMS will be 
employed to assemble Space Station and satellite. 

NASA/JSC is currently developing an early robot, the extra- vehicular 
astronaut (EVA) retriever, to perform as an aid to an EVA/man maneuverable unit 
(MMU) astronaut and retrieve the astronaut if circumstances dictate such an opera- 
tion. 


In the development of the space vision systems cost effectiveness, speed, small 
size, lightweight, high reliability and flexibility, and ease of operation must be con- 
sidered. Microwave techniques of time-domain imaging is one of the tools to 
achieve the above purposes. 

Robotic technology has been widely utilized by NASA to explore space; 
however, one of the most challenging problems to deal with is robotic vision. Vision 
is needed in fixed locations to keep track of location of objects, to pick and place 
machines/systems, and to monitor performance of non-automatic systems. In- 
deed, vision provides intelligence and flexibility in locating automated and non- 
automated machines. 

As the robotics era dawns in space, vision will provide the key sensory data 
needed for multi-faceted intelligent operations. In general, the 3D scene/object 
description along with location, orientation, and motion parameters will be needed. 
Sensor complements may include both active and passive microwave and optical 
types with multi-function capability. The fusion of the information from these sen- 
sors to provide accurate parameters for robots provides by far the greatest chal- 
lenge in vision. Furthermore, the compression, storage, and transmission of the 
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information associated with multi-sensor capability require novel algorithms and 
hardware for efficient operation. 

The vision requirements for space robotics are characterized by environmen- 
tal factors and tasks the robot has to perform. The natural space environment con- 
sists of intense light and dark periods. At a normal Space Station altitude, the 
sunlight intensity will fluctuate above 60 minutes of extreme brightness and 30 
minutes of darkness. Furthermore, due to the absence of atmosphere, light is not 
diffused/scattered. The ubiquity of white surfaces intensifies the problem of relying 
on photometric data for object identification/discrimination. A secondary source of 
concern effecting vision is the absence of gravity. For free flying and tethered ob- 
jects may be found due to the lack of disturbance caused by aerodynamic and 
gravitational forces. 

The unpredictable nature of maintenance and repair tasks creates a problem 
in the development of the capability/design of space robots. The vision capabilities 
must be adaptive/versatile to accommodate these uncertainties 

The quest for the highest resolution microwave imaging and the principle of 
time-domain imaging has been the primary motivation for recent developments in 
time-domain techniques. In the last decade, modern techniques in sampling 
devices and advent fast pulse generators has brought new technology to the prac- 
ticing radar engineers in the measurements of picoseconds at a greatly reduced 
cost. With the present technology, fast time varying signals can now be measured 
and recorded both in magnitude and in-phase. 

Conventional radar ranging has been enhanced by modern technology’s ability 
to extract details about the scattering object. 

V. Conclusions 

In this paper various research approaches to time-domain imaging have been 
reviewed with particular emphasis on the mathematical and physical approxima- 
tions made in them. 

Yeung and Evan dealt with the Kirchhoff-Huygens Principles, which deal with 
far field scattering whereas Boerner dealt with both far field and near field scatter- 
ing. We feel that the far and near field scattering technique is a better approxima- 
tion, particularly when applied to a scattering object that is not a perfect conductor. 

Based on the theoretical and experimental verification provided by these 
authors it is necessary that future research be done so as to go beyond the exist- 
ing approximations of Boerner, et al. For instance, a non-trivial generalization, 
should apply to conductivity and non-conducting targets. 
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ABSTRACT 1 

Two objectives were addressed during a 10 week 1988 NASA/ASEE summer 
faculty fellowship at the Johnson Space Center Atmospheric Reentry Materials 
Structures Evaluation Facility (ARMSEF). These objectives were the evaluation of 
mass spectrometry for the measurement of atomic and molecular species in an arc jet 
environment, and the determination of atomic recombination coefficients for RCG 
(reaction cured glass) coated HRSI (high temperature surface insulation) materials 
subjected to simulated reentry conditions. 

Evaluation of mass spectrometry for the measurement of atomic and molecular 
species provided some of the first measurements of point compositions in an arc jet 
tunnel environments. A major objective of this project centered around the sampling 
residence time. A three staged vacuum sampling system pulled the molecules and 
atoms from the arc jet to a quadrupole ionization mass spectrometer in 400 
milliseconds. Conditions investigated included a composition survey across the nozzle 
exit at 3 cm z-distance from the nozzle exit for 3 different currents. Also, a point 
composition survey was taken around a shock created by the the presence of a blunt 
body. 


Results for the composition survey across the nozzle exit showed concentration 
gradients for oxygen and its related secondary species: oxygen atoms, nitric oxide, and 
carbon dioxide. This gradient was more pronounced at higher currents and locations 
closer to the nozzle exit. The gradients were skewed from the centerline with 
maximum concentration occurring at a theta of 1.0 degree. A composition survey 
across the shock showed nitric oxide concentrations dropped while water 
concentrations increased about 6 fold. Other species also showed small changes. By 
using the existing copper pitot tube, little disturbance occurred to the arc jet 
operations. 

Recombination measurements on RCG coated HRSI materials showed low catalytic 
recombination coefficients for both nitrogen and oxygen atoms. The experiments were 
performed on a 5.08 cm radius HRSI hemispherical model (LI-900) coated with a 
reaction cured glass coating. Gas mixtures investigated included pure nitrogen, 11%, 
14% and 23% oxygen in nitrogen. Flowrates and currents ranged from 0.036 to 0.056 
kg/sec, and 300 to 550 amps respectively. The results showed that the nitrogen 
recombination coefficient is relatively constant at 0.014 over a temperature range of 
1350 to 1550 K. These results compare favorably with those reported by Scott (1981) 
and Kolodziej and Stewart (1987). The analyses of air runs were not as conclusive. 
Experimental repeats within a run were good, however, repeats outside of a run 
showed a variation in surface temperatures of 60 K. Oxygen concentration appears to 
have an influence, however, more measurements are recommended. 


1 This report is dedicated to Dr. Frederic A. Wierum - a pioneer in Arc Jet 
Diagnostics 
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INTRODUCTION 


This report is divided into two major sections. The first section will discuss 
mass spectrometry of atomic and molecular species in an arc jet. The second section 
will discuss atomic recombination measurements on RCG coated HRSI materials. 

SECTION 1 

EVALUATION OF MASS SPECTROMETRY FOR ARC JET SPECIES IDENTIFICATION 


BACKGROUND 

Quadrupole ionization mass spectrometry is a well developed analytical technique 
used for the determination of gas compositions via mass separation. The principle 
centers around an ionization filament which emits electrons. These free electrons 
strike a fraction of molecules and atoms as they pass through the ionization cage. 
These molecules and atoms lose an electron as a result of this collision and a charged 
positive ion results. Under the presence of a rf voltage superimposed on a dc voltage 
established by a two sets of oppositely charged rods (quadrupoles), certain charged 
ions can pass through the quadrupoles and strike a detector (a faraday cup in this 
instance). By ramping a potential across the quadrupoles a "scan" can be performed 
which will measure the m/e (mass to charge ratio) present in the mass spectrometer 
at any given time. Ions created in the process can also fragment and thus a parent 
molecule like a nitrogen will have a signal at 14 and 28 plus a small amount at 29 for 
the NN 1 ^ isotope. Further details about mass spectrometry may be found in a book 
by McLafferty. 

Mass spectrometry has been attempted in arc jet environments before with 
limited success. A major limitation was residence time. Real time responses were on 
the order of 30 minutes - a response time much too long. A mass spectrometer 
sampling system developed by R.J. Willey at Northeastern University has enable rapid 
response of an order of 1 second for gas samples from an atmospheric environment. 
This rapid response was achieved by a staged vacuum system. Therefore an earlier 
drawback could be eliminated if the sampling system could be adapted to an arc jet 
environment. 

Another concern in mass spectrometry as applied to arc jet environments is the 
insurance that what is measured is what appears in the stream. The presence of high 
atom concentrations in the stream present a possible sampling problem because of 
atomic recombination in the probe and sample line. Recombination reactions can 
occur in both the gas phase and along the probe walls as the gas sample travels from 
the source to the mass spectrometer. A preliminary study completed at Northeastern 
University showed that homogeneous gas phase recombination reactions could be 
virtually eliminated by the running at pressures below 20 microns at a residence time 
of 1 second. Wall reactions could be eliminated by using a specially designed probe 
constructed an inert wall material (quartz). 



EXPERIMENTAL 


A block diagram of the mass spectrometer system is shown in Figure 1. The 
initial study used the existing pitot tube (constructed of copper) as the sample probe. 

The pitot tube was tied to the mass spectrometer by a 2 m 0.188 I.D. Resistoflex 
tubing (Teflon™ lined) and 1.5 m 1/4" O.D. Eastman Nylon T ^ tubing. Four valves 
separated 3 vacuum stages and the mass spectrometer chamber. These valves were a 
solenoid valve which allowed remote sampling capability, 2 Nupro SS-4H bellow vacuum 
valves and one Nupro SS-4BG-TW bellow vacuum valve. The vacuum pumps were 
standard mechanical vacuum pumps manufactured by Sargent Welch, Duo Seal and 
Cenco. The mass spectrometer chamber was evacuated by a Balzer’s turbo molecular 
pump Model # TPH050 and this pump was backed by an Edward’s Model E2M-1 high 
vacuum pump. 

The mass spectrometer used was an UTI AMX-100. The mass spectrometer has a 
m/e range of 1 to 100. Thus all light gases below a molecular weight of 100 can be 
observed. Signals from the mass spectrometer were interfaced to an IBM-XT computer 
through a Cybord A/D and D/A converter (Isaac Model 91-1). Programming written 
at Northeastern University allowed for the collection of mass scans or selected peak 
scanning as a function of time. 

Experimental conditions investigated were: 

Currents: 400, 800, 1100 amps 

Pitot z-distances from the nozzle: 3 cm, 23 cm, 43 cm 

Radial Positions from the centerline: +/- 7.9, 6.32, 4.74, 3.16, 1.58,0 degrees 

Gases: Air, Nitrogen, and Argon 

RESULTS AND DISCUSSION 


Shown in Figure 2 is the mass spectrometer response at m/e- 30 (Nitric Oxide) 
for a composition survey at 3 current levels across an arc normal to the nozzle exit 
plane and located 3 cm from the nozzle exit. Several phenomena are evident. As the 
current is increased, nitric oxide concentrations increase, in this case from 2% to 11% 
at the centerline. Further, the results showed more nitric oxide at the centerline 
then at the edges of the nozzle thus indicating a concentration gradient across the 
nozzle. Gradients were also observed for oxygen molecules, oxygen atoms, and carbon 
dioxide. 

Also shown in Figure 2 is a pitot pressure survey. Note that it does not show 
as high as % maxima. Thus the peaks seen at 800 and 1100 amps for nitric oxide 
suggest that centerline concentrations for this species are higher. Another source for 
nitric oxide may be the reaction of nitrogen atoms with oxygen atoms in the probe. 
However as is discussed below, nitrogen atoms are reaching the mass spectrometer and 
it is currently suspected that the nitric oxide measured is coming from the arc tunnel. 
The skewness shown in the figure may represent an arc column slightly out of line. 
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Figure 1. Block Diagram for the Mass Spectrometer System. 
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Figure 2. Example of Nitric Oxide Concentration fomComposition Survey across the 
Nozzle Exit Conducted at Three Currents. 
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Figure 3 shows the corrected mass spectral response at m/e » 14 (N), 16 (O), 28 
(N 2 ), 30 (NO), and 32 (O 2 ). Corrected response means that the signals for each of 
the m/e ratios have had parent molecule fragments subtracted. For example, the m/e 
signal at 14 will have a contributions from N 2 and NO. Thus the 14 signal was 
corrected by subtracting the contribution due to N 2 and NO. Note in this figure how 
O 2 trends drop slightly as current is increased, however, it is a significant response. 
The source of the O 2 is suspected to be from wall recombination of oxygen atoms. 
Using a wall recombination coefficient of oxygen of 0.06 for copper lined sections of 
the probe system and 0.0005 for the inert wall probe areas results in an approximation 
that 70% of the oxygen atoms may be lost to wall recombination. It is estimated that 
this figure can be lowered to 30% by simply using the quartz sample probe. 

High nitrogen atom concentrations are somewhat puzzling. The homogeneous 
reaction models predicted that in the presence of O atoms, N atoms would convert to 
nitric oxide in the probe. However, the m/e signal at 14 is significant at about 14% 
of the total mass spectral response even after the signal is adjusted for fragments 
from N 2 . Thus, nitrogen atoms do reach the mass spectrometer. 

An indication of nitric oxide concentrations around a shock is shown in Figure 4. 
Figure 4 shows the nitric oxide concentration decreasing dramatically inside a shock. 
An explanation for this observation is that nitric oxide dissociates across the shock to 
nitrogen and oxygen atoms. If the source of nitric oxide was probe reactions, one 
would expect that nitric oxide concentration would increase across a shock. 
Spectrographic work by Willey (1987) showed the primary emitter in the free stream 
to be nitric oxide. Thus conditions may exist where nitric oxide concentrations are 
above 1% in the arc jet. Mappings for other species showed the definite shock 
boundary, however, their concentrations were relatively constant on both sides of the 
shock with the exception of water. Water increased about 4 to 6 fold on the 
backside of the shock. This result suggests that some moisture is being pulled in 
from the ejector system. Further information about this project will be available in a 
separate report. 

CONCLUSIONS 

1. Sampling times of less than a second can be achieved. 

2. Nitric oxide is a good probe molecule. 

3. Concentration gradients occur across the nozzle and are strongest for nitric oxide. 

RECOMMENDATIONS 

1. Install the quartz sample probe and 1/2" teflon or nylon lines throughout the 
sampling system so a non-catalytic sampling system can be compared to the existing 
pitot tube sampling system. 

2. Verify the presence and percent concentrations of nitric oxide observed at higher 
currents spectrographically. 

3. Install pyrex orifices in place of valves and use inert valve bodies for locations 
where valves are necessary (solenoid valve and mass spect chamber valve). 
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Mol* Fraction 



Figure 3. Mole Fractions of Oxygen, Oxygen Atoms, Nitrogen Atoms, and Nitric 
Oxide as a Function of Pitot Position and Current. 
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Figure 4. Mass Spectrometer Signal Counts at m/e*30 (NO) for a Shock Survey 
around a Blunt Body. 
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SECTION 2 


CATALYTIC MEASUREMENTS ON HRSI AND RCC MATERIALS 

BACKGROUND 

This objective involved the determination of atomic recombination rates on HRSI 
(high temperature reuseable surface insulation) and RCC (reinforced carbon carbon) 
materials at various gas compositions, flow rates, and enthalpies. 

Knowledge of catalytic responses on thermal protection systems are of 
tremendous importance in heat shield design. The difference between a non catalytic 
surface and fully catalytic surface influences the rate of heat transfer to the surface 
during a reentry manuever, and therefore, the difference in temperature measured at 
the surface. The ideal surface would be non catalytic i.e. atoms would strike the 
surface and reflect off. On the other hand, on a fully catalytic surface all atoms 
which strike the surface recombine to the ground state, and the energy released in 
this process would be transferred to the surface. Thus the surface can heat up 
tremendously. Demonstration of the catalytic effect on surfaces was reported for 
STS- 2 which had 2 catalytic points on its thermal protection system (Stewart, 1931) 
and subsequent analysis predicted surface temperature increases of 270 K at these 
catalytic points (Scott, 1982). 

Data analysis followed an approached developed by Scott (1981,83). The analysis 
begins by using simple flow conservation equations for an inviscid stream tube and 
frozen flow relations across the shock wave (Scott, 1983). Shock properties are 
determined using the measured quantities of enthalpy, stagnation pressure, and mass 
flowrate along with the calculated effective area taking into account the boundary- 
layer effect. The effective area was determined using the NATA code (Bade and Yo$, 
1976). Then the Fay and Riddell (1958) relation as modified by Goulard (1958) and 
further refined by Scott (1981) was used to generated curves for the ratio of heating 
rate to fully catalytic heating rate versus recombination coefficient. Since 
recombination can occur for both nitrogen and oxygen atoms the approach requires 
separating the influence of nitrogen before analyzing results for air. This is done by 
running tests in pure nitrogen conditions and obtaining the nitrogen recombination 
coefficient in terms an Arrhenius expression which would serve as an input into the 
air analysis. Further details about the method can be found in Scott (1981, 83). 


EXPERIMENTAL 

Hemispherical test articles made out of HRSI (LI-900) RCG (reaction coated 
glass) were place in the 5-MW arc jet chamber at JSC. The arc column consisted of 
5 heater packs attached to a 2.25 inch throat which was connected to a 15° solid 
angle nozzle with an exit diameter of 15 inches. Further details about the facility 
made be found in Rochelle et al. (1983). The test articles were placed in the right 
arm while a comparable metallic hemispherical model instrumented with slug 
calorimeters was placed in the left. The slug calorimeter served as a reference 
heating rate measurement. The heating rate was determined by measuring the rise in 
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temperature divided by the change in time during the period the calorimeter was 
inserted into the stream. The constant for the calorimeter was 0.930 BTU/ft 2 /°F The 
recombination coefficients for the calorimeter were assumed constant at 0.6 for 
nitrogen and 0.15 for oxygen. 

The procedure was: 

1. Arc on. 

2. Run for 2 to 5 minute period to allow gas flows and arc conditions to stabilize. 

3. Swing in the left arm (calorimeter) into the stream for a period of 1 to 2 seconds. 

4. Swing in the right arm (the model) into the stream wait until the stagnation point 
thermocouple reads constant, take a test point. 

5. Swing in the pitot tube into the stream, take a test point to obtain the stagnation 
pressure. 

6. Change to the next current setting. 

7. Repeat steps 3 to 6 until full ranges have been covered. 

The range of conditions studied were: 


Nitrogen 

0.038 kg/sec (0.08 lb/sec) 

at 400 to 550 amps 

11% (mass) O 2 in N 2 

0.038 kg/sec (0.08 lb/sec) 

at 300 to 400 amps 

14% O 2 in N 2 

0.056 kg/sec (0.14 lb/sec) 

at 300 to 340 amps 

23% O 2 in N 2 (Air) 

0.038 kg/sec (0.08 lb/sec) 

at 350 to 450 amps 


The sequence of testing was: 

7/12/88 Air - 12 points followed by Nitrogen - 1 point 

7/14/88 Nitrogen - 7 points followed by 14% 02-3 points 

7/15/88 Air - 3 points followed by 11% O 2 -6 points followed by N 2 - 3 points 

The experimental procedure included several repeat points in order to check 
repeatability of the results. 


26-9 



RESULTS AND DISCUSSION 


Table 1 presents conditions tested and resultant heating rates and temperatures 
measured. Figure 5 shows that good repeatability could be achieved within a run, 
however, between runs temperatures for air conditions varied by about 60 K. Another 
observation made is that in the presence of 11% oxygen, the model temperatures were 
hotter than the Air-1 case. Heating rates also followed a similar trend. Figure 6 
shows the heat flux normalized by the square root of the static pressure. In this 
figure one sees that the trends are uniform, however, an unexplained difficultly is . 
why the measured normalized heating rate was about 2 times the value predicted by 
the Fay-Riddell correlation (shown as the straight line). The Fay-Riddell prediction 
was also checked through the NATA code which showed it value to be slightly higher 
than calculated by the catalysis program primarily because of higher predicted shock 
properties of temperature and pressure. It should also be noted that the NATA code 
calculated a Lewis number of 1 .27 while the catalysis code calculated a Lewis number 
of less than 1 at around 0.9 (dependent on the experimental point). Lewis numbers 
less than 1 result when the total atomic mass fraction is high. This signifies that the 
ratio of diffusion transport to the surface is less then thermal transport away from 
the surface. 

Figure 7 shows the recombination coefficient determined for nitrogen. The curve 
is relatively flat representing a process that has low activation energy (gas phase 
recombination reactions also have low activation energies). The magnitude of the 
results compare favorably with those reported by Scott (1981) (his values ranged from 
0.014 to 0.02 in the similar temperature range). The results also compare with 
Kolodziej and Stewart (1987) in which they reported values from 0.008 to 0.016. Both 
papers reported relatively low activation energies which also compare with these 
results. 

Figure 8 shows the recombination coefficient determined for oxygen. These 
curves show the difficultly of data repeatability between runs. Further, one would 
expect the results for the 14% O 2 to fall between Air-2 and the 11% O 2 if the Air-1 
data could be thrown out for an obvious reason. Attempts to explain the day to day 
discrepancies included verification of voltage measurement, current measurement and 
water flow rates. The voltage measurement check involved connecting a Simpson 
meter across the supply lines. The Simpson meter measurement checked with the 
value measured by the Norbitrol. Current calibrations were done several times with 
no major discrepancies found. Finally, the water turbine meter position was changed 
so that no sudden flow disturbance occurred upstream. No flow change was observed 
for this correction. In a review of arc jet data, the only major discrepancies was 
with the oxygen supply pressure and actual flowrate. However, the product of these 
two were the same for both air runs as it should be. Further, the column pressures 
were very similar indicating that the arc jet column conditions were identical. So 
what explanations remain to be explored? Surface change is the most likely 
candidate. Perhaps emissivity changes or perhaps nitrogen changes. Another 
explanation may be thermocouple relocation, however, two independent nitrogen runs 
were repeatable. Another possibility may be the model location was not repeatable 
and this should be explored further. In summary, further testing is required. 
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TABLE 1 


SUMMARY OF CONDITIONS STUDIED FOR CATALYTIC RECOMBINATION 
ON RCG COATED HRSI MATERIAL 


Run 

Number 


2-749-1 

2-749-2 

2-749-3 

2-749-4 

2-749-5 

2-749-6 

2-749-7 

2-749-8 

2-749-9 

2-749-10 

2-749-11 

2-751-1 

2-751-2 

2-751-3 

2-751-4 

2-751-5 

2-751-6 

2-751-7 

2-751-8 

2-751-9 

2-751-10 

2-753-1 

2-753-2 

2-753-3 

2-753-4 

2-753-5 

2-753-6 

2-753-7 

2-753-8 

2-753-9 

2-753-11 

2-753-12 

2-753-13 


Gases 

Air 


Air 


Air 


Air 


Air 


Air 


Air 


Air 


Air 


Air 


n 2 


n 2 


n 2 


n 2 


n 2 


n 2 


n 2 


n 2 


14% 

o 2 

14% 

o 2 

14% 

o 2 

Air 


Air 


Air 


11% 

02 

11% 

o 2 

11% 

o 2 

11% 

o 2 

11% 

o 2 

11% 

o 2 

n 2 


n 2 


n 2 



Flowrate 

kg/sec 


0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0368 

0.0363 

0.0381 

0.0381 

0.0381 

0.0381 

0.0381 

0.0390 

0.0368 

0.0558 

0.0563 

0.0563 

0.0359 

0.0359 

0.0359 

0.0368 

0.0368 

0.0368 

0.0368 

0.0395 

0.0372 

0.0368 

0.0368 

0.0368 


Current 

amps 


350 
376 
403 
426 

450 

451 
425 

401 
376 
349 

402 
400 
451 
501 
524 
550 

500 

400 
301 
326 
341 

351 
374 

401 
300 
326 
351 
376 
400 
326 
400 

501 
550 


Enthalpy 

MJ/kg 


7.31 
7.75 
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Figure 6. Normalized Heat Flux Measured as a Function of Total Enthalpy. 
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Figure 7. Nitrogen Recombination Coefficient as a Function of Temperature - 
Arrhenius Plot. 



0.64 0.66 0.66 0.67 0.68 0.69 0.7 

1000/Tw, 1/K 

Figure 8. Oxygen Recombination Coefficient as a Function of Temperature - for 
Various Run Conditions. 
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After the test program was completed, the surface of the slug calorimeter was 
analyzed by energy dispersive X-Ray spectroscopy (EDAX). The results showed that 
the silver coating agglomerated on the surface exposing the copper surface 
underneath. Further elemental analysis showed that the small clumps of were 
primarily composed of silver with a trace of oxygen and copper. The valley areas 
were primarily composed of copper with some oxygen, however, the oxygen level was 
not high enough to represent cuprous oxide or cupric oxide on the surface. 


CONCLUSIONS 

1. HRSI - RCG coated thermal protection materials are relatively inert to nitrogen 
atom recombination. 

2. Repeatability within runs was possible, however, outside of runs, results were not 
as repeatable. 

3. The scatter of data observed between runs for oxygen recombination experiments 
cannot be attributed to differences in voltage, current, or water flowrate measurement 
errors. 

3. The silver coated calorimeters changed surface structure morphology during the 
test program with silver agglomerates forming on the surface. 


RECOMMENDATIONS 

1. Continue experimental measurements. Temperature ranges between 1600 and 1800 
K should be investigated. 

2. Attempt catalytic measurements on RCC (Reinforced Carbon Carbon) materials. 

3. Determine why the surface temperature differed by about 60 K for identical arc 
heater conditions between the two air runs. 

4. Include emissivity measurements in all catalysis test programs. 

5. Insure that all arc jet measurements are reliable before testing. 
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